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SECTION  I 

INTRODUCTION  AND  SUMMARY 


MODELING  AND  IDENTIFICATION  REQUIREMENTS 

The  regulation  accuracy  and  speed  of  response  of  any  control  system  is 
directly  proportional  to  its  bandwidth.  The  increasing  demand  for  reliable 
engine  controls  with  tighter  surge  margins  requires  higher-bandwidth  control 
systems.  Thus,  there  is  a need  to  model  engine  dynamics  out  to  higher 
frequencies  than  has  been  required  in  the  past.  It  is  not  sufficient  to  provide 
large  gain  and  phase  margins  to  provide  for  model  uncertainty.  These 
margins  reduce  the  bandwidth  and  therefore  the  regulation  accuracy.  Some 
performance  will  have  to  be  sacrificed  to  allow  for  engine  variations  due  to 
tolerances  and  wear.  However,  performance  should  not  be  sacrificed  for 
inaccurate  modeling  practice. 

This  report  presents  a modeling  and  identification  procedure  which  is  simple 
and  inexpensive.  In  addition,  it  provides  both  high-fidelity  models  and 
dynamic  state  identification  which  gives  an  indication  of  the  physical  phe- 
nomena which  are  occurring.  The  procedure  has  been  used  for  many  years. 
The  accuracy  of  the  technique  has  been  considerably  improved  by  the  use  of 
two-channel  Fourier  filtering  and  computerized  state  identification  of  engine 
frequency  responses. 

The  optimal  control  system,  which  is  described  in  Volumes  I and  II  of  this 
report,  demonstrated  improved  regulation  of  the  engine  at  steady  state  and 
during  transients.  This  was  accomplished  by  linear  regulation  of  pressure 
and  temperature  boundaries  as  well  as  rotor  speed.  Each  of  these  regulators 
used  fuel  control  bandwidths  higher  than  standard  practice.  In  implementa- 
tion of  the  optimal  control  system,  the  bandwidth  had  to  be  reduced  to  allow 
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margin  for  differences  between  the  model  and  the  engine.  This  provided  the 
incentive  to  study  the  engine  dynamics  at  slightly  higher  frequencies. 

It  is  important  to  have  a practical  technique  for  accurately  determining  engine 
dynamics.  If  nominal  engine  dynamics  were  accurately  known,  as  well  as 
engine- to- engine  variations  in  the  dynamics,  the  proper  amount  of  gain  margin 
could  be  allowed  for  the  variations.  Then  excess  performance  would  not  be 
sacrificed  because  of  lack  of  knowledge  of  the  nominal  dynamics.  These  con- 
siderations will  become  more  important  when  control  systems  are  designed 
for  disturbance  insensitivity  in  the  future.  Disturbances  such  as  augmentor 
lightoff  and  inlet  instability  will  require  higher  bandwidth  control  systems. 

The  primary  objective  of  this  modeling  program  was  to  develop  a modeling 
procedure  which  would  be  very  practical.  The  procedure  outlined  in  this 
report  meets  that  objective.  The  procedure  is  inexpensive  due  to  the  small 
amount  of  engine  running  time  required.  Only  0.  5 hour  is  required  to  oscil- 
late each  engine  actuator  over  the  frequency  range.  The  nine  engine  responses 
were  recorded  on  tape  simultaneously.  This  was  possible  because  of  the  tape 
synchronizing  signals  used  by  the  servoanalyzer.  The  oscillations  are  very 
small,  on  the  order  of  2 percent  of  actuator  deflection,  so  that  the  engine  is 
not  endangered  during  the  test. 

The  accuracy  of  the  results  was  improved  in  three  ways.  First,  the  actuator 
dynamics  were  separated  from  the  engine  dynamics  with  the  two- channel 
capability  of  the  analyzer.  The  Fourier  filtering  capability  allowed  analysis 
at  low  signal-to-noise  levels  and  provided  describing  functions  for  nonlinear 
dynamics  such  as  dead  zone  and  hysteresis.  Finally,  the  automation  of  the 
identification  procedure  provided  very  close  matching  of  the  models  to  the 
frequency  responses. 

The  limitation  of  the  procedure  is  that  it  is  a linear  analysis  of  a nonlinear 
process.  Several  operating  points  must  be  evaluated  to  determine  the  varia- 
tions in  the  linear  models.  If  the  nonlinearities  are  dominant  even  in  the 
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small  perturbations,  such  as  with  hysteresis  and  saturation,  their  describing 
functions  for  sinusoidal  inputs  are  not  always  useful.  However,  linear  models 
are  required  for  all  types  of  control  synthesis.  Therefore,  this  procedure 
should  be  very  useful  even  though  nonlinear  models  are  the  ultimate  objective. 

Experimental  data  of  engine  dynamics  are  very  limited  for  two  reasons. 

First,  engine  test  time  is  very  expensive,  so  that  only  the  highest  priority 
- tests  are  performed.  Second,  the  dynamics  are  well  known  at  low  frequency, 
within  the  bandwidth  of  current  control  systems.  Higher-bandwidth  systems 
will  increase  the  priority  of  dynamic  tests  in  the  future,  at  least  until  the 
dynamics  become  well  known  at  slightly  higher  frequencies.  The  technique 
described  in  this  report  should  prove  to  be  very  cost  effective,  since  the 
engine  test  time  can  be  reduced  to  a few  hours.  The  primary  impact  this 
technique  has  on  engine  testing  is  the  requirement  for  electrical  readout  on 
actuator  positions  and  engine  sensors. 

The  two- channel  Fourier  analyzer  is  a key  item  to  making  the  procedure 
practical.  This  is  because  it  can  analyze  the  amplitude  ratio  and  phase  shift 
between  any  two  signals  in  a very  noisy  environment.  The  computerized 
identification  procedure,  which  determines  the  dynamic  states  and  transfer 
functions  from  the  frequency  responses,  is  very  useful  for  determining 
practical  models  of  the  engine  dynamics.  These  models  provide  insight  to 
the  physics  of  the  engine  which  can  be  used  as  guidelines  in  future  analytical 
modeling  efforts. 

Frequency  response  analysis  is  limited  to  the  study  of  small  perturbations 
around  a nominal  operating  point.  This  linear  analysis  is  required  for  con- 
trol system  design,  but  is  not  a substitute  for  the  nonlinear  modeling  of  the 
engine.  The  frequency  response  analysis  is  also  limited  to  perturbations 
around  steady- state  operating  points,  since  the  engine  must  remain  at  the 
nominal  operating  point  for  a few  minutes  at  a time.  The  modeling  of 
transient  conditions  must  be  accomplished  by  approximating  transient  loads 
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on  the  rotor.  However,  for  control  design  purposes,  frequency  response 
measurements  can  be  obtained  from  a hybrid  simulation  at  any  operating 
condition  and  compared  with  the  engine  responses  at  steady  state  to  verify 
the  results. 


FORM  OF  THE  MODEL 

A dynamic  model  was  obtained  for  the  response  of  five  engine  variables  to 
the  four  engine  actuators.  These  20  linear  transfer  functions  were  obtained 
at  three  operating  points.  The  actuator  transfer  functions  were  separated 
from  the  engine  dynamics.  Thus,  there  were  64  transfer  functions  evaluated 
in  all.  Figures  1,  2,  and  3*  show  these  transfer  functions  combined  into  a 
dynamic  transfer  matrix  for  each  operating  point.  It  is  clear  from  the  matrix 
that  there  is  considerable  interaction  between  the  control  variables,  since 
each  actuator  has  an  effect  on  most  of  the  responses.  However,  the  effect  of 
the  inlet  guide  vane  (IGV)  and  bleed  (BLD)  variables  is  considerably  less  than 
the  fuel  flow  (Wf)  and  exhaust  area  (Ag).  This  fact  is  obscured  in  the  transfer 
function  by  the  units  chosen  for  the  DC  gain.  Fuel  flow  is  expressed  in  pounds 
per  hour  and  exhaust  area  in  square  inches,  while  IGV  and  BLD  are  expressed 
in  percent  of  full  scale  (i.e. , Pg/BLD  is  psi /full- scale  deflection). 

A state  space  model  of  the  engine  can  be  obtained  from  the  transfer  matrix. 
The  dynamic  states  are  the  poles  of  the  transfer  functions.  The  order  of  the 
state  space  model  is  the  lowest- order  denominator  that  can  be  factored  out  of 
the  matrix  (i.  e. , add  up  all  the  unique  poles.  If  they  appear  more  than  once, 
they  are  counted  only  once).  Before  one  could  determine  the  lowest- order 
state  space  model,  the  roots  would  have  to  be  analyzed  to  determine  which 
ones  are  appearing  in  several  responses.  It  is  not  alwayr  clear,  as  the 


*To  avoid  interrupting  the  continuity  of  the  text,  all  referenced  figures  and 
tables  are  gathered  at  the  end  of  their  respective  section  or  appendix. 
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errors  in  the  procedure  cause  states  to  appear  at  slightly  different  frequen- 
cies in  different  responses.  Not  all  of  the  states  have  to  be  included  for 
control  design  purposes.  High-frequency  terms  can  be  left  off.  After 
truncating  the  terms,  frequency  response  plots  can  be  compared  with  the 
experimental  data  to  ensure  that  accuracy  is  maintained  out  to  the  desired 
bandwidth  of  the  control  system. 


RESULTS  OF  ENGINE  RESPONSE  MEASUREMENTS 

Engine  frequency  responses  reveal  much  information  about  engine  dynamics 
in  addition  to  that  needed  for  modeling  purposes.  In  the  paragraphs  that 
follow,  some  of  the  interesting  facts  picked  out  of  the  data  are  discussed. 
Much  of  this  information  would  not  be  available  were  it  not  for  the  Fourier 
filtering  capability  of  the  instrumentation.  This  noise-rejection  capability 
allowed  accurate  responses  to  be  analyzed  out  in  frequency  to  the  limit  of 
the  actuator  response.  The  fuel  flow  responses  were  measured  out  to  100  Hz. 
The  geometry  actuators  proved  to  be  much  faster  responding  than  originally 
modeled. 

It  should  be  noted  that  fuel  flow  dynamics  cannot  be  separated  from  engine 
dynamics  as  desired.  This  is  due  to  the  fact  that  fuel  flow  is  effected  by 
pressure  drop  across  the  fuel  nozzle.  The  dynamics  of  the  combustor  pres- 
sure appear  in  both  the  actuator  and  engine  responses.  This  phenomenon  is 
described  in  detail  in  Appendix  A.  This  suggests  a need  for  flow  feedback 
in  engine  control  systems. 

The  exhaust  area  actuator  was  much  faster  than  expected.  It  responsed  well 
out  to  8 Hz.  The  hysteresis  and  dead  zones  added  only  30  degrees  of  phase 
shift  at  low  frequency.  Of  course,  the  nonlinearity  in  the  actuator  will  cause 
it  to  respond  differently  at  other  amplitude  levels.  The  amplitude  level  for 
This  data  was  ±5  percent  of  full  scale.  This  control  variable  has  a large 
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effect  on  the  engine  response.  Therefore,  it  should  prove  to  be  an  important 
control  in  a multivariable  system.  There  is  considerably  more  phase  shift 
in  the  combustor  pressure  response  to  exhaust  area  than  there  is  to  fuel  flow, 
however.  If  the  optimal  control  system  were  to  be  redesigned,  the  exhaust 
area  would  be  added  as  a second  dynamic  control  variable. 

As  mentioned  above,  the  compressor  bleeds  and  inlet  guide  vanes  had  little 
effect  on  the  engine  responses.  As  noted  in  the  transfer  matrix,  some  of  the 
transfer  functions  are  shown  as  zero.  The  small  effect  is  obscured  in  the 
transfer  functions  because  of  the  units  chosen  for  these  variables.  The  gain 
in  the  transfer  functions  appears  high  because  it  is  in  terms  of  output/full- 
scale  deflection.  If  fuel  flow  and  exhaust  area  were  to  be  expressed  in  these 
units,  the  gain  would  be  much  higher.  The  combustor  pressure,  Pg,  response 
to  bleeds  has  very  little  phase  shift,  which  means  that  the  effect  is  immediate. 
However,  the  magnitude  of  the  effect  is  small,  presumably  due  to  the  small 
size  of  the  bleed  openings.  Thus,  the  bleeds  and  inlet  guide  vanes  are  not 
very  effective  dynamic  control  variables  for  this  engine. 

One  of  the  most  interesting  results  of  the  frequency  response  tests  is  the 
second-order  dynamic  response  of  spool  speed.  This  second  state  appearing 
in  the  spool  speed  response  causes  considerable  phase  shift  at  frequencies 
within  the  bandwidth  range  of  most  current  engine  control  systems.  The  fuel 
flow  response  shows  the  two  first-order  lags.  The  first  appears  as  expected 
at  3 radians  per  second  and  the  second  at  77  radians  per  second  (at  the  high- 
speed operating  point).  The  exhaust  area  response  shows  a second  lag  in 
spool  speed  also,  but  it  moves  out  in  frequency  at  low  speed.  The  additional 
phase  shift  in  spool  speed  response  caused  considerable  difficulty  when  the 
high-bandwidth  optimal  controllers  were  run  on  the  engine.  Therefore,  the 
engine  tests  were  run  with  reduced  gains  to  allow  for  additional  gain  and 
phase  margin. 
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Another  very  interesting  result  revealed  in  the  frequency  responses  is  the 
unexpectedly  large  time  delay  in  the  engine  pressure  responses.  This  time 
delay  ranges  from  11  to  14  milliseconds.  It  causes  very  rapid  phase  shift  at 
frequencies  above  10  Hz.  This  is  above  the  frequency  range  of  speed  control 
bystems.  However,  pressure  disturbance  control  systems  in  the  future  will 
have  to  consider  this  time  delay. 

The  turbine  outlet  temperature,  T5,  response  provided  accurate  dynamic 
data  out  to  50  Hz,  even  though  the  thermocouple  has  very  slow  response. 

The  added  gain  reduction  and  phase  shift  did  not  prevent  temperature 
response  models  from  being  obtained.  This  implies  that  fast  temperature 
control  is  possible  with  adequate  lead  compensation  for  the  thermocouple. 
However,  the  signal  noise,  which  was  rejected  by  Fourier  filtering,  will  limit 
the  amount  of  lead  that  can  be  applied. 


COMPARISON  OF  EXPERIMENTAL  AND  ANALYTICAL  MODELS 

One  of  the  main  objectives  of  this  program  was  to  compare  the  engine  model, 
used  for  optimal  control  design,  with  experimental  data.  Engine  frequency 
response  data  were  plotted  with  frequency  responses  from  the  linearized 
NASA  component  model  and  the  analog  computer  model  at  APL.  These  data 
are  presented  in  Section  II.  The  following  paragraphs  point  out  some  of  the 
interesting  comparisons. 

The  fuel  flow  responses  are  of  primary  concern,  as  fuel  flow  is  usually  the 
only  dynamic  control  variable.  The  fuel-metering  actuator  is  modeled  accu- 
rately. The  spool  speed  response  is  modeled  accurately  at  low  frequency. 

As  mentioned  above,  there  is  considerably  more  phase  shift  in  the  engine 
data  above  2 Hz.  The  experimental  model  has  a second  first-order  lag  at 
about  8 Hz.  This  error  in  the  model  meant  that  the  bandwidth  of  the  optimal 
controller  had  to  be  reduced  during  the  engine  tests.  This  caused  a corres- 
ponding decrease  in  regulation  accuracy. 
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The  pressure  responses  are  accurate  in  the  component  model  except  for  the 
large  time  delay.  The  analog  model  has  considerable  phase  error  in  the 
AP/P  and  P5  responses  at  low  frequency.  The  temperature  responses  have 
considerable  phase  error  above  7 Hz.  In  addition,  the  analog  model  tempera- 
ture gain  begins  falling  off  rapidly  at  0.  3 Hz,  which  is  about  a decade  slower 
than  the  engine  data. 

The  exhaust  area  actuator  is  a decade  faster  than  the  models.  The  gain  is 
flat  out  to  5 Hz,  as  compared  with  0.5  Hz.  This  difference  is  very  important 
in  deciding  whether  to  use  the  exhaust  area  a:  a dynamic  control  variable. 

The  optimal  controller  did  not  use  the  exhaust  area  because  it  was  assumed 
to  be  too  slow.  At  5 Hz,  the  exhaust  area  can  be  used  effectively.  If  the 
optimal  controller  were  redesigned,  the  exhaust  area  would  be  added  as  a 
second  control  variable.  This  would  have  a considerable  effect  on  the  con- 
troller performance. 

The  spool  speed  response  to  exhaust  area  has  similar  second-order  dynamics 
as  with  fuel  flow.  The  phase  shift  increases  from  that  of  the  model  above 
3 Hz.  The  combustor  pressure  responses  to  exhaust  area  are  well  modeled. 
However,  the  turbine  discharge  pressure,  P&,  varies  considerably  between 
the  three  models.  responds  faster  and  with  greater  amplitude  on  the 
engine  than  in  the  component  model.  The  analog  model  is  1.5  decades  slower. 

The  bleed  and  IGV  actuators  are  much  faster  than  the  models.  However,  this 
is  not  an  important  factor,  since  these  actuators  are  not  very  effective, 
dynamically.  The  models  show  these  control  variables  to  be  even  less 
effective.  Spool  speed  response  to  bleed  has  more  phase  shift  than  the  model 
above  0. 1 Hz.  The  response  to  bleed  is  1.5  decades  slower  in  the  models. 
The  data  shows  very  little  phase  shift  in  pressure  response. 
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STATE  IDENTIFICATION  AND  PHYSICAL  INTERPRETATION 

Dynamic  engine  models  were  identified  from  the  engine  frequency  response 
data.  This  procedure  proved  to  be  very  successful  in  that  accurate  pole  and 
zero  locations  were  identified  that  coincided  well  with  engine  physics.  This 
is  an  established  practice  in  the  control  field.  However,  the  technique  was 
improved  by  the  accurate  measurement  techniques  and  computerized  statis- 
tical identification.  Therefore,  the  procedure  promises  to  be  a practical 
technique  for  future  use.  The  identification  algorithm  is  discussed  in  detail 
in  Appendix  B.  The  models  are  discussed  in  detail  in  Section  III.  The 
transfer  functions  are  listed  in  Figures  1,  2,  and  3,  with  some  of  them 
shortened  for  clarity. 

An  attempt  was  made  to  interpret  the  physical  meaning  of  the  dynamic  states 
in  Section  III.  Many  of  the  interpretations  are  well  known  from  comparison 
with  existing  models.  But,  the  higher-frequency  terms  do  not  have  obvious 
interpretations.  The  authors  have  made  an  attempt  to  classify  these  addi- 
tional dynamics.  However,  this  task  is  more  appropriately  accomplished  by 
engine  manufacturers  who  have  more  experience  with  engine  physics.  An 
objective  of  this  program  was  to  provide  data  to  allow  engine  modelers  to 
perfect  their  analytical  procedures. 

The  significant  differences  between  the  engine  and  the  component  model  are: 
(1)  the  additional  first-order  lag  in  the  spool  dynamics  at  8 Hz;  (2)  the  0.  011- 
second  time  delay  in  the  pressure  response;  and  (3)  the  location  of  the  second- 
order  pole  in  the  exhaust  actuator.  The  additional  spool  root  could  be  asso- 
ciated with  gas  dynamics  in  the  turbine  or  compressor.  The  time  delay  is 
assumed  to  be  occurring  in  the  fuel  combustion.  The  dynamics  of  the  exhaust 
actuator  are  understood  fairly  well. 
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Figure  2.  J85- 13  Dynamic  Trar 
at  85  Percent  Maximi 
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SECTION  n 

FREQUENCY  RESPONSE  DATA -COMPARISON  OF  ENGINE, 

ANALOG  MODEL,  AND  LINEARIZED  NASA  COMPONENT  MODEL 

INTRODUCTION 

Frequency  response  data  experimentally  obtained  from  a J85  engine  in  a 
test  cell  at  APL  are  presented  in  this  section  and  compared  with  similar  data 
obtained  from  two  analytic  models,  the  linearized  NASA  component  model 
and  the  APL  analog  model.  The  Bode  plots  presented  show  that  both  analytic 
models  adequately  represent  the  low-frequency  dynamics  of  the  J85  engine 
below  2 Hz.  However,  significant  high-frequency  effects  above  2 Hz  asso- 
ciated with  gas  dynamics  identified  from  the  engine  responses  are  not  in- 
cluded in  either  analytic  model.  Details  of  the  engine-analytic  model  com- 
parison are  discussed  in  the  following  paragraphs. 

A BAFCO  servoanalyzer  was  used  to  obtain  the  J85  frequency  response  plots. 
This  instrument  is  designed  to  perform  frequency  response  analysis  of  com- 
plicated servomechanisms  such  as  the  J85  engine  through  the  implementation 
of  Fourier  analysis.  Since  it  is  a two- channel  analyzer,  it  can  measure  the 
dynamics  between  any  two  outputs;  hence,  it  is  not  restricted  to  input-output 
pairs.  This  feature  allows  separation  of  actuator  dynamics  from  engine 
dynamics.  Operation  of  the  analyzer  is  shown  in  diagram  form  in  Figure  4 
and  briefly  discussed  here.  A complete  description  of  the  analyzer  is  in- 
cluded in  Appendix  A. 

Th^analyzer  produces  a sinusoidal  voltage  signal  with  time-dependent  fre- 
quency (i.  e. , the  frequency  varies  logarithmically  with  time)  which  is  used 
to  drive  one  of  the  engine  actuators.  Responses  from  two  engine  sensors 
are  fed  back  into  the  analyzer  which  contains  the  necessary  electronics  to 
compute  the  amplitude  ratio  and  phase  shift  between  the  two  sensor  signals. 
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This  information  is  recorded  as  a function  of  frequency  by  two  x - y plotters, 
thereby  producing  the  Bode  frequency  response  plot  of  sensor  2 with  respect 
to  sensor  1.  These  input  and  output  variables  are  described  schematically 
in  Figure  4. 


The  BAFCO  servoanalyzer  was  used  in  this  manner  to  obtain  frequency 
responses  of  various  engine  parameters  with  respect  to  the  four  engine  con- 
trols: fuel  flow,  exhaust  area,  compressor  bleed,  and  inlet  guide  vane.  Four 
frequency  sweeps,  one  for  each  control  variable,  were  performed  at  each  of 
three  operating  points:  engine  speed  (N/N  ) = 70  percent,  85  percent,  and 
95  percent.  Representative  data  obtained  from  these  tests  are  presented  in 
Figures  5 through  29.  Complete  documentation  of  the  data  is  included  in 
References  1 through  3.  The  frequency  responses  included  in  these  figures 
are  discussed  in  the  subsections  that  follow,  in  the  following  order: 


Actuator 

responses 


Engine  responses 
to  fuel  flow 


Ref. 

Figure 

Response 

N/N 

max 

5 

Fuel  flow/ fuel  command 

95% 

6 

Exhaust  area  /exhaust  command 

95% 

7 

Compressor  bleed/bleed 
command 

95% 

8 

Inlet  guide  vane /inlet  guide 
vane  command 

95% 

9 

Spool  speed /fuel  flow 

95% 

10 

Spool  speed/fuel  flow 

85% 

11 

Spool  speed /fuel  flow 

70% 

12 

Compressor  discharge  pressure/ 
fuel  flow 

95% 

13 

Compressor  discharge  pressure/ 
fuel  flow 

85% 

14 

Compressor  discharge  pressure/ 
fuel  flow 

70% 

15 

Turbine  discharge  pres  sure /fuel 
flow 

95% 

18 


Eef. 

N/N 

Figure 

Response 

IT 

Engine  responses 
to  fuel  flow 

16 

Turbine  discharge  pressure/ 
fuel  flow 

85% 

(continued) 

17 

Turbine  discharge  pressure/ 
fuel  flow 

70% 

18 

Mach  number  sensor/fuel  flow 

95% 

19 

Turbine  discharge  temperature/ 
fuel  flow 

95% 

Engine  responses 

20 

Spool  speed/exhaust  area 

95% 

to  exhaust  area 

21 

Compressor  discharge  pressure/ 
exhaust  area 

95% 

22 

Turbine  discharge  pressure/ 
exhaust  area 

95% 

23 

Turbine  discharge  temperature/ 
exhaust  area 

95% 

Engine  responses 

24 

Spool  speed/ compressor  bleed 

95% 

to  compressor  bleed 

25 

Compressor  discharge  pressure/ 
compressor  bleed 

95% 

26 

Turbine  discharge  temperature/ 
compressor  bleed 

95% 

Engine  responses  to 

27 

Spool  speed /inlet  guide  vane 

95% 

inlet  guide  vane 

28 

Compressor  discharge  pressure/ 
inlet  guide  vane 

95% 

29 

Turbine  discharge  temperature/ 
inlet  guide  vane 

95% 

Actuator  inputs  and  sensor  measurements  represented  in  these  data  are 
identified  below.  A more  complete  description  of  the  actuators  and  sensors 
is  included  in  Appendix  A. 
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Actuator  Input  Description 


1) 

2) 

3) 

4) 


Fuel  flow  command,  u^  - Request  to  fuel  valve. 

Exhaust  area  command,  u^  = Request  to  exhaust  nozzle  actuator 

Compressor  bleed  command,  Ug^Q  = Request  to  compressor 
bleed  actuator. 

Inlet  guide  vane  command,  u^^  = Request  to  inlet  guide  vane 
actuator. 


Sensor  Output  Description 


1)  Fuel  flow,  Wj  = Fuel  flow  into  combustion  chamber.  This 
signal  was  recorded  as  the  pressure  differential  across  the 
fuel  nozzle,  i.  e. , fuel  nozzle  pressure  minus  compressor 
discharge  pressure,  and  corrected  to  actual  fuel  flow  in 

lb  /hr  with  a steady-state  calibration. 

2)  Exhaust  area,  Ag  = Effective  cross-sectional  area  of  exhaust 
nozzle.  This  signal  was  recorded  as  the  feedback  voltage 
(calibrated  in  inches  squared)  from  a mechanical  potenti- 
ometer positioned  on  the  nozzle  drive  mechanism. 

3)  Compressor  Bleed,  BLD  = Effective  area  of  compressor 
bleeds.  The  scale  is  nondimensionalized  in  the  sense  that 
1 corresponds  to  fully  open  bleeds  and  0 corresponds  to 
fully  closed  bleeds.  A potentiometer  located  on  the  actuator 
mechanism  was  used  to  record  this  signal. 

4)  Inlet  Guide  Vane,  IGV  = Incidence  angle  of  inlet  guide  vanes. 
The  nondimensional  scale  is  constructed  with  0.0 
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corresponding  to  the  high-speed  position  of  the  inlet  guide 
vanes  and  1.  0 corresponding  to  low-speed  position.  A 
potentiometer  which  measures  actuator  movement  was  used 
to  record  this  signal. 

5)  Spool  Speed,  N = Angular  frequency  of  rotor  shaft.  A sensor 
which  measures  elapsed  time  per  revolution  of  the  rotor  was 
used  to  obtain  this  signal.  The  sensor  does  not  contain  any 
dynamics  in  the  frequency  range  tested. 

6)  Compressor  Discharge  Pressure,  Pg  = Static  pressure  at  the 
compressor  discharge.  This  signal  was  measured  with  a 
static  pressure  tap  embedded  in  the  wall  of  the  engine  slightly 
behind  the  compressor  outlet  guide  vanes. 

7)  Turbine  Discharge  Pressure,  Pg  = Total  pressure  at  turbine 
discharge.  The  P&  sensor  consists  of  a system  of  five  total 
pressure  probes  spread  around  the  engine  and  in  back  of  the 
turbine  discharge.  A single  signal  is  obtained  by  averaging 
the  outputs  of  the  five  probes. 

AP 

8)  Mach  Number  Sensor,  p-  = Total  minus  static  pressure 
divided  by  total  pressure  at  compressor  discharge.  This 
signal  was  obtained  from  a special  sensor  built  by  Bendix. 

All  subtraction  and  division  necessary  to  obtain  the  AP/Pg 
signal  is  performed  in  the  sensor  which  is  located  behind  the 
compressor  discharge. 

9)  Turbine  Discharge  Temperature,  Tg  = Temperature  at  tur- 
bine discharge.  The  Tg  sensor  is  composed  of  19  individ- 
ual thermocouples  coupled  in  parallel.  The  thermocouples 
are  spaced  around  the  engine  a few  inches  behind  the  turbine 
discharge. 
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Also  included  in  Figures  5 through  29  are  frequency  response  measurements 
obtained  from  the  two  analytic  models,  the  linearized  NASA  component  model 
and  the  A PL  analog  model.  The  linearized  NASA  component  model  is  the 
analytic  model  which  was  used  to  synthesize  optimal  controllers  for  the 
engine.  State  variables  associated  with  the  model  are  identified  in  Figure  30. 
A complete  description  of  the  model  is  included  in  Reference  4.  The  APL 
analog  model  is  described  in  Reference  5. 

A digital  computer  program  was  used  to  obtain  frequency  response  data  from 
the  linearized  NASA  component  model.  The  program  computes  the  amplitude 
ratio  and  phase  shift  between  an  input-output  pair. 

The  BAFCO  servoanalyzer  was  used  to  obtain  frequency  response  measure- 
ments from  the  APL  analog  model.  Inlet  guide  vane  and  compressor  bleed 
data  are  not  presented  for  the  analog  model,  since  the  model  does  not  contain 
representations  of  these  two  controls. 

Steady-state  data  defining  the  three  operating  points  examined  for  this  pro- 
ject are  presented  in  Tables  1,  2,  and  3.  Corresponding  data  obtained  from 
the  linearized  NASA  component  model  and  the  APL  analog  model  are  also 
listed. 


ACTUATOR  RESPONSES 

Bode  frequency  response  plots  for  the  four  engine  actuators  (fuel  valve, 
exhaust  nozzle,  compressor  bleed,  and  inlet  guide  vane)  are  presented  in 
Figures  5 through  8.  Also  presented  in  these  figures  are  frequency 
responses  of  the  actuator  models  included  in  the  two  analytic  models,  the 
linearized  NASA  component  model  and  the  APL  analog  model.  Comparison 
of  the  results  supports  two  observations:  (1)  the  fuel  valve  actuator  is  accu- 
rately represented  in  the  analytic  models,  and  (2)  the  engine  geometry  actu- 
ators exhibit  higher  bandwidth  than  their  counterparts  in  the  analytic  models. 
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The  frequency  response  of  the  fuel  valve  actuator  at  95  percent  of  maximum 
spool  speed  is  presented  in  Figure  5.  Both  the  linearized  NASA  component 
model  and  the  APL  analog  model  fuel  valves  are  accurate  representations 
of  the  engine  fuel  valve  in  the  low- frequency  region,  i.  e. , the  frequency 
responses  agree  up  to  about  4.  0 Hz.  There  are  significant  differences 
above  4.  0 Hz.  These  include  experimental  fuel  valve  dynamics,  centered  at 
7.  0 Hz,  and  high-frequency  phase  rolloff.  These  results  are  also  valid  at 
the  other  two  operating  points,  70  percent  and  85  percent  of  maximum  spool 
speed. 

The  frequency  responses  of  the  engine  exhaust  area  actuator  and  the  analytic 
models  are  presented  in  Figure  6.  They  show  the  engine  actuator  to  be  of 
higher  bandwidth  and  to  have  a different  phase  response  than  the  two  analytic 
models.  The  principal  differences  are:  (1)  the  .mponent  model  gain  rolls- 
off  much  sooner  and  faster  than  the  engine  actuator  gain,  anti  (2)  the  engine 
actuator  has  -30  degrees  phase  shift  at  low  frequency,  whereas  the  compo- 
nent model  has  little  or  no  phase  shift.  A nonlinear  hysterisis  effect  in  the 
engine  exhaust  actuator  is  responsible  for  the  -30  degrees  of  phase  shift  at 
low  frequency.  Since  the  component  model  is  a linear  model,  this  nonlinear 
effect  is  not  duplicated  in  the  component  model  response. 

The  nonlinear  hysterisis  effect  of  the  engine  actuator  is  included  in  the 
analog  model,  but  the  effect  is  too  pronounced:  the  analog  actuator  model 
has  about  20  degrees  more  phase  shift  than  the  engine  actuator.  This  dis- 
crepancy could  be  minimized  by  reducing  the  deadband  uncertainty  in  the 
nozzle  actuator  simulation  in  the  analog  model. 

Frequency  responses  of  the  compressor  bleed  actuator  and  inlet  guide  vane 
actuator  are  contrasted  with  those  of  the  linearized  NASA  component  model 
in  Figures  7 and  8.  Representations  for  the  analog  model  are  not  included 
in  these  figures,  since  compressor  bleed  and  inlet  guide  vane  effects  are  not 
included  in  the  analog  model. 


23 


The  main  difference  between  the  frequency  responses  of  the  engine  com- 
pressor bleed  and  inlet  guide  vane  actuators  and  the  frequency  responses  of 
the  component  model  actuators  is  that  the  engine  actuators  are  of  consider- 
ably higher  bandwidth  than  the  component  model  actuators.  This  incompati- 
bility also  accounts  for  the  observed  differences  in  phase  response. 


ENGINE  RESPONSE  TO  FUEL  FLOW 

Frequency  response  plots  of  spool  speed,  N,  compressor  discharge  pres- 
sure, Pg,  turbine  discharge  pressure,  P&,  compressor  Mach  number, 
AP/Pg,  and  turbine  discharge  temperature,  Tg,  for  oscillations  in  fuel  flow 
are  presented  in  Figures  9 through  19.  Comparison  of  these  engine  re- 
sponses with  similar  responses  obtained  from  the  linearized  NASA  compo- 
nent model  and  the  A PL  analog  model  leads  to  the  following  conclusions: 

1)  The  spool  speed  and  turbine  discharge  temperature  frequency 
responses  of  the  engine  agree  very  well  with  the  responses  of 
the  two  analytic  models  up  to  2.  0 Hz.  This  is  above  the  nor- 
mal bandwidth  for  speed  control  loops  of  1.  0 Hz. 

2)  The  pressure  and  Mach  number  responses  of  the  engine  agree 
very  well  with  the  responses  of  the  analytic  models  up  to 

10.  0 Hz  except  for  the  turbine  discharge  pressure  response  of 
the  analog  model. 

3)  High-frequency  dynamics  present  in  the  engine  responses, 
primarily  caused  by  time  delay,  are  not  represented  in  the 
analytic  models. 

4)  The  phase  shift  of  the  turbine  discharge  pressure  response  of 
the  analog  model  does  not  agree  with  either  the  engine  or  com- 
ponent model  results. 
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These  conclusions  are  discussed  in  the  following  paragraphs. 


Spool  speed  frequency  responses  are  presented  for  three  operating  points, 

70  percent,  85  percent,  and  95  percent  maximum  spool  speed,  in  Figures  9, 
10,  and  11.  Examination  of  these  responses  shows  that  both  the  component 
model  and  the  analog  model  contain  good  approximations  to  the  fuel  flow 
effects  of  the  engine  in  the  low-frequency  region,  up  to  2.  0 Hz.  Above  2.  0 
Hz  the  phase  response  of  the  engine  differs  considerably  from  the  phase  re- 
sponse of  the  two  models:  the  engine  phase  shift  is  much  greater  than  that  of 
either  model.  This  behavior  indicates  that  the  engine  contains  significant 
high-frequency  spool  dynamics  which  are  not  included  in  either  analytic  model. 

Compressor  discharge  pressure,  Pg,  responses  at  the  three  operating  points, 
are  presented  in  Figures  12,  13,  and  14.  Agreement  between  the  engine  data 
and  the  responses  of  the  component  model  and  analog  model  is  exhibited  out 
to  about  10  Hz.  The  principal  differences  between  the  engine  data  and  the 
analytic  models  are:  (1)  the  engine  data  contains  a time  delay  of  about  10  to 
15  milliseconds  which  is  not  represented  in  the  analytic  models,  and  (2)  the 
frequency  responses  indicate  that  the  engine  contains  significant  dynamics  in 
the  5.  0 to  8.  0-Hz  frequency  range  which  are  not  included  in  the  analytic 
models. 

Turbine  discharge  pressure,  Pg,  frequency  responses  presented  in  Figures 
15,  16,  and  17,  substantiate  most  of  the  conclusions  drawn  from  the  com- 
pressor discharge  pressure  responses.  Agreement  between  the  engine  data 
and  the  component  model  is  observed  out  to  about  10  Hz.  These  engine  re- 
sponses also  show  the  time  delay  of  about  10  to  15  milleseconds  and  the  pres- 
sure dynamics  in  the  5.  0 to  8.  0-Hz  frequency  range  which  were  noted  in  the 
compressor  discharge  pressure  responses.  However,  one  significant  differ- 
ence between  the  Pg  responses  and  the  Pg  responses  should  be  noted:  the 
analog  simulation  of  Pg  agrees  very  well  with  engine  data,  but  the  analog 
simulation  of  Pg  does  not  agree  with  the  corresponding  engine  data.  The 
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phase  shift  response  of  the  analog  Pg  response  is  incorrect  throughout  the 
frequency  range.  Correction  of  this  incompatibility  would  enhance  the  use- 
fulness of  the  analog  model. 

The  engine  Mach  number  sensor,  AP/Pg,  frequency  response  at  95  percent 
maximum  spool  speed  is  presented  in  Figure  18  and  compared  with  the  fre- 
quency response  of  the  analog  model  simulation.  No  data  are  presented  for 
the  component  model  because  it  does  not  include  a model  of  the  sensor.  The 
principal  difference  between  the  responses,  shown  in  the  figure,  is  that  the 
analog  response  has  less  phase  shift  than  the  engine  response.  Similar 
results  were  obtained  at  the  other  two  operating  points,  70  percent  and  85  per- 
cent maximum  spool  speed. 

Turbine  discharge  temperature,  Tg,  frequency  responses  of  the  engine  and  the 
two  analytic  models  at  95  percent  maximum  spool  speed  are  presented  in 
Figure  19.  The  results  show  agreement  between  the  engine  response  and  the 
responses  of  the  two  models  out  to  about  2.  0 Hz.  Beyond  2.  0 Hz  the  phase 
shift  of  the  engine  response  is  much  greater  than  the  phase  shift  of  the  models, 
indicating  that  the  engine  contains  some  high-frequency  dynamics  which  are 
not  identified  in  the  models.  This  observation  substantiates  the  conclusions 
drawn  from  the  spool  speed  responses  discussed  previously. 


ENGINE  RESPONSE  TO  EXHAUST  AREA 

Bode  plots  of  spool  speed,  N,  compressor  discharge  pressure,  Pg,  turbine 
discharge  pressure,  Pg,  and  turbine  discharge  temperature,  Tg,  for  oscilla- 
tions in  exhaust  area,  Ag,  are  presented  in  Figures  20  through  23.  Compari- 
son of  the  engine  responses  shown  in  these  figures  with  corresponding  re- 
sponses obtained  from  the  linearized  NASA  component  model  and  the  A PL 
analog  model  suggests  the  following  conclusions: 
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1)  DC  gain  levels  do  not  agree  very  well  between  the  engine 
responses  and  the  two  analytic  models.  This  discrepancy  is 
simply  a calibration  problem. 

2)  Except  for  differences  in  DC  gain  levels,  the  spool  speed 
responses  and  turbine  discharge  temperature  responses  of 
the  two  analytic  models  agree  fairly  well  with  the  corre- 
sponding engine  responses  in  the  frequency  range  below 

2.  0 Hz. 

3)  The  compressor  discharge  pressure,  P^,  response  of  the 
engine  is  accurately  simulated  by  the  two  analytic  models. 

4)  Neither  the  analog  model  nor  the  component  model  accurately 
represents  the  engine  turbine  discharge  pressure,  P&, 
response. 

These  conclusions  are  discussed  in  the  following  paragraphs. 

The  spool  speed  response  of  the  engine  at  95  percent  maximum  spool  speed  is 
compared  with  the  corresponding  responses  of  the  two  analytic  models  in 
Figure  20.  Except  for  differences  in  DC  gain  level,  the  model  responses  are 
seen  to  agree  with  the  engine  frequency  response  in  the  frequency  range  below 
2.  0 Hz.  Beyond  2.  0 Hz  the  engine  phase  shift  drops  considerably  below  the 
phase  responses  of  the  two  models,  indicating  that  the  engine  contains  some 
high-frequency  dynamics  which  are  not  included  in  the  models. 

Most  of  the  mismatch  in  DC  gain  levels  can  be  attributed  to  the  exhaust  area 
calibration  incompatibility  between  the  engine  and  the  models.  The  two 
exhaust  nozzle  simulations  and  the  engine  exhaust  area  sensor  all  have  differ- 
ent nozzle  area  calibrations.  A method  for  correcting  the  nozzle  area  cali- 
brations of  the  two  models  to  agree  with  the  engine  calibration  was  not  identi- 
fied because  engine  nozzle  area  could  not  be  measured  directly.  Different 
DC  gain  levels  are  characteristic  of  all  of  the  exhaust  actuator  frequency  data. 
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Compressor  discharge  pressure,  Pg,  response  of  the  engine  and  the  two 
models  at  95  percent  maximum  spool  speed  are  shown  in  Figure  21.  These 
responses  also  show  the  DC  gain  mismatch  discussed  above.  Other  than  that, 
the  frequency  responses  of  the  two  analytic  models  closely  approximate  the 
frequency  response  of  the  engine.  This  result  is  also  valid  for  the  frequency 
responses  at  the  other  operating  points,  70  percent  and  85  percent  maximum 
spool  speed. 

The  turbine  discharge  pressure,  Pg,  frequency  responses  are  presented  in 
Figure  22.  They  show  that  the  P5  simulations  included  in  the  models  are 
not  as  accurate  as  the  Pg  simulations.  The  gain  and  phase  of  the  component 
model  Pg  response  have  the  same  basic  shape  as  the  corresponding  gain  and 
phase  of  the  engine  response,  but  the  component  model  gain  is  about  10  deci- 
bels lower  and  the  component  model  phase  response  shows  about  30  degrees 
more  phase  shift.  Neither  the  gain  nor  phase  of  the  analog  model  Pg  re- 
sponse matches  the  engine  response. 

The  turbine  discharge  temperature,  Tg,  response  of  the  engine  is  compared 
with  the  corresponding  responses  of  the  two  analytic  models  in  Figure  23. 
These  responses  show  that  the  component  model  approximates  the  engine 
response  well,  especially  in  the  frequency  range  below  2.  0 Hz.  The  analog 
model  also  reasonably  approximates  the  engine  response;  however,  the  analog 
model  gain  response  is  quite  low  and  the  phase  response  shows  more  phase 
shift  than  the  engine. 


ENGINE  RESPONSES  TO  COMPRESSOR  BLEED 

Sample  Bode  frequency  response  plots  of  spool  speed,  compressor  discharge 
pressure,  and  turbine  discharge  temperature  for  oscillations  in  compressor 
bleed  are  presented  in  Figures  24,  25,  and  26.  The  plots  include  frequency 
data  representing  the  engine  and  the  linearized  NASA  component  model;  the 
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analog  model  is  not  represented,  since  the  simulation  does  not  contain  com- 
pressor bleed  effects.  Comparison  of  the  engine  frequency  responses  with 
the  component  model  responses  supports  the  following  conclusions: 

1)  The  shape  of  the  gain  responses  of  the  component  model 
agrees  with  the  shape  of  the  corresponding  engine  responses. 
However,  there  are  some  significant  differences  in  DC  gain 
levels. 

2)  Phase  responses  of  the  component  model  do  not  accurately 
represent  the  engine  phase  responses.  The  engine  frequency 
responses  exhibit  more  phase  shift  than  the  component  model 
responses  throughout  most  of  the  frequency  range. 

Individual  Bode  plots  are  discussed  in  the  following  paragraphs. 

The  spool  speed,  N,  frequency  response  of  the  engine  at  95  percent  maximum 
spool  speed  is  compared  with  the  corresponding  frequency  response  of  the 
component  model  in  Figure  24.  Agreement  between  the  engine  and  component 
model  gain  responses  is  very  good,  within  3 decibels,  but  the  engine  phase 
response  shows  much  more  phase  shift  at  high  frequencies  than  the  component 
model  phase  response,  indicating  that  the  engine  contains  high-frequency 
dynamics  which  are  not  included  in  the  linear  model. 

Engine  and  component  model  frequency  responses  of  Pg  and  Tg  at  95  percent 
maximum  spool  spee^  are  presented  in  Figures  25  and  26.  Both  plots  show 
the  engine  to  have  30  degrees  more  phase  shift  at  low  frequency  than  the 
component  model.  Gain  responses  of  the  engine  and  component  model  are 
similar  in  shape,  but  the  DC  gain  levels  are  off  by  more  tnan  P decibels. 
Frequency  responses  of  these  variables  obtained  at  the  other  ,.vo  operating 
points,  70  percent  and  85  percent  maximum  spool  speed,  further  substantiate 
these  observations. 
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ENGINE  RESPONSES  TO  INLET  GUIDE  VANES 


Frequency  response  plots  of  spool  speed,  compressor  discharge  pressure, 
and  turbine  discharge  temperature  for  oscillation  of  the  inlet  guide  vanes  are 
presented  in  Figures  27,  2 8,  and  2 9.  Engine  frequency  response  data  and 
linearized  NASA  component  model  frequency  response  data  are  compared  in 
the  plots  at  one  operating  point,  95  percent  maximum  spool  speed.  Analog 
model  data  are  not  included,  since  the  analog  simulation  does  not  contain 
inlet  guide  vane  effects.  Comparison  of  the  engine  and  component  model  fre- 
quency responses  supports  the  following  conclusions: 

1)  The  spool  speed  and  F^  responses  of  the  component  model 
are  approximate  representations  of  the  engine  frequency 
responses.  The  shape  of  the  gain  curves  agree,  and  the 
phase  responses  agree  within  30  degrees  in  the  frequency 
range  below  2. 0 Hz. 

2)  The  T&  response  of  the  component  model  is  a poor  repre- 
sentation of  the  corresponding  engine  response.  The  gain 
responses  differ  significantly  in  shape  and  the  phase  re- 
sponses differ  by  as  much  as  140  degrees  in  the  frequency 
range  tested. 

These  conclusions  are  explained  in  the  following  paragraphs. 

Engine  and  component  model  spool  speed  frequency  responses  are  presented 
in  Figure  27.  These  data  show  agreement  typically  within  3 decibels  in  the 
gain  responses,  but  the  phase  shift  of  the  engine  above  0.  5 Hz  is  significantly 
greater  than  the  phase  shift  of  the  component  model.  The  greater  phase  shift 
of  the  engine  data  indicates  that  the  engine  contains  high-frequency  dynamics 
which  are  not  represented  in  the  component  model.  As  previously  noted,  a 
similar  conclusion  was  made  concerning  the  spool  speed  responses  to  oscilla- 
tions in  exhaust  area  and  compressor  bieeds. 
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The  Pg  responses  presented  in  Figure  2 8 show  that,  except  for  a different 
DC  gain  level  and  a constant  phase  error  of  about  20  degrees,  the  component 
model  is  a good  enough  approximation  for  design  purposes  to  the  engine  fre- 
quency response,  considering  the  nonlinearity  of  the  actuator. 

The  frequency  response  data  of  Figure  2 9 show  that  inlet  guide  vane  effect  on 
Tg  is  not  correctly  simulated  in  the  component  model.  Neither  the  gain  nor 
phase  responses  of  the  component  model  agree  with  the  corresponding  re- 
sponses of  the  engine. 
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E1(0=E1  SIN  (w1t+  0X> 
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AMPLITUDE  RATIO  = |^| 

PHASE  SHIFT  = ^ 


Figure  4.  Frequency  Response  Analysis 
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Figure  5.  Fuel  Flow/ Fuel  Command,  W./u, 


0 VC-  II-  881-  **l-  8*1-  »U-  til-  *»«-  lit-  Otc- 


•K 

00*001 

I 

1 

i 

1 2 1 

19*26 

I 

4 

t 

1 

y 

1 S9*t| 

ZZ*S« 

I 

1 

I 

l *****  1 

C*»*6Z 

I 

4 

1 

I 

9 

I U'll 

9S*CZ 

1 

t 

I 

1 **•*.! 

CI*0V 

1 

4 

t 

i 

9 

I toO*"l 

0 1 * C9 

1 

t 

I 

• 

1 0C  ** 

Cto*0S 

I 

4 

I 

9 

I 1 9*  b 

i l**»S 

1 

8 I ‘ 

i ■ 

• 

V 1 8t»*i 

2 1 * OS 

I 

1 

i i 

I 6C*i 

2to*9to 

1 

4 

8 t 

i 

9 

to 

I ton** 

bb’ll 

1 

t 

I tot** 

l%*6C 

1 

4 

8 t 

• • 1 

y 

to 

I ZW*s 

Z«*9C 

1 

1 

I 

I tto’s 

Sl’toC 

1 

4 

a i* 

1 

9 to 

1 CQ*^ 

29*  IC 

I 

t 

I 

I 99*  * 

62*62 

I 

4 

t 

I 

9 

1 2t**> 

21  *Z2 

1 

• 

8 

t 

to 

I 00** 

2!  *S2 

1 

t 

I 

I OZ’t 

92*C2 

l 

• 

! 

9 

I Cto*» 

toS*  12 

I 

1 

t 

1 8|*« 

S6*6 1 

I 

•<f  8 

1 

y to 

I *fc'  € 

09*01 

I 

X 

I 

1 * * c 

n*zi 

1 

d 

X 

1 

9 

1 ef< 

S8*SI 

I 

X 

1 

I *C*c 

©9  * to  l 

I 

X 

I 

I 9 1 * < 

6S*CI 

I 

• d e 

t 

9 to 

I oov 

6S*  2 1 

I 

I 

I 

I 9H*  1 

99*11 

l 

4 

I 

1 

9 

1 _ II  M 

08*01 

I 

t 

1 

i-2  bs- 1 

00*01 

I 

• 

1 

1 

I z*>*  1 

92*6 

I 

1 

I 

1 /EM 

«S*b 

I 

4 

1 

1 9 

i ^2*  1 

tot’Z 

I 

X 

I 

1 Zl*  i 

9C*Z 

I 

X 

1 

I HC*  1 

10*9 

I 

. d 8 

• 

Olto 

I 0«*  « 

IC  *9 

1 

0 

1 

to 

I if 

90*S 

I 

I 

I 

I 98’ 

|9*S 

1 

• e 

• 

VI 

I 00* 

10*S 

1 

d 

I 

9 1 

I *Z* 

*1 9**» 

I 

i 

t 

V 1 

1 «<<• 

oe*9 

I 

: 

l 

u 1 l«‘ 

0to*C 

1 

• 8 

• 

to  I 

1 bS* 

69*c 

1 

: 

I 

0 I is* 

I9*C 

I 

. 8 

• 

to  1 

5 1 OS’ 

9 1 * C 

I 

d 

« 

0 I 

1 /-• 

Ct  * 2 

I 

i 

1 

I l** 

1 Z *2 

1 

• 8 

• 

to  I 

I 0*»* 

18*2 

1 

i 

1 

i a* 

CC*2 

I 

i 

I 

I tot* 

SI  *2 

I 

X 

l 

1 2C* 

00*2 

1 

• to 

• 

V 1 

1 62* 

SH’I 

I 

X 

1 

1 Z2* 

I z • « 

I 

d 

to 

1 

I sz* 

bS*  1 

I 

i 

X 

I C7* 

Z*>*  1 

I 

l 

I 

1 if 

9t*  I 

I 

• 8 

•v 

I 

I 02* 

92*1 

1 

1 

1 

1 b\ • 

Z I * l 

1 

t 

I 

» 0 zr 

HO*  | 

1 

1 

l 

I-^bl- 

00*1 

I 

• 

• 

I 

l SI  * 

C6* 

1 

1 

1 

1 *1  • 

90  * 

I d 

9 

I 

CO 

1 Cl* 

t>Z* 

1 

1 

I 

< 

* a co 

1 zr 

*z* 

i 

1 

1 

£ 

1 ii* 

09* 

I 

• 8 

to 

I 

I 01* 

C9  * 

i 

8 

to 

I 

Z 

I <»0* 

06* 

I 

i 

l 

< 

• o< 

I M>* 

9S* 

i 

8 

to 

1 

<3 

1 «»* 

OS* 

I 

i 

I 

1 Z«* 

9«»* 

i 

e 

to 

1 

oi 

1 Z 0* 

C*»* 

1 d 

9 

1 

>-  z 

J 000 

1 *»0* 

0*»* 

i 

8 

to 

I 

U) 

* 

H 

I *»1* 

zc* 

i 

t 

1 

I so* 

toC* 

I 

8 

to 

1 

-J 

1 so* 

2C  * 

I 

* 

1 

jjL 

&! 

SujS 

I so* 

62* 

I 

I 

Id 

I 

b 

i 

to 

9 

1 

I 

I 

I 

1 

1 too* 

I *0* 

I *0* 

1 CO* 

z2* 

S2* 

Ci* 

22* 

I 

t 

I 

UQo 

I CO* 

02* 

I 

l 

I 

zgS 

1 CO* 

01* 

I 

1 

I 

Si  5 

I iC* 

Zl* 

Id 

9 

1 

ujO< 

I CO* 

91* 

i 

t 

I 

I If 

SI* 

I 

1 

1 

l 20* 

to|* 

I 

l 

I 

1 !»• 

Cl* 

Id 

9 

I 

1 a- 

21* 

I 

I 

1 

1 ^ 20* 

II* 

o*si  s**’ 

»•* 

s*z- 

0*si-  S*22- 

>:ti- 

S*ZC-  >K i to9  Zh 

J3s/av« 

ADNinMBi  *SA  3S*hb  on*  M*9  JC  iOld 


34 


Figure  6.  Exhaust  Area  /Exhaust  Command, 
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Figure  7.  Compressor  Bleed /Bleed  Command 
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Figure  8.  Inlet  Guide  Vane /Inlet  Guide  Vane  Command,  IGV/u. 
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Figure  9.  Spool  Speed/Fuel  Flow,  N/W^--95  and  100  Percent  Maximum  Speed 


• H- 

It* 

HI-  Wt- 

•at- 

tu-  iu-  aa<-  »«•  •«- 

• 

* 

• 

• 

i 

• 

• 

] 

; 

1 

---o i*p  (i  **9i 

aa*aaai 

1 

i 

i 

I oH‘l*l 

ll'W 

1 

i 

i 

l * is*vci 

94 *4  Si 

I 

4 

i 

i 

» 1 I*>H 

(CH4 

1 

i 

t 

i aa*iti 

9V*SC4 

1 

4 

f 

t 

o i tvvai 

H*IH 

1 

i 

i 

i i**«ai 

H*l« 

1 

4 

t 

i 

l aa'u 

vt*vas 

1 

t 

i 

I cr*t 

4I*I9S 

I 

4 

i 

> i 

•1*194 

1 

i 

t 

1 tv*w 

9|*V«9 

1 

i 

t V 

I 

49**49 

1 

i 

t 

1 *fl* 

!!*••£ 

1 

i 

1 V 

1 aa*»t 

•i’tvc 

i 

t 

I 91*99 

SV*|9C 

I 

i 

1 0 

1 ic*«>9 

£1*9 II 

1 

i 

1 

1 !•*** 

V9*li4 

I 

i 

1 V 

i art* 

C**l4? 

I 

i 

t 

l ta'oc 

•1*144 

1 

i 

1 « 

1 19*4C 

C9*1C4 

1 

i 

t 

i wn 

9VSI4 

I 

i 

1 

1 94  Me 

(•*••! 

i 

§4 

i 

t V* 

1 |9*6l 

94*901 

t 

t 

1 

(l*lil 

i 

t 0 

i 

•9*0SI 

i 

t 

1 vC*»l 

•4*991 

t 

i y 

1 (9*14 

9»*SC1 

1 

• 

t 

1 V 

i »a*'/ 

•V*S4I 

1 

i 

i 0 

1 V9*"l 

•S*  9 1 1 

i 

i 

1 _6l*il 

••*491 

i 

l 0 

1-5  41*  •'■-t 

•9*991 

I 

i 

i 

1 -«*i**l 

19*19 

1 

i 

1 0 

1 99*l| 

44*99 

i 

1 

1 *<**<! 

(9*64 

i 

1 0 

1 I4*H 

V9*C4 

I 

i 

1 

1 99*.'| 

ci*  yv 

1 

44 

t 

1 0 V 

l *a*<  i 

• rev 

t V 

1 91*6 

(9*99 

1 

t 

1 

l IV* 

Cl*  99 

l 

1 0* 

1 96*4 

11*09 

I 

» 

1 

1 6 1 *1 

49*V9 

1 

44 

• 

i 

1 0* 

1 9t»*>- 

16*49 

1 

t 

1 

1 91  *- 

U'H 

1 

44 

i 

1 • V 

1 49*9 

4V*  VI 

1 

i 

1 

1 19*S 

Sl*9l 

1 

44 

• 

t 

1 • » 

1 (6*s 

4V*ll 

1 

i 

1 

1 V***- 

64*64 

1 

i 

i u 

1 Cl  **» 

11*44 

1 

4 

i 

i • * 

1 *#** 

11*94 

1 

t 

1 0 

1 94*1 

vl*C4 

1 

t 

1 

1 (9*1 

99*»1 

1 

i 

1 V 

1 91*1 

99*61 

1 

1 

t 

1 • * 

1 9fc*c 

99*«l 

1 

i 

1 

1 44** 

11*41 

I 

i 

1 9 

1 49*< 

9V*9 1 

1 

i 

1 

l 91*. 

•V  * 9 | 

1 

t 

1 *J 

1 9I*« 

•9*C  1 

1 

• 

t 

1 • V 

1 99*i 

69*11 

1 

4 

i 

1 0 

O l V9* 1 

WII 

k*J  1 44*1 

09*01 

1 

t 

1 9 

sc  i-°  »s’i 

00*01 

1 

i 

O I 49-  1 

V*  * 6 

1 

i 

1 

Si  4i*  * 

• 9*M 

1 

4 

i 

IV 

1 94*1 

96*4 

1 

i 

I 

1 41*1 

Vt*4 

1 

i 

9 

1 *9*1 

19*9 

1 

4 • 

i 

• i* 

1 C9*l 

II  *v 

1 

a 

i 

• 

1 16* 

99*9 

1 

i 

V l 

1 vy* 

|9*9 

1 

a • 

i 

• V 

1 0»* 

10*9 

i 

1 

1 94  * 

9V*9 

1 

a 

i 

9 VI 

1 99* 

01*9 

1 

i 

1 

1 IV* 

96*C 

1 

a • 

i 

• V 1 

1 69* 

6 V * I 

1 

i 

1 

1 9%* 

|9  * C 

1 

u • 

9V  1 

1 04* 

Vl*l 

I 

i 

i 

1 49* 

C6*4 

1 

t 

l 

1 C 9 * 

14*1 

1 

a • 

V 1 

1 99* 

19*4 

1 

i •. 

1 4 C * 

CI  *1 

I 

i 

1 

1 91* 

91*1 

1 

i 

1 

1 41* 

00*4 

I • 

• 

V 1 

1 64* 

99*  1 

i o 

1 

1 4 4’ 

14*1 

I 

1 

1 94* 

99*1 

g 

1 

1 14* 

49*  I 

1 

i 

t 

1 c4* 

V(*l 

1 t • 

a 

1 

i ac* 

94*  1 

I 

i 

1 

1 6l  * 

41*1 

l • 

1 9 

1 

1 © *«* 

-0*1 

1 

1 

1 £1 

1-  VI* 

00*  1 

1 <1  « ata  1 41  * 

16* 

1 

t 

> £1 

1 9 | * 

Vta* 

I 

I 

1 

1 (!* 

64  * 

1 

i §1  •«* 

1 41* 

94* 

• V* 

1 %4 

10* 

i at* 

tv* 

1 ■ > 

* 1 * 

t 

| #.** 

• 9* 

1 

I • 

I 

1 V 

; .i  ssfi 

I ui  * cocoeo  J 

99* 
69* 
V9  * 

• 1 * 

, K 

1 If 

19* 

1 

1 

1 -J 

1 99* 

09* 

i a 4 

1* 

1 o 

1 Vf  * 

4t  * 

i 

1 

1 C 

-1  1 49* 

91* 

41  * 

* 

R 1 so* 

i 

1 

Si 

1 1 99* 

44* 

• IV 

O 1 9t* 

o I 99* 

94* 

1 

i 

14* 

1 4 

10 

‘ 3? 

< l Ci* 

14* 

1 

1 

If 

! -a  ! 

91* 

1 

1 

i 

1 (9* 

4 1* 

10 

i 

1 ta* 

VI* 

1 

i 

1 49- 

91* 

1 

i 

1 49* 

9 1* 

1 

i 

1 49* 

(I* 

10 

i 

1 49* 

41* 

1 

1 

i 

1 rV49* 

II* 

• ••4  • 

•l  a 

•*« 

••a  a**i- 

9*49-  >MV*T  /►* 

339/uve 

ONinnJMJ  *9*  J9VM*  GN*  MVJ  <0  101^ 


38 


Figure  10.  Spool  Speed/Fuel  Flow,  N/W, — 85  Percent  Maximum  Speed 
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Figure  11.  Spool  Speed/Fuel  Flow,  N/W^ — 70  Percent  Maximum  Speed 
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Figure  12.  Compressor  Discharge  Pressure/Fuel  Flow,  P3/Wf — 
95  and  100  Percent  Maximum  Speed 
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Figure  15.  Turbine  Discharge  Pressure /Fuel  Flow,  Ps/W^ 
95  and  100  Percent  Maximum  Speed 
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Figure  16.  Turbine  Discharge  Pressure/Fuel  Flow,  Ps/Wp 
85  Percent  Maximum  Speed 
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Figure  18.  Mach  Number  Sensor/ Fuel  Flow, 
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Figure  19.  Turbine  Discharge  Temperature /Fuel  Flow, 
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Figure  20.  Spool  Speed /Exhaust  Area, 
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Figure  21.  Compressor  Discharge  Pressure/ Exhaust  Area, 


CO 


in 

Pu 

rt 

v 

u 

< 

-*-> 

CO 

3 

rt 


W 

<u 

3 

CO 

CO 

V 

PU 

0) 


335>/U*H 


*3*3003* J ”SA  3S¥M<I  ONV  M*‘J  30  i03d 


50 


Figure  22.  Turbine  Discharge 
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Figure  24.  Spool  Speed/Compressor  Bleed,  N/BLD 
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Figure  25.  Compressor  Discharge  Pres  sure /Compressor  Bleed, 
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Figure  26.  Turbine  Discharge  Temperature/ Compressor  Bleed,  T5/BLD 
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Figure  27.  Spool  Speed  /Inlet  Guide  Vane,  N/IGV 
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Figure  28.  Compressor  Discharge  Pres  sure /Inlet  Guide  Vane,  P,/IGV 
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Figure  29.  Turbine  Discharge  Temperature /Inlet  Guide  Vane,  T-/IGV 
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Figure  30.  Linearized  Component  Model 


Table  1.  Steady-State  Data--Engine  Test 


Response 

Units 

Operating  Point,  N/Nmax 

70  Percent 

85  Percent 

95  Percent 

Wf 

ib/hr 

630.  0 

910.0 

1504.0 

A8 

in2 

162.  0 

158.0 

115.0 

IGV 

Nondimens  ional 

0.  943 

0.600 

0.049 

BLD 

Nondimensionai 

0.983 

0.638 

0. 115 

N 

rpm 

11,319.  0 

13,  976.0 

15,279.0 

P3 

psi 

34.4 

53.2 

74.8 

P5 

psi 

15.  8 

18.0 

24.1 

4P/P3 

psi/psi 

0. 125 

0.  152 

0.  161 

T4 

°R 

1331.0 

1377. 0 

1632.0 

T5 

°R 

990.0 

957.  0 

1141.0 

APfn 

psi 

135.6 

146.  8 

198.2 

Table  2.  Steady-State  Data --Component  Model 


Response 

Units 

Operating  Point,  N/Nmax 

70  Percent 

85  Percent 

100  Percent 

Wf 

lb/hr 

738.  0 

1052. 0 

2927.  0 

A8 

in2 

132.  0 

162.  0 

91.  5 

IGV 

Nondimensionai 

1.  0 

0.  619 

0 

BLD 

Nondimensionai 

1.  0 

0.  724 

0 

N 

rpm 

11,  550.  0 

14,025.0 

16,  500. 0 

P3 

psi 

39.  9 

60.4 

101.  5 

P5 

psi 

17.  2 

20.  3 

36.  6 

T4 

°R 

1426.  0 

1498.  0 

2216. 0 

T5 

°R 

1194. 0 

1179.  0 

1796.0 
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SECTION  m 

MULTIVARIABLE  DYNAMIC  ENGINE  MODEL 


INTRODUCTION 

Complete  transfer  function  models  corresponding  to  the  Bode  frequency  re- 
sponse data  obtained  from  the  APL  J85  engine  are  presented  and  discussed 
in  this  section.  Consistency  and  physical  significance  of  the  models  are 
analyzed  in  the  following  paragraphs.  The  results  demonstrate  that  multi- 
variable  dynamic  engine  models  can  be  successfully  identified  from  experi- 
mental frequency  response  measurements. 

Transfer  function  models  were  identified  from  the  Bode  frequency  response 
plots  with  a computer  algorithm  called  TFNS.  The  algorithm  is  programmed 
to  calculate  the  transfer  function  which  best  approximates  the  experimental 
frequency  response  measurements.  Polynomials  in  the  transfer  function  are 
computed  iteratively  so  that  the  square  of  the  difference  between  the  experi- 
mental frequency  response  and  frequency  response  of  the  transfer  function 
approximation  is  minimized.  That  is,  if  the  experimental  frequency  response 
is  denoted  by  G(ju)  and  the  frequency  response  of  the  transfer  function  approxi- 
mation is  denoted  by  G (ju),  the  computer  algorithm  calculates  the  poly- 

a 

nomials  in  Ga(ju)  so  that  the  following  error  is  minimized: 

00 

E^J  |G(jw)  - Ga(ju)|2  dw 


A detailed  discussion  of  the  TFNS  program  is  included  in  Appendix  B.  The 
trane  '?r  function  models  identified  with  the  TFNS  program  are  listed  in 
Tables  4 through  23.  Models  are  presented  for  three  engine  speed  operating 
points,  N/Nmax  = 70  percent,  85  percent,  and  95  percent.  The  transfer 
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functions  included  in  these  data  are  listed  below  in  the  same  order  ir.  which 
they  are  discussed  in  the  subsections  that  follow.  Actuator  inputs  arc  sensor 
measurements  represented  in  these  transfer  functions  are  described  1 
Section  II. 


Ref. 

Table 


Transfer  Function 


Actuator  4 

transfer  functions 

5 


6 


7 


Fuel  flow  8 

transfer  functions 

9 


10 


11 

12 


Exhaust  area  13 

transfer  functions 

14 


15 


16 

17 


Fuel  flow/fuel  command,  W^/u^ 

Exhaust  area/exhaust  area  command, 

VUAg 

Compressor  bleed/compressor  bleed 
command,  BLD/ug^j) 

Inlet  guide  vane/inlet  guide  vane  command, 

igv/uicv 

Spool  speed /fuel  flow,  N/Wf 
Compressor  discharge  pressure /fuel  flow, 

Vwf 

Turbine  discharge  pressure/fuel  flow, 

p5/w. 

A P 

Mach  number  sensor /fuel  flow,  /W. 

3 1 

Turbine  discharge  temperature /fuel  flow, 

Vwt 

Spool  speed/exhaust  area,  N/Ag 
Compressor  discharge  pres  sure /exhaust  area, 

*VA8 

Turbine  discharge  pres  sure /exhaust  area, 

*VA8 

, Ap 

Mach  number  sensor /exhaust  area,  -5—  /A« 

P3  8 

Turbine  discharge  temperature /exhaust 
area,  T5/AQ 


62 


Ref. 

Table 

Transfer  Function 

Compressor  bleed 

18 

Spool  speed/compressor  bleed,  N/BLD 

transfer  functions 

19 

Compressor  discharge/compressor  bleed, 

p3/bld 

20 

Turbine  discharge  temperature /compressor 
bleed,  Tg/BLD 

Inlet  guide  vane 

21 

Spool  speed/inlet  guide  vane,  N/IGV 

transfer  functions 

22 

Compressor  discharge  pres  sure /inlet 
guide  vane,  Pg/IGV 

23 

Turbine  discharge  temperature /inlet  guide 
vane,  Tg/IGV 

The  quality  of  these  transfer  function  models  is  summarized  in  the  following 
statements; 

• Engine  dynamics  in  the  0.  05  to  100- Hz  frequency  range  are 
identified  in  the  models. 

• Consistency  is  demonstrated  between  corresponding  transfer 
functions  at  different  operating  points. 

• Consistency  is  demonstrated  between  different  transfer  functions 
at  the  same  operating  point. 

These  statements  are  briefly  discussed  in  the  following  paragraphs. 

The  transfer  functions  presented  in  Tables  4 through  23  contain  a compre- 
hensive, accurate  description  of  the  J85  engine.  All  significant  actuator 
dynamics  and  gas  dynamics  in  the  frequency  range  0.  05  to  100  Hz  are  accu- 
rately represented  in  the  data.  In  the  past  it  has  beenpossible  to  experiment- 
ally measure  engine  dynamics  only  over  a very  limited  frequency  range  due 
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to  the  difficulty  in  computing  accurate  amplitude  and  phase  shift  information 
from  noisy  sensor  signals.  This  difficulty  was  circumvented  with  the  BAFCO 
two- channel  servoanalyzer  which  is  unique  in  its  ability  to  accurately  identify 
frequency  response  data  from  noisy  signals. 

Examination  of  corresponding  transfer  functions  at  different  operating  points 
shows  the  results  to  be  consistent  with  respect  to  the  relative  positions  of  the 
poles  and  zeros.  As  an  example,  consider  the  transfer  function  models  of 
Tg/Ag  at  the  three  operating  points,  N/Nmax  = 70  percent,  85  percent,  and 
95  percent.  These  three  transfer  functions  are  all  first  over  second  order. 
Poles  and  zeros  in  the  transfer  functions  are  listed  below: 

Tc/A0  Transfer  Function  Models 
0 o 


Association 

N'Nma*  1 70* 

85% 

95% 

Zeros  (rad/sec) 

-8.75 

CO 

H 

1 

-42.4 

f -0.361 

-1.01 

eg 

r-H 

1 

Poles  (rad/sec)  < 

1 

1 -3.67 

-8. 14 

-13.0 

The  data  clearly  show  the  location  of  the  poles  and  zeros  to  shift  consistently 
between  operating  points.  That  is,  the  poles  and  zeros  for  the  N/N^^  = 

70  percent  operating  point  are  all  located  at  lower  frequencies  than  the  cor- 
responding poles  and  zeros  for  the  85  percent  operating  point.  Similarly, 
the  poles  and  zeros  for  the  85  percent  operating  point  are  all  located  at 
lower  frequencies  than  the  corresponding  poles  and  zeros  for  the  95  percent 
operating  point.  This  shows  consistency  with  the  ohysics  of  turbine  engines. 

Consistency  between  different  transfer  functions  at  the  same  operating  point 
is  also  exhibited  by  the  data.  Dynamic  states  are  properly  identified  in  sev- 
eral transfer  functions.  For  example,  consider  the  three  transfer  functions 

N/W-,  N/Aq,  and  P0/A0  for  the  N/N,„„  = 95  percent  operating  point.  The 
i o o o max 

poles  in  these  transfer  functions  are  listed  below: 
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Comparison  of  Transfer  Function  Models  for 
N/Nmax  = 95  Percent 

Association  N/W^  N/Ag  Pg/Ag 

Poles  (rad/sec)  -2.  86  -3.38  -3.41 

-76.7  -53.9  -45.8 

-260.0  

These  data  show  a pole  located  at  approximately  -3.  0 radians  per  second 
common  to  all  three  transfer  functions.  Another  pole  located  at  about  -50 
radians  per  second  is  common  to  both  the  N/Ag  and  Pg/Ag  transfer  functions. 
Possibly  the  pole  located  at  -76,7  radians  per  second  in  the  N/W^  transfer 
function  also  corresponds  to  the  -50  radians  per  second  pole  identified  in  the 
other  two  transfer  functions.  Some  thermodynamic  states  should  appear  in 
several  output  responses  in  this  way. 

In  many  cases  the  state  variables  identified  in  the  transfer  function  models 
can  easily  be  interpreted  in  terms  of  physical  engine  parameters.  The  most 
obvious  of  these  relationships  are  summarized  below: 

• Spool  inertia  shows  up  in  all  the  transfer  functions. 

• A time  delay  has  been  identified  in  the  P0/W»,  Pc/W,,  and 

Ap.  3 f 5 f 

transfer  f motions. 

0 1 

• Tg  thermocouple  time  constants  are  contained  in  the  Tg 
response, 

• Tailpipe  gas  dynamics  show  up  in  the  fuel  flow  respones, 
as  well  as  in  some  of  the  other  responses. 

All  of  the  transfer  functions  contain  a low-frequency  pole  between  about  0.  7 
and  3.  5 radians  per  second  which  is  associated  with  spool  inertia.  Approxi- 
mate time  constants  are  identified  below  for  each  of  the  three  operating 
points  modeled: 
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Spool  Inertia  Time  Constants  in  Seconds 
Association  = 70%  85%  95% 

IIlclX  

Time  constant  1.2  0.55  0.36 

Range  of  data  1.4  to  0.5  0.85  to  0.27  0.65  to  0.20 

These  data  were  computed  by  averaging  the  results  for  each  operating  point. 

The  Pg/Wf,  Pg/Wj.  and-~/Wf  transfer  functions  contain  Pade  approxima- 
tions corresponding  to  a time  delay  of  between  11  and  15  milliseconds. 
Although  this  time  delay  is  believed  to  be  associated  with  the  fuel  combustion 
process,  the  specific  physical  interpretation  of  the  delay  is  not  clear  at  this 
time.  This  delay  is  considerably  larger  than  the  combustion  time  delay  pre- 
dicted by  existing  engine  models.  The  time  delay  data  are  summarized 
below: 


Combustion  Time  Delay  in  Seconds 
Association  N/N  = 70%  85%  95% 

Delay  time  0.015  0.014  0.011 

Identification  of  these  time  delays  from  the  transfer  function  data  is  dis- 
cussed later  in  this  section  under  the  heading  "Engine  Transfer  Functions 
for  Fuel  Flow. " 

The  Tg  transfer  functions  contain  a set  of  poles  in  the  0.  35  to  1. 90- radian 
per  second  frequency  range  which  represents  Tg  thermocouple  dynamics. 
Approximate  time  constants  are  identified  below  for  each  of  the  three 
operating  points  modeled: 
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Tg  Thermocouple  Time  Constants  in  Seconds 

Association  N/Nmax  = 85%  95% 

Time  constant  2.71  1.21  0.60 

Range  of  data  2.  53  to  2.  85  0.  87  to  2. 52  0.  53  to  0.  70 

These  data  were  computed  by  averaging  the  results  for  each  operating  point. 

The  fuel  flow  transfer  functions  also  contain  a consistent  set  of  poles  which 
represent  the  fundamental  resonance  in  the  engine  tailpipe.  The  resonance 
is  characterized  by  a natural  frequency  of  between  35  and  50  radians  per 
second  with  a damping  coefficient  around  0. 2. 


ACTUATOR  TRANSFER  FUNCTIONS 

Transfer  function  models  corresponding  to  the  four  engine  actuators  are  pre- 
sented in  Tables  4 through  7.  Poles  and  zeros  for  the  fuel  valve  are  identi- 
fied out  to  about  50  Hz.  The  geometry  actuators,  exhaust  area,  compressor 
bleeds,  and  inlet  guide  vanes  are  lower  bandwidth  than  the  fuel  valve  and, 
consequently,  the  poles  and  zeros  for  these  actuators  are  identified  only  out 
to  about  10  Hz. 

Comparison  of  corresponding  actuator  transfer  functions  at  different  operat- 
ing points  shows  that  the  actuators  respond  more  quickly  at  spool  speeds  near 
maximum  than  they  do  at  lower  speeds  (i.  e. , the  natural  frequencies  of  the 
actuators  increase  with  increasing  spool  speed).  This  behavior  is  related  to 
the  hydromechanical  design  of  the  actuators.  The  actuators  respond  faster 
at  higher  spool  speeds  because  the  hydraulic  pressure  power  source  in- 
creases with  spool  speed. 
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Transfer  function  models  of  the  fuel  valve  at  the  three  operating  points  tested 

are  listed  in  Table  4.  Both  the  95  percent  and  85  percent  N/N  models 

max 

contain  six  orders  of  numerator  and  denominator  dynamics;  the  70  percent 

N/N  model  contains  five  orders  of  numerator  and  demoninator  dynamics, 
max 

Poles  associated  with  the  models  are  identified  below. 


W^/Uj  Dynamics 


N/Nma* 

= 70% 

85% 

95 

% 

Association 

Frequency 

Damping 

Frequency 

Damping 

Frequency 

Damping 

Tailpipe 

dynamics 

34.2 

0.502 

44.9 

0.206 

49.  4 

0.242 

Actuator 

1 148.0 

0.208 

137.  0 

0.  145 

134.  0 

0.  172 

dynamics 

| -156.0 

— 

239.  0 

0.  179 

248.  0 

0.  172 

These  models  contain  representations  of  tailpipe  resonance  and  actuator  dy- 
namics. Ideally,  the  models  would  only  contain  actuator  dynamics,  but  be- 
cause fuel  flow  was  obtained  by  measuring  the  pressure  differential  across 
the  fuel  nozzle,  (APfuel  nozzle  4 Pfuel  nozzle  - Pg),  these  models  contain 
dynamics  identified  with  the  Pg  response  in  addition  to  actuator  dynamics. 

A detailed  analysis  of  the  effect  of  P3  dynamics  on  the  measurement  of  fuel 
flow  is  included  in  Appendix  A. 


Resonance  in  the  tailpipe  section  is  represented  by  a complex  pole  with 
natural  frequency  around  40  radians  per  second  and  a damping  ratio  of  0. 25. 
The  natural  frequency  of  the  pole  is  seen  to  increase  with  increasing  spool 
speed  because  the  velocity  of  the  airflow  in  the  tailpipe  increases  with  spool 
speed.  Similar  representations  of  tailpipe  dynamics  are  included  in  the 


P3/W(. 


P5/Wf. 


A P 

and^-/W.  models  discussed  in  the  following  suosection. 
*3  1 


Except  for  the  70  percent  model,  fuel  valve  actuator  dynamics  are  repre- 
sented as  a pair  of  complex  poles,  one  with  a natural  frequency  of  140  radians 
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per  second  and  a damping  ratio  of  0.  18,  and  the  other  with  a natural  fre- 
quency of  240  radians  per  second  and  a damping  ratio  of  0, 175.  This  repre- 
sentation of  actuator  dynamics  agrees  with  the  results  obtained  from  bench 
tests  of  the  fuel  valve  which  also  showed  the  valve  to  be  a fourth- order  de- 
vice characterized  by  two  complex  poles,  one  with  a natural  frequency  of  20 
Hz  and  a damping  ratio  of  0.  20,  and  the  other  with  a natural  frequency  of  30 
Hz  and  a damping  ratio  of  0.  80.  " 

Exhaust  actuator  models  are  listed  in  Table  5.  These  models  show  the 
exhaust  actuator  to  be  a second- order  device  with  a natural  frequency  around 
33  radians  per  second  and  a damping  ratio  of  about  0.  40.  The  natural  fre- 
quency increases  with  increasing  spool  speed  from  27.2  radians  per  second 
at  the  70  percent  N/Nmax  operating  point  to  38.  2 radians  per  second  at  the 
95  percent  N/Nrr)ay  operating  point. 

Transfer  function  models  of  the  compressor  bleeds  and  inlet  guide  vane  actu- 
ators presented  in  Tables  6 and  7 show  both  of  these  actuators  to  be  essen- 
tially first-order  devices.  Time  constants  for  the  actuators  are  summarized 
below. 


Bleed  and  IGV  Time  Constants  in  Seconds 


Association 

N /Nmav  = ™% 

max 

85% 

95% 

Compressor  bleed 

0.  032 

0.  029 

0.  028 

time  constant 

Inlet  guide  vane 

-- - 

0.066 

0.  050 

time  constant 

No  data  are  presented  for  the  inlet  guide  vanes  at  N/ Nmax  = 70  percent  be- 
cause the  IGV  actuator  has  little  effect  on  engine  responses  at  this  speed. 
The  time  constants  listed  above  show  the  BLD  and  IGV  actuators  to  be  much 
faster  than  expected;  the  ELD  actuator  was  previously  modeled  with  a time 
constant  of  0.  5 second  and  the  IGV  actuator  was  previously  modeled  with  a 
time  constant  of  0.  2 second. 
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ENGINE  TRANSFER  FUNCTIONS  FOR  FUEL  FLOW 


Transfer  function  models  of  spool  speed  N,  compressor  discharge  pressure, 
P3<  turbine  discharge  pressure,  Pg,  Mach  number  sensor,  AP/Pg,  and 
turbine  discharge  temperature,  Tg,  with  respect  to  fuel  flow,  W^,  are  pre- 
sented in  Tables  8 through  12.  Engine  dynamics  in  the  0.05  to  100- Hz  fre- 
quency range  are  identified  in  the  models,  including  representations  of  spool 
inertia,  gas  dynamics,  sensor  dynamics,  and  combustion  time  delay.  Some 
of  the  important  features  of  these  models  are  summarized  below  and  discussed 
in  the  following  paragraphs: 


• Spool  inertia  is  consistently  identified  in  the  N/W^,  Pg/W^,  and 

a rs 

/Wf  transfer  functions. 

3 1 

• Pade  approximations  to  a combustion  time  delay  are  included 

in  the  Pg/Wf,  Pg/Wf,  and  pP/Wf  models. 

3 

• The  Tg/Wf  transfer  functions  contain  Tg  thermocouple  sensor 
dynamics. 

• Tailpipe  gas  dynamics  are  identified  in  the  N/W,,  PQ  / W„, 

AP  x d x 

Pg/Wf,  and  p /Wf  models. 


Table  8 contains  transfer  function  models  of  N/W^  for  the  three  operating 
points  tested,  70  percent,  85  percent,  and  95  percent  maximum  spool  speed. 
Two  sets  of  poles  are  consistently  identified  in  these  models.  The  first  set 
of  poles  occurs  at  relatively  low  frequency,  around  2.  0 radians  per  second, 
and  is  associated  with  spool  inertia.  The  second  set  of  poles  which  occurs  at 
higher  frequency,  around  60  radians  per  second,  and  is  associated  with  gas 
dynamics  in  the  engine  tailpipe.  Locations  of  the  poles  for  these  three  models 
are  listed  below: 
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N/W^  Dynamics 


Root  Association 

70% 

85% 

95% 

Spool  inertia 

-0.714 

-1.  17 

-2.86 

Tailpipe  dynamics 

-42.8 

-60.  1 

-76.7 

^ mmrnm 

-260.0 

Notice  that  the  location  of  the  poles  varies  consistent'y  with  spool  speed. 
This  result  is  predicted  by  engine  thermodynamics. 

Transfer  function  models  of  Pg/W^  are  presented  in  Table  9.  These  models 
contain  representations  of  spool  inertia,  gas  dynamics,  and  combustion  time 
delay.  All  three  transfer  functions  contain  seven  poles.  Associations  with 
engine  and  sensor  dynamics  are  outlined  below: 


P3/Wf  Dynamics 


7 0% 

85 

01 

10 

95% 

Association 

Frequency 

Damping 

Frequency 

Damping 

Frequency 

Damping 

Spool  inertia 

-0.  800 

--- 

-1.  39 

— 

- 2.83 

— 

Tailpipe 

dynamics 

34.  0 

0.041 

38.  0 

0.  162 

49.6 

0.  021 

Combustion 
time  delay 

|274.  0 

0.  892 

149.  0 

0.812 

290.  0 

0.  853 

l 479.  0 

0.  328 

423.  0 

0.  133 

542.  0 

0.  039 

The  lowest- frequency  pole  at  each  operating  point  is  associated  with  spool 
inertia.  These  roots  agree  with  the  spool  inertia  roots  identified  from  the 
N/W^  transfer  function  models. 

The  second  lowest-frequency  pole  at  each  operating  point  is  a complex  pair 
which  is  associated  with  gas  dynamics  in  the  engine  tailpipe  section.  The 
resonant  frequency  increases  with  increasing  spool  speed  because  the  air- 
flow velocity  increases  with  increasing  spool  speed.  A representation  of 
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these  dynamics  (not  as  accurate)  was  also  identified  from  the  N/Wf  transfer 
functions. 

The  two  high-frequency  complex  pairs  at  each  operating  points  are  part  of  a 
Pade  approximation  of  a time  delay  related  to  the  combustion  process.  The 
Pade  approximation  is  usually  made  up  of  two  or  more  pairs  of  poles  and 
zeros  which  have  the  same  natural  frequency  and  equal  damping  ratios  with 
opposite  algebraic  signs.  For  example,  consider  the  Pade  approximation 
included  in  the  95  percent  maximum  spool  speed  Pg/Wf  transfer  function. 

In  the  frequency  domain  the  Pade  approximation  is: 


CM 

m oo 

2 . 2fc0,  754)S  J 
282  1 

1— 

.1528 

2 , 2(-0.  056)S  , , 
528  J 

[m] 

2 , 2(0.  853)S  . 

+ 290  1 

1 S 
1542J 

2 . 2(0. 039)S  . , ‘ 
+ 542  1 

This  approximation  is  composed  of  two  pole- zero  pairs.  The  zeros  in  each 
pair  are  roughly  identical  to  the  poles  except  for  the  algebraic  sign  on  the 
damping  term.  This  characteristic  identifies  the  function  as  a Pade  approxi- 
mation to  a time  delay. 

Similar  Pade  approximations  are  included  in  the  P^/W^.  models  at  the  other 
two  operating  points,  70  percent  and  85  percent  maximum  spool  speed.  The 
delay  times  in  seconds  corresponding  to  these  approximations  are  listed 
below: 


Combustion  Time  Delays  in  Seconds 
Association  70%  85%  95% 

Delay  time  0.015  0.014  0.011 

Thus,  the  delay  time  decreases  with  increasing  spool  speed. 
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Transfer  function  models  for  Pg/W^  are  presented  in  Table  10.  Engine  dy- 
namics represented  in  these  models  are  identified  below: 


Pg/Wf  Dynamics 


L 

70% 

85% 

95% 

A'-sociation 

■SI 

liiiiutii 

Spool  inertia 

1 

0.291 

-3.  13 
44.7 

0.239 

-8.82 

45.7 

0.268 

Gas  dynamics 

[ 

0.734 

-81.7 
-130.  0 

mm 

123.0 

0.866 

Combustion  time 
delay 

| 

0.  335 

-495.  0 
574.  0 

685.0 

0.741 

The  lowest-frequency  pole  at  each  operating  point  is  believed  to  be  associated 
with  spool  inertia,  although  the  roots  identified  here  are  higher  in  frequency 
than  the  spool  inertia  roots  identified  in  the  N/Wf  and  Pg/Wf  transfer 
functions. 

The  poles  located  in  the  intermediate  frequency  range,  50  to  130  radians  per 
second,  represent  engine  gas  dynamics.  Two  resonance  conditions  are 
included  in  these  models.  The  complex  pair  at  about  45  radians  per  second 
with  a damping  ratio  around  0.  26  is  identified  with  the  fundamental  resonance 
in  the  tailpipe  section.  The  other  complex  pair  at  about  100  radians  per 
second  with  a damping  ratio  around  0,  80  represents  either  a secondary 
resonance  in  the  tailpipe  section  or  a resonance  in  a smaller  volume  in  the 
engine  such  as  the  burner  or  compressor  volume.  It  should  be  noted  that  the 
higher- frequency  resonance  is  not  accurately  identified  in  the  N/Nmay  = 85 
percent  model;  it  is  represented  by  two  real  roots  (natural  frequencies  of 
-81.  7 and  -130  radians  per  second)  instead  of  by  one  complex  pair.  This 
inconsistency  is  apparently  related  to  experimental  error  in  the  raw  fre- 
quency data. 
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The  high-frequency  poles  at  each  operating  point  are  part  of  Pade  approxi- 
mations to  combustions  time  delays.  Although  the  poles  and  zeros  associated 
with  these  Pade  approximations  are  not  identical  to  the  poles  and  zeros 
associated  with  the  Pade  approximations  identified  in  the  Pg/W^  models,  the 
time  delays  modeled  are  the  same.  Thus,  at  70  percent  maximum  spool 
speed  the  delay  is  0.  015  second,  at  85  percent  the  delay  is  0.  014  second, 
and  at  95  percent  the  delay  is  0.  Oil  second. 

Ap , 

Mach  number  sensor,  -^~i Wf,  transfer  function  models  are  listed  in  Table  11. 

3 

The  dynamic  representations  included  in  these  models  are  similar  to  those 
included  in  the  Pg/Wf  and  Pg/Wf  models.  Pole  identifications  are  enumer- 
ated below: 


AP 

p— /Wf  Dynamics 


1 

70% 

85% 

95% 

Association 

Frequency 

Damping 

Frequency 

Damping 

Frequency 

Damping 

Spool  inertia 

-0.716 

— 

- 1.49 

— 

-3.50 

— 

Tailpipe 

| 34. 9 

0.  138 

63.6 

0.  188 

49.7 

0.  004 

dynamics 

( -58.  4 

— 

-31.2 

— 

-60.  4 

mm 

Combustion  time 

( 171.0 

0.  567 

180.0 

0.  541 

240.0 

delay 

( 345.0 

0.062 

— 

— 

— 

E9 

The  lowest- frequency  pole  at  each  operating  point  is  associated  with  spool 
inertia.  These  roots  agree  very  well  with  the  spool  inertia  roots  identified 
from  the  N/W^  and  Pg/W^.  transfer  function  models. 

The  next  two  poles  are  associated  with  gas  dynamics.  Tailpipe  dynamics 
are  represented  by  the  complex  pair  at  about  50  radians  per  second  with  a 
damping  ratio  around  0. 15.  This  representation  roughly  agrees  with  the 
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results  presented  for  the  N/Wp  Pg/W^,  and  Pg/W^.  transfer  functions.  A 
specific  interpretation  of  the  other  pole  is  not  clear. 

The  high-frequency  poles  are  associated  with  the  combustion  time  delay. 
These  poles  are  part  of  Pade  approximations  equivalent  to  the  approximations 
included  in  the  Pg/W^  and  Pg/W^  models.  Thus,  at  70  percent  maximum 
spool  speed  the  delay  is  0.  015  second,  at  85  percent  the  delay  is  0.  014  second 
and  at  95  percent  the  delay  is  0.  011  second. 

Transfer  function  models  for  Tg/W^  are  presented  in  Table  12.  The  dynamic 
states  included  in  these  models  are  identified  below: 


Tg/Wf  Dynamics 


70% 

85% 

95% 

Association 

Damping 

Damping 

Tj.  sensor 

-0.  395 

— 

-0.749 

— 

— 

— 

Spool  inertia 

00 

o> 

CM 

1 

— 

-4.31 

— 

-6.  19 

— 

Gas 

l -94. 0 

i — 

— 

-214.0 

— 

-160.  0 

— 

dynamics 

-360.0 

--- 

--- 

- - “ 

The  lowest-frequency  pole  at  each  operating  point  in  these  models  represents 
the  response  characteristic  of  the  Tg  thermocouple  sensor.  This  pole  is 
missing  from  the  95  percent  maximum  spool  speed  model  because  frequency 
data  were  not  obtained  at  a low  enough  frequency  at  this  operating  point  to 
identify  the  root. 

The  next- lowest-frequency  poles  are  associated  with  spool  inertia.  These 

poles  are  located  at  a higher  frequency  than  the  corresponding  poles  in  the 

A P 

N/W^,  Pg/Wj,  and  p /W^.  'vodels.  This  result  is  consistent  with  the  spool 

3 

inertia  representation  in  the  1>r/Wf  model. 

o i 
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The  high" frequency  poles  have  not  been  identified  specifically  but  are  be- 
lieved to  be  gas  dynamics.  It  should  be  noted  that  the  low-bandpass  quality 
of  the  Tg  sensor  reduces  the  signal-to-noise  ratio  at  high  frequency,  increas- 
ing the  identification  error. 


ENGINE  TRANSFER  FUNCTIONS  FOR  EXHAUST  AREA 

Transfer  function  models  of  spool  speed,  N,  compressor  discharge  pressure, 
P3,  turbine  discharge  pressure,  Pg,  Mach  number  sensor,  Ap/Pg,  and  tur- 
bine discharge  temperature,  Tg,  with  respect  to  exhaust  area,  Ag,  are  pre- 
sented in  Tables  13  through  17.  Representations  of  spool  inertia,  gas  dy- 
namics, and  Tg  sensor  dynamics  in  the  0.  5 to  10- Hz  frequency  range  are 
included  in  the  models.  Important  features  of  the  models  are: 

• Consistency  between  the  exhaust  area  transfer  functions 
and  the  fuel  flow  transfer  functions  is  demonstrated. 

• Spool  inertia  is  consistently  identified  in  the  N/A„,  P„/A„, 

Ap  o o o 

and  „ / A0  transfer  functions. 

P3  8 

• The  Tg/Ag  models  include  the  Tg  thermocouple  sensor 
dynamics. 

Transfer  function  models  of  N/Ag  are  presented  in  Table  13.  Spool  speed  is 
modeled  by  two  lags  in  series,  one  with  a break  frequency  around  2.  0 Hz  and 
the  other  with  a break  frequency  around  10  Hz.  The  first  of  these  two  lags 
represents  spool  inertia.  Specific  location  of  the  roots  associated  with  this 
lag  are  listed  below: 


Poles  Associated  with  Spool  Inertia 


Model 

70% 

85% 

95% 

N/A8 

-0.698 

-2.68 

-3.  38 

N/W, 

-0.714 

-1.  17 

-2.86 
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Spool  inertia  roots  identified  from  the  N/W^.  transfer  functions  discussed  in 
the  previous  section  are  also  shown  in  this  table.  Comparison  of  the  spool 
inertia  roots  identified  from  the  two  models  shows  the  results  to  be  consistent. 
The  models  agree  very  well  for  the  70  percent  and  95  percent  maximum  spool 
speed  operating  points.  Agreement  at  the  85  percent  operating  point  is  not  as 
good. 

The  second  lag,  with  a break  frequency  at  54  radians  per  second,  included  in 
the  95  percent  N/Ag  model  represents  gas  dynamics  in  the  engine  tailpipe. 

This  root  is  also  identified  in  the  95  percent  N/W^  model.  Corresponding 
roots  are  not  identified  in  the  70  percent  and  85  percent  N/Ag  models  because 
the  frequency  data  could  not  be  analyzed  at  high  frequencies. 

Transfer  function  models  of  Pg/Ag  are  presented  in  Table  14.  These  models 
contain  representations  of  spool  inertia  which  are  consistent  with  the  results 
presented  in  the  N/Ag  models.  The  poles  included  in  the  Pg/AQ  transfer 
functions  are  identified  below: 

Pg/Ag  Dynamics 
Association  70% 

Spool  inertia  -0.864 

Tailpipe  dynamics 

The  low-frequency  poles  at  each  operating  point  are  identified  with  spool 
inertia.  These  roots  agree  with  the  corresponding  roots  identified  in  the 
N/Ag  models. 

The  higher- frequency  poles  included  in  the  85  percent  and  95  percent  maxi- 
mum spool  speed  models  are  associated  with  engine  gas  dynamics  in  the  tail- 
pipe. A set  of  roots  at  about  the  same  frequencies  was  also  identified  in  the 
Pg/W^  models  discussed  in  the  previous  section. 


85%  95% 

-2.59  -3.41 

-26.9  -45.8 
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It  should  be  noted  that  the  time  delays  identified  in  the  Pg/W^  models  are  not 
included  in  the  Pg/Ag  models  because  the  Ag  frequency  test  was  only  per- 
formed for  frequencies  below  10  Hz.  In  order  to  measure  the  time  delays 
which  are  included  in  the  models,  frequency  response  data  must  be  mea- 
sured in  the  10  to  100- Hz  frequency  range. 

Pg/Ag  transfer  function  models  are  identified  in  Table  15.  Engine  dynamics 
identified  from  these  models  are  listed  below: 


Pg/Ag  Dynamics 


Association 

70% 

85% 

95% 

Spool  inertia 

— 

-4.  42 

-5.28 

Gas  dynamics 

-33.2 

— 

-97.2 

The  low-frequency  poles  identified  at  each  operating  point  are  associated  with 
spool  inertia.  These  poles  do  not  agree  with  the  spool  inertia  roots  identified 
in  the  N/Ag  and  Pg/Ag  models;  the  roots  presented  here  are  located  at  higher 
frequencies  than  the  corresponding  roots  in  the  N/Ag  and  Pg/Ag  models. 

This  discrepancy  was  also  observed  in  the  comparison  of  spool  inertia  roots 
identified  from  the  Pg/Wf,  Pg/W^,  and  N/W^  models. 

The  high-frequency  poies  included  in  the  70  percent  and  95  percent  Pg/Ag 
models  are  associated  with  engine  gas  dynamics.  Poles  located  at  roughly 
the  same  frequency  are  also  included  in  the  Pg/W^  models. 

Transfer  function  models  of  (AP/Pg)/Ag  are  listed  in  Table  16.  These  models 
contain  representation  of  only  one  state,  spool  inertia.  Locations  of  the 
poles  in  these  models  are  listed  below: 
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Association 


AP 

~=~  /A„  Dynamics 
*3  0 

70%  85%  95% 

Spool  inertia  -1.98  -3.74  -4.90 

These  roots  correspond  more  closely  with  spool  inertia  representations  of 
the  P5/Ag  models  than  they  do  with  the  spool  inertia  representations  included 
in  the  N/Ag  and  Pg/Ag  models. 

Tg/Ag  transfer  function  models  are  presented  in  Table  17.  Consistent 
representation  of  spool  inertia  and  Tg  thermocouple  dynamics  are  contained 
in  these  models.  Specific  root  locations  are  identified  below: 

Tg/Ag  Dynamics 


Association 

70% 

85% 

95% 

Tg  sensor 

-0.  361 

-1.01 

-1.42 

Spool  inertia 

-3.67 

-8. 14 

-13.0 

The  low-frequency  roots  are  associated  with  the  Tg  thermocouple  sensor. 
Location  of  these  poles  agrees  with  the  data  from  the  Tg/W^  transfer  function 
models. 


The  higher- frequency  poles  are  identified  as  spool  inertia.  As  was  the  case 

with  the  fuel  flow  data,  the  spool  inertia  roots  are  at  considerably  higher 

AP 

frequency  in  these  models  than  in  the  N/Ag,  Pg/Ag  or  p-  /Ag  models. 

These  roots  are  also  larger  than  the  spool  inertia  roots  the  Tg/W^  models. 
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ENGINE  TRANSFER  FUNCTIONS  FOR  COMPRESSOR  BLEED 

Transfer  function  models  of  spool  speed,  N,  compressor  discharge  pressure, 
P3,  and  turbine  discharge  temperature  Tg,  with  respect  to  compressor  bleed, 
BLD,  are  presented  in  Tables  18,  19,  and  20.  Spool  inertia,  gas  dynamics, 
and  Tg  sensor  dynamics  in  the  frequency  range  0. 5 to  10  Hz  are  represented 
in  the  models.  In  general,  the  bleed  models  presented  in  this  section  are  not 
as  accurate  as  the  fuel  flow  or  exhaust  area  models  presented  in  the  previous 
sections  due  to  the  relative  insensitivity  of  the  measured  variables  to  the 
bleed  control.  Important  features  of  the  bleed  models  are  summarized  below: 

• Consistency  with  the  fuel  flow  and  exhaust  area  models  is 
demonstrated. 

• Consistent  identification  of  spool  inertia  is  included  in  the  N/BLD 
and  the  Pg/BLD  models. 

• The  Tg/BLD  models  contain  Tg  sensor  characteristics. 

Transfer  function  models  of  N/BLD  are  presented  in  Table  18.  These  models 
contain  representation  of  only  one  state,  spool  inertia.  Location  of  the  poles 
associated  with  this  state  are  listed  below: 

Poles  Associated  With  Spool  Inertia 


Model 

70% 

85% 

95% 

N/BLD 

-0.  936 

-1.  14 

-1.53 

N/Wf 

-0.714 

-1.  17 

-2.86 

N/A8 

-0.698 

-2.68 

-3.38 

Also  listed  in  the  table  are  spool  inertia  roots  identified  from  the  N/W^  and 
N/Ag  models.  Comparison  of  the  data  in  the  table  shows  the  N/BLD  models 
to  be  consistent  with  the  N/W^  and  N/Ag  models.  The  models  agree  very 
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well  at  the  70  percent  maximum  spool  speed  operating  point.  At  the  85  per- 
cent point  the  N/BLD  and  N/W^  model  agree  very  well,  but  the  root  in  the 
N/Ag  root  is  somewhat  high.  Model  agreement  is  not  as  good  at  the  95  per- 
cent operating  point;  the  roots  differ  by  as  much  as  1.  3 radians  per  second. 

Models  of  P3/BLD  are  presented  in  Table  19.  Representations  of  both  spool 
inertia  and  gas  dynamics  are  included  in  these  models.  The  poles  are  iden- 
tified below: 

P3/BLD  Dynamics 

Association  70%  85%  95% 

Spool  inertia  -1,06  -1,98  -1.96 

Gas  dynamics  -80.1  -40.  1 -76.3 

Neither  spool  inertia  nor  gas  dynamics  are  represented  as  accurately  in  these 
data  as  in  the  data  corresponding  to  the  Pg/Wf  and  Pg/Ag  models.  The  spool 
inertia  root  associated  with  the  85  percent  Pg/BLD  model  appears  to  be  too 
large.  This  is  also  the  case  with  the  gas  dynamics  roots  identified  in  the 
70  percent  and  95  percent  models. 

Tg/BLD  transfer  function  models  are  listed  in  Table  20.  The  poles  which 
have  been  identified  in  these  models  are  listed  below: 

Te/BLD  Dynamics 
0 

Association 

Te  sensor 
o 

Thermocouple  time  constants  identified  in  these  models  agree  with  the 
results  presented  for  the  Tg/W^.  and  Tg/Ag  transfer  functions.  The  high- 
frequency  root  identified  in  the  85  percent  and  95  percent  models  is  probably 
incorrectly  located. 


70%  85%  95% 

0.351  -1.15  -1.88 

-159.0  -416.0 
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ENGINE  TRANSFER  FUNCTIONS  FOR  INLET  GUIDE  VANE 


Transfer  function  models  of  spool  speed,  N,  compressor  discharge  pres- 
sure, Pg,  and  turbine  discharge  temperature  Tg,  with  respect  to  inlet  guide 
vane,  IGV,  are  presented  in  Tables  21,  22,  and  23.  These  models  contain 
representations  of  spool  inertia  and  Tg  thermocouple  dynamics  in  the  0. 5 to 
10- Hz  frequency  range. 

Models  are  presented  for  only  the  85  percent  and  95  percent  maximum  spool 
speed  operating  points.  At  the  70  percent  operating  point  the  effect  of  the 
IGV  control  on  the  engine  is  insignificant.  IGV  frequency  response  data  was 
not  measured  at  this  operating  point  . 

The  IGV  models  presented  here  are  not  as  accurate  as  the  fuel  flow  and 
exhaust  area  models  presented  earlier  in  this  section  because  the  engine  is 
not  as  sensitive  to  the  IGV  control  as  it  is  to  the  fuel  flow  and  exhaust  area 
controls.  Thus,  some  of  the  poles  identified  in  the  IGV  models  are  more 
subject  to  error  due  to  the  extremely  low  signal-to-noise  ratio  which  charac- 
terizes the  IGV  frequency  responses.  Fortunately,  the  extraneous  dynamics 
are  easily  identified  by  comparing  the  IGV  models  with  the  fuel  flow,  exhaust 
area,  and  compressor  bleed  models  previously  discussed. 

The  significant  features  of  the  IGV  models  are: 

• Lbw-frequency  dynamics  (below  2.  0 Hz)  are  consistently 
identified  between  the  IGV  models  and  the  W^,  Ag,  and  BLD 
models. 

• Spool  inertia  states  are  included  in  the  N/IGY,  Pg/IGV,  and 
Tg/IGV  models. 

• The  T5/IGV  models  contain  Tg  thermocouple  dynamics. 
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Transfer  function  models  of  N/IGV  are  presented  in  Table  21.  These  models 
contain  spool  inertia  dynamics  which  are  consistent  with  the  N/Wf,  N/Ag, 
and  N/BLD  models.  Spool  inertia  roots  of  all  four  models  are  listed  below: 


Poles  Associated  With  Spool  Inertia 


Model 

70% 

85% 

95% 

N/IGV 

— 

-1.  81 

-2.70 

N/Wf 

-0.714 

-1.  17 

-2.86 

n/a8 

-0.  698 

-2.68 

-3.38 

N/BLD 

-0.  936 

-1.  14 

-1.53 

Agreement  between  the  N/IGV  models  and  the  other  three  models  is 
demonstrated. 

The  N/IGV  model  at  the  95  percent  maximum  spool  speed  operating  point  also 
contains  a second  pole  located  at  19  radians  per  second.  Comparison  of  this 
pole  with  the  poles  identified  in  the  N/W^,  N/Ag,  and  N/BLD  models  suggests 
that  the  pole  is  incorrectly  identified.  A second  pole  does  exist  in  the  N/Wf 
and  N/Ag  models  for  the  95  percent  operating  point,  but  it  is  at  a higher 
frequency,  around  60  radians  per  second. 

Models  of  Pg/IGV  are  presented  in  Table  22.  The  only  clearly  identifiable 
dynamics  in  these  models  are  the  representations  of  spool  inertia.  All  of  the 
poles  for  these  models  are  listed  below: 

Pg/IGV  Dynamics 
Association  85% 

Spool  inertia  -2.  01 

-17.  8 
-1570.  0 


95% 

-3.26 

-400.  0 
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The  poles  associated  with  spool  inertia  agree  very  well  with  the  spool  inertia 
roots  identified  in  the  N/IGV  models.  This  relationship  is  consistent  with  the 
results  obtained  for  the  fuel  flow,  exhaust  area,  and  compressor  bleed  models. 

The  85  percent  maximum  spool  speed  model  also  contains  a second  pole 
located  at  17.  8 radians  per  second.  This  pole  is  believed  to  be  associated 
with  tailpipe  gas  dynamics,  although  it  is  located  at  a lower  frequency  than 
the  corresponding  poles  in  the  85  percent  P^/W^,  P^/Ag,  and  P^/BLD 
models. 

The  high-frequency  poles  included  in  both  the  85  percent  and  95  percent  maxi- 
mum spool  speed  Pg/IGV  models  are  extraneous.  They  are  located  at  too 
high  a frequency  *o  be  correctly  measured  on  frequency  responses  in  the  0.  5 
to  10- Hz  frequency  range. 

Tg/IGV  frequency  response  models  are  presented  in  Table  23.  These  models 
contain  representations  of  spool  inertia  and  Tg  thermocouple  dynamics  which 
agree  with  the  results  previously  discussed  for  the  Tg/Wf,  Tg/Ag,  and  Tg/ 

BLD  models.  In  addition,  the  Tg/IGV  models  also  contain  a set  of  poles  in 
the  3 to  9-Hz  frequency  rang0  which  is  tentatively  associated  with  gas  dynamics. 
The  poles  included  in  these  models  are  identified  below: 

Tg/IGV  Dynamics 


Association 

85% 

95% 

T,-  thermocouple 

-0.390 

-1.90 

Spool  inertia 

-4.59 

CO 

o 

1 

Gas  dynamics 

-17.9 

-58.5 
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Table  4.  Fuel  Valve  Transfer  Functions 


N/N 
' max 

Fuel  Nozzle  Pressure/Fuel  Flow  Command, 

70% 

1.0 

Hr-‘l 

[| 

lifs),+aTtlr*15+1][r 

b*‘l 

85% 

1.0 

[!«of  * + '] 

[tii)’  + aT!iPs  + 1l 

(1*1* 

[| 

S l2  2(0. 206)S  ,1  r 

44.9/  + 44.9 

95% 

1.0 

[Mil**  “=*#“♦>' 

IflAf 

[1 

it 
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Table  5.  Exhaust  Actuator  Transfer  Functions 


Table  6.  Bleed  Actuator  Transfer  Functions 
N/Nmax  Bleed  Position /Bleed  Command,  BLD/uBLD 
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Table  7.  Inlet  Guide  Vane  Actuator  Transfer  Functions 


Nm  Qv 

max 

IGV  Position/IGV  Command,  IGV/uIQV 

85% 

‘•“[ifa-1! 

[im  + 1] 

95% 

1,0  [l89  * *1 

[ro + •] 

Table  8.  Fuel  Flow  Transfer  Functions 
Spool  Speed/ Fuel  Flow 


N'Nmax 

Spool  Speed/ Fuel  Flow,  N/Wf 

70% 

12- 65  [its  - 

■ s + ii  r s . . + ji 

IOT4  + XJ  L4O  XJ 

85% 

6*  36  1.82.  3 “ *] 

.U7  + l]  left + x] 

95% 

2-86  [ifs-1] 

1 
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Table  10.  Fuel  Flow  Transfer  Functions — Turbine  Di 
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Table  11.  Fuel  Flow  Transfer  Functions- -Mach  No, 
Sensor /Fuel  Flow 


90 


Table  12.  Fuel  Flow  Transfer  Functions --Turbine 
Discharge  Temperature /Fuel  Flow 


N'Nmax 

Turbine  Discharge  Temperature /Fuel  Flow,  Tg/W^ 

70% 

-0.  328[0  450  -l]  [122  *] 

[o.395  + 1][2.98  ' ']  [94.4  hl] 

85% 

«-  *«  * *]  [ifs  - -K 1] 

95% 

[09  + *]  [i§3 

Table  13.  Exhaust  Area  Transfer  Functions --Spool 
Speed/Exhaust  Area 


N/[W 

Spool  Speed /Exhaust  Area,  N/Ag 

70% 

3.84 

[o.m + 

85% 

2.99 

[2.68  + 1 ] 

95% 

11.0 

[jTSFl  + >]  [5CT+  *] 
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Table  14.  Exhaust  Area  Transfer  Functions — Compressor 
Discharge  Pressure  /Exhaust  Area 


N/N 

max 

Compressor  Discharge  Pres  sure /Exhaust  Area,  P3/AQ 

70% 

0.0601 

liO-1] 

I s 

Xll 

10.864  ' XJ 

85% 

0.0512 

L 13.  9 ' A [*173  _1] 

[ 2^59  + *]  [26.9  +1] 

95% 

0 . 044  5 

[or+1][40+1J 

Table  15.  Exhaust  Area  Transfer  Functions --Turbine 
Discharge  Pressure /Exhaust  Area 


n/n 

max 

Turbine  Discharge  Pressure /Exhaust  Area,  Pg/Ag 

70% 

-0-0175  [o.?19  + 1] 

[24.5  + *] 

[ 2.^27  + -1]  [33^2  + X] 

85% 

.0.°363  [r|5  + l] 

[o5+1] 

V 

95% 

-0'100  (o^Ti  + 1] 

[Si"+  *]  [sr? + 1 
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Table  16.  Exhaust  Area  Transfer  Functions --Mach  No. 
Sensor /Exhaust  Area 


Table  17.  Exhaust  Area  Transfer  Functions --Turbine 
Discharge  Temperature /Exhaust  Area 

N/Nmax  Turbine  Discharge  Temperature /Exhaust  Area,  Tg/Ag 


Table  18.  Bleed  Transfer  Functions — Spool 
Speed/ Bleed  Position 


N/N 

max 

Spool  Speed /Bleed  Position,  N/BLD 

70% 

-402  [tW-1] 

[ 0.936  ' 1) 

85% 

'505[ir?  - ‘] 

[iV>] 

95% 

134[5S27-‘] 

1 

rH 

+ 

eo 

rH 

Table  19.  Bleed  Transfer  Functions--Compressor 
Discharge  Pressure/Bleed  Position 


N/N 

max 

Compressor  Discharge  Pres  sure /Bleed  Position,  Pg/BLD 

70% 

-4.63 

85% 

”**'  85  [ 5.  94  + *]  [81.  6 + *] 

[&*  + '}[& T+‘] 

95% 

-4'54[  1S2  +1J[77S.  -1! 

f S . 1 f S .1 

Ll.  96  1J  L76TT+  *] 
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Table  20.  Bleed  Transfer  Functions --Turbine  Discharge 
Temperature /Eleed  Position 


N'Nma* 

Turbine  Discharge  Temperature /Bleed  Position,  Tg/BLD 

70% 

83*4  lor  ■ *] 

[os!  + *] 

85% 

Lrfe+>] 

95% 

19-' 

Table  21.  Inlet  Guide  Vane  Transfer  Functions - 
Spool  Speed /IGV  Position 


N/N 

max 

Spool  Speed/IGV  Position,  N/IGV 

85% 

457  LsV1. 

95% 

631 

.50  -1. 

[tt+1] 

■msu 
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Table  22.  Inlet  Guide  Vane  Transfer  Functions  --Compressor 
Discharge  Pres  sure /IGV  Position 


N/N 
' max 

Compressor  Discharge  Pressure/IGV  Position,  Pg/IGV 

85% 

°*  280  Lo.  079  “ *]  [29.  6 + *]  [48. 1 - !] 

[oi  + ;]  [n75  + *1  [iFT5  + *] 

95% 

°-  886  [o.?75-1][3iS4-1] 

Table  23.  Inlet  Guide  Vane  Transfer  Functions --Turbine 
Discharge  Temperature /IGV  Position 


N/N 

m&x 

Turbine  Discharge  Temperature /IGV  Position,  Tg/IGV 

85% 

"18-3  [ 6:376  -1 

1 [ iOT  ' 1 

H* h" 

.] 

[o.390  * l]  [ 

jV*][ 

95% 

wm 

Olrai 

] 

[05+1][7S08  + 1] 

[«£*+»] 
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APPENDIX  A 

DESCRIPTION  OF  INSTRUMENTATION 


This  appendix  describes  the  equipment  used  to  measure  the  frequency  response 
of  the  engine.  A block  diagram  of  the  experimental  setup  is  shown  in  Figure 
A-l.  The  equipment  includes  engine  sensors,  interface  electronics,  an  ana- 
log computer,  an  FM  tape  recorder,  and  a two- channel  servoanalyzer. 

The  BAFCO  servoanalyzer  produces  a sinusoidal  test  signal  which  drives  the 
engine  through  one  of  the  actuators.  A system  of  actuator  and  engine  sensors 
measures  the  responses  of  interest.  Signals  from  the  sensors  are  routed  to 
their  destinations  through  the  interface  electronics  which  serve  as  a link 
between  the  engine  and  the  frequency  analysis  equipment.  The  sensor  signals 
are  also  scaled  in  the  interface  electronics.  All  of  the  sensor  signals  are 
sent  to  the  analog  computer  from  the  interface  element.  The  DC  portion  of 
the  signals  is  removed  in  the  computer  and  the  remaining  AC  component  is 
amplified.  The  amplified  AC  components  are  relayed  to  the  tape  recorder 
where  they  are  recorded  along  with  three  reference  signals  produced  by  the 
BAFCO  servoanalyzer. 

In  addition  to  the  signals  which  are  recorded  on  tape,  one  input-output  pair, 
spool  speed /actuator  command,  is  analyzed  directly  to  monitor  the  frequency 
response  test.  These  two  signals  are  also  recorded  on  tape  to  facilitate 
checking  out  the  instrumentation  in  the  tape  playback  mode.  In  the  playback 
mode  the  BAFCO  servoanalyzer  is  connected  directly  to  the  tape  recorder  so 
that  the  recorded  signals  can  be  analyzed  one  pair  at  a time.  Instrumentation 
is  discussed  in  more  detail  in  the  following  paragraphs. 
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ENGINE  ACTUATORS 


All  four  engine  actuators  (fuel  valve,  exhaust  area,  compressor  bleeds,  and 
inlet  guide  vanes)  are  position  control  mechanisms  of  the  type  represented 
by  the  following  simplified  block  diagram: 


POSITION  

ACTUATOR 

ACTUATOR 

REQUEST 

J 

\m  \ 

POSITION 

The  actuators  are  open-loop  devices  with  loops  closed  in  the  interface  elec- 
tronics (except  exhaust  area).  Inputs  to  the  actuators  are  position  requests 
which  are  adjusted  through  the  interface  electronics.  The  DC  content  of  the 
request  signals  is  manually  selected  through  four  potentiometers  (one  for 
each  actuator)  located  on  the  operator  console.  A summing  junction  is  also 
pr  ovided  so  that  a sinusoidal  perturbation  signal  from  the  BAFCO  servo- 
analyzer  can  be  added  directly  to  the  DC  position  requests.  The  individual 
actuators  are  discussed  in  the  following  paragraphs. 


Fuel  Valve 


A schematic  diagram  of  the  fuel  system  is  shown  in  the  upper  half  of  Figure 
A -2.  The  system  is  represented  by  two  components,  a fuel  valve  and  a fuel 
line  which  delivers  fuel  from  the  valve  to  the  spray  nozzles.  Dynamic  char- 
acteristics of  these  components  are  considered  in  the  following  paragraphs. 
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The  fuel  valve  contains  two  stages,  a metering  valve  and  a bypass  valve. 
Position  feedback  is  provided  around  both  stages.  Valve  characteristics 
established  from  bench  tests  of  the  actuator  are  identified  as: 

Natural  Frequency  Damping  Ratio 
Metering  valve  30  Hz  0. 80 

Bypass  valve  20  Hz  0.20 

i 

The  fuel  valve  is  connected  to  the  engine  with  approximately  8 feet  of  fuel 
line  varying  in  diameter  from  0.  5 inch  to  0.  75  inch.  Fluid  inertia  and  com- 
pressibility effects  were  analyzed  to  estimate  the  dynamic  characteristics  of 
this  line.  The  results  show  that  the  fuel  line  does  not  contain  any  dynamics 
in  the  frequency  range  below  50  Hz. 

Fuel  flow  dynamics  are  coupled  to  engine  dynamics  by  the  effect  of  Pg.  In 
the  experimental  tests,  fuel  flow  was  obtained  by  measuring  the  pressure 
differential  across  the  nozzle  openings,  since  the  bandwidth  of  the  flowmeter 
was  much  less  than  the  metering  valve.  Nozzle  pressure  differential  is 
defined  as 

APf„  ■ pfn  - p3  <A-» 

where 

AP^n  = Fuel  nozzle  pressure  differential 
P^n  = Pressure  in  the  fuel  line 
P3  = Compressor  discharge  pressure 

Fuel  line  pressure  was  measured  with  a wall  static  pressure  tap  in  the  fuel 
line  located  about  6 inches  downstream  of  the  fuel  valve.  Compressor 
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discharge  pressure  was  measured  with  a wall  static  tap  located  just  aft  of 
the  compressor  outlet  guide  vanes.  These  sensors  are  described  in  detail 
in  the  following  sections.  Neither  sensor  contains  any  dynamics  in  the  fre- 
quency range  below  100  H*. 


The  outputs  of  the  two  sensors  were  algebraically  combined  in  the  A PL  ana- 
log computer  to  obtain  AP^.  Nozzle  pressure  differential  is  related  to  actual 
fuel  flow  by  the  relation 


CdAN 


(A-2) 


o 

where  is  nozzle  area,  gc  is  the  gravitational  constant  (gc  / 32.2  ft/sec  ) 
and  p is  density  of  the  fuel.  The  discharge  coefficient,  C^,  in  this  expres- 
sion was  determined  from  steady-state  fuel-flow  calibration  data. 


The  dynamic  characteristics  of  both  the  fuel  valve  and  the  fuel  line  were 
measured  together  in  the  frequency  response  tests.  A block  diagram  of  the 
procedure  is  shown  in  the  lower  half  of  Figure  A-2. 


Since  Pg  was  used  in  the  fuel  flow  measurement,  the  transfer  function  models 
of  the  fuel  system  contain  some  engine  dynamics  in  addition  to  actuator  and 
line  dynamics.  The  transfer  function  for  the  fuel  system  is 


APfn<s> 

Uf(s) 


N^sJCDJsJDJs)  - N„(s)N,(s)] 
~ D1(s)D2<s>D3<s> 


(A -3) 


where 

N^(s)  represents  the  numerator  dynamics  of  the  fuel  valve 
D^(s)  represents  the  denominator  dynamics  of  the  fuel  valve 
Ngis),  Dgis)  represent  fuel  line  dynamics 
Ng(s),  Dg(s)  represent  engine  dynamics  in  the  Pg  response 
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Thus,  the  poles  associated  with  the  fuel  system  transfer  function  are  the 
product  of  the  poles  associated  with  the  three  individual  components:  fuel 

m 

valve,  fuel  line,  and  Pg  response.  This  explains  why  some  engine  dynamics 
appear  in  the  fuel  system  transfer  functions  presented  in  Section  III. 

The  fuel  line  dynamics  are  negligible  below  50  Hz;  therefore,  the  transfer 
functions  are  dominated  by  metering  valve  and  engine  dynamics: 


4P,n<s> 

uf(s) 


N1(s)CDg(s)  - N3(s)] 
DjUJDgfs) 


(A -4) 


if  frequency  (u)  £ 50  Hz. 

Similarly,  the  transfer  function  models  of  engine  responses  with  respect  to 
fuel  flow  contain  some  fuel  system  dynamics  in  addition  to  engine  dynamics. 
For  example,  the  transfer  function  representation  for  Pg  is 


P3(S) 


N2(s)N3(s) 

APfn(s)  = D2(s)D3(s)  - N2(s)N3(s) 


(A-5) 


Thus,  fuel  line  dynamics  described  by  N2(s)  and  D2(s)  are  included  in  this 
model.  Even  though  the  line  dynamics  can  be  neglected,  the  transfer  function 
models  obtained  are  still  not  the  desired  result:  N3(s)/Dg(s).  That  is. 


p3(s> 

AP.  (s) 
in 


n3(s) 

D3(s)  - Ng(s) 


(A-6) 


if  u £ 50  Hz. 


The  poles  associated  with  Pg  dynamics,  Dg(s),  are  shifted  by  the  numerator 
dynamics,  Ng(s). 
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Although  the  individual  transfer  function  models  obtained,  AP^/u^  and 
P3^Pfn'  n0*  ^ave  the  desired  form»  the  model  for  P3 /uf  obtained  by 
multiplying  AP^/u^  and  P3 together  is  valid.  That  is, 

P3  P3  *Pfn 

uf  APfn  uf 


N2(s)N3(s)  N1(s)[D2(s)D3(s)  - N2(s)N3(s)] 

D2(s)D3(s)  - N2(s)N3(s)  * D^sjD^slD^s) 

N1(s)N2(s)N3(s) 

D1(s)D2(s)Dg(s) 


It  is  this  model,  Pg/u^.,  not  the  individual  models,  Pg/AP^n  and  AP^/u^, 
which  is  used  for  controller  synthesis.  However,  one  should  keep  in  mind 
that  the  engine  responses  to  fuel  flow  have  their  high-frequency  poles  shifted 
by  pressure  dynamics. 


Exhaust  Actuator 

The  exhaust  actuator  is  a mechanical  clutch-brake  device  geared  to  the  rotor 
shaft.  The  main  body  of  the  actuator  is  mounted  on  a pad  near  the  front  of 
the  engine.  A screw  linkage  connects  the  actuating  mechanism  with  the 
nozzle  apparatus  at  the  rear  of  the  engine. 

Input  to  the  actuator  is  requested  nozzle  area.  A position  feedback  loop 
around  the  actuator  is  built  into  the  hardware. 

Nozzle  area  is  measured  by  a linear  potentiometer  mounted  on  the  actuator 
assembly  which  senses  actuator  movement.  A steady-state  calibration  was 
made  to  convert  the  potentiometer  readings  in  voltage  to  nozzle  area  units. 
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Compressor  Bleed  Actuator 


The  bleed  actuator  is  a hydraulically  powered  valve  which  controls  the 
amount  of  air  vented  out  of  the  third,  fourth,  and  fifth  stages  of  the  compres- 
sor. Pressurized  fuel  delivered  by  an  auxiliary  pump  is  used  to  drive  the 
actuator. 

Bleed  area  is  measured  by  a linear  potentiometer  affixed  to  the  actuator 
linkages.  The  potentiometer  measures  actuator  movement.  A steady-state 
calibration  was  performed  to  convert  the  potentiometer  readings  from  vol- 
tage units  to  nondimensional  units.  In  the  latter  system  of  units,  1 
corresponds  to  fully  open  bleeds  and  0 corresponds  to  closed  bleeda 


Inlet  Guide  Vane  Actuator 

The  IGV  actuator  is  a hydraulic  actuator  which  regulates  the  incidence  angle 
of  the  compressor  inlet  guide  vanes.  Fluid  pressure  to  drive  the  actuator  is 
provided  by  the  same  pump  which  drives  the  bleed  actuator. 

IGV  angle  is  measured  by  a linear  potentiometer  which  measures  actuator 
displacement.  The  potentiometer  readings  are  converted  to  a nondimensional 
system  of  units,  with  1 corresponding  to  the  nominal  low-speed  IGV  setting 
and  0 corresponding  to  the  nominal  high-speed  IGV  setting. 


ENGINE  SENSORS 
Spool  Speed  Sensor 

Spool  speed  was  measured  with  an  electronic  sensor  which  measures  elapsed 
time  per  revolution  of  the  rotor  shaft.  The  sensor  is  composed  of  two  parts, 
a gear  wheel  with  magnetic  teeth  and  a 2. 4 7 -MHz  oscillator.  The  gear  wheel 
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is  connected  to  the  rotor  shaft  to  turn  at  a fraction  of  the  rotor  speed.  A 
magnetic  pickoff  attached  to  the  engine  and  positioned  over  the  gear  wheel 
counts  revolutions  of  the  gear  wheel  by  detecting  magnetic  teeth  on  the  gear. 
Rotor  speed  is  determined  by  counting  the  number  of  oscillator  cycles  per 
revolution  of  the  gear  wheel  and  performing  the  simple  calculation: 

N = oscillator  cycles /revolution  x 2'47  MHz/sec  x gear  ratio 

This  sensor  does  not  contain  any  dynamics  in  the  frequency  range  tested, 

0.  04  to  100  Hz. 


Compressor  Discharge  Pressure  Sensor 

A static  pressure  tap/ pressure  transducer  sensor  was  used  to  measure  com- 
pressor discharge  pressure.  The  static  tap  is  embedded  in  the  wall  of  the 
engine  about  2 inches  behind  the  compressor  outlet  guide  vanes.  A piece  of 
tubing  0.  25-inch  in  diameter  and  about  4.  5 inches  long  connects  the  tap  to 
the  pressure  transducer.  The  fundamental  resonance  frequency  of  the  con- 
nection line  is  about  650  Hz,  and  the  response  of  the  pressure  transducer  is 
flat  to  beyond  10,  000  Hz. 


Turbine  Discharge  Pressure  Sensor 

Turbine  discharge  pressure  was  measured  with  a system  of  five  total  pres- 
sure probes  spread  around  the  circumference  of  the  engine  behind  the  turbine 
outlet.  Readings  of  the  individual  probes  are  averaged  to  determine  the 
turbine  discharge  pressure.  The  five  probes  are  connected  to  a 0. 125-inch- 
diameter  pressure  manifold  which  encircles  the  engine  and  is  connected  to  a 
pressure  transducer.  Thus,  the  lengths  of  the  lines  between  the  pressure 
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probes  and  the  pressure  transducer  vary  between  1 foot  and  3 feet.  The 
fundamental  resonant  frequency  in  the  connection  lines  is  estimated  at  200 
Hz.  The  response  of  the  pressure  transducer  is  flat  to  beyond  10,  000  Hz. 


Mach  Number  Sensor 

f*T  ” PS 

A special  sensor  built  by  Bendix  was  used  to  measure  — p at  the  com- 

rT 

pressor  exit.  The  main  body  of  the  sensor  is  connected  to  static  and  total 
pressure  taps  located  at  the  compressor  exit  by  lines  which  are  0.  25-inch 
in  diameter  and  about  a foot  long.  The  resonant  frequency  of  these  lines  is 
estimated  at  400  Hz.  A bench  test  of  the  sensor  showed  its  response  to  be 
flat  out  to  about  20  Hz. 


Turbine  Discharge  Temperature  Sensor 

Turbine  discharge  temperature  is  measured  by  a system  of  19  thermocouples 
spread  around  the  circumference  of  the  engine  about  8 inches  behind  the  tur- 
bine exit.  The  individual  thermocouples  are  connected  in  parallel  so  that  sin 
average  temperature  reading  is  obtained. 


Fuel  Nozzle  Pressure 

A static  pressure  tap  embedded  in  the  wall  of  the  line  connecting  the  fuel 
valve  with  the  fuel  nozzles  is  used  to  measure  fuel  nozzle  pressure.  The  tap 
is  positioned  about  6 inches  downstream  of  the  metering  valve.  A 4 -inch 
length  of  0.  2 5 -inch-diameter  tubing  connects  the  pressure  tap  with  a pressure 
transducer.  The  fundamental  resonance  frequency  of  this  tubing  is  well  be- 
yond the  maximum  frequency  tested,  100  Hz.  The  frequency  response  of  the 
transducer  is  flat  to  beyond  10,  000  Hz. 
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INTERFACE  ELECTRONICS 


The  interface  electronics  package  serves  as  a link  between  the  test  instru- 
ments and  the  engine.  The  package  consists  mainly  of  amplifier  circuits 
which  perform  two  functions:  (1)  actuator  commands  from  the  test  equipment 
are  scaled  to  be  compatible  with  the  actuator  hardware;  and  (2)  sensor  out- 
puts from  the  engine  are  scaled  to  fall  in  the  voltage  range  0 to  5.  0 volts. 

This  equipment  does  not  affect  the  measurement  of  engine  dynamics  in  the 
frequency  range  tested,  0.  04  to  100  Hz. 


BAFCO  SERVOANALYZER 

A BAFCO  servoanalyzer  was  used  to  measure  the  frequency  responses  of  the 
engine.  This  instrument  is  capable  of  obtaining  Bode  frequency  response 
■Acts  (amplitude  ratio  and  phase  shift  versus  frequency)  either  directly  from 
the  engine  or  from  tape-recorded  data.  Operation  of  the  analyzer  is  de- 
cribed  in  the  following  paragraphs. 

The  normal  or  on-line  operation  of  the  servoanalyzer  is  shown  in  the  upper 
half  of  Figure  A-3.  In  this  mode  of  operation  the  servoanalyzer  measures 
the  frequency  response  between  a pair  of  engine  sensor  outputs,  labelled 
Ej(t)  and  Eg(t).  The  servoanalyzer  produces  a sinusoidal  test  signal  with 
time -dependent  frequency  (the  frequency  of  the  test  signal  is  swept  either 
linearly  or  logarithmically  with  time)  which  drives  the  engine  through  one  of 
the  actuators.  Responses  from  the  two  sensors  are  fed  into  the  servoanalyzer, 

E2 

where  amplitude  ratio,  | -g—  | , and  phase  shift,  ~ • of  x - y plotters 

which  record  the  data  as  a function  of  frequency. 

In  addition  to  on-line  operation,  the  servoanalyzer  can  also  be  used  to  obtain 
Bode  frequency  data  from  tape-recorded  signals.  This  operating  mode  is 
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called  the  tape  mode  and  is  shown  in  the  lower  half  of  Figure  A-3.  Opera- 
tion in  this  mode  is  identical  to  on-line  operation  with  the  following  excep- 
tions : 

• In  the  tape  mode  the  signals  to  be  analyzed,  E^t)  and  E2(t), 
are  obtained  from  playback  of  the  data  tape  instead  of  directly 
from  the  engine  sensors. 

• To  maintain  the  precise  relationship  to  the  drive  signal, 
three  additional  reference  signals  are  needed  to  synchronize 
the  BAFCO  servoanalyzer  computations.  These  three  signals 
are  recorded  on  the  data  tape  at  the  same  time  as  the  sensor 
signals,  E^(t)  and  E2(t),  are  recorded. 

Both  the  on-line  and  tape  modes  were  used  to  measure  the  engine  frequency 
responses.  During  the  actual  engine  tests,  all  sensor  signals  were  recorded 
on  magnetic  tape,  and,  in  addition,  the  signals  from  one  sensor  pair  were 
analyzed  on-line  to  monitor  the  results.  Frequency  responses  from  the 
other  signals  were  obtained  later  by  playing  back  the  data  tape  with  the 
analyzer  operating  in  the  tape  mode. 

The  computations  performed  in  the  servoanalyzer  are  shown  in  the  block 
diagram  of  Figure  A-4.  First,  the  fundamental  components  of  the  incoming 
sensor  signals,  E2(t)  and  E^t),  are  identified  by  Fourier  filters.  These 
filters  remove  noise  and  harmonics  from  the  sensor  signals  by  dividing  each 
signal  into  an  in-phase  component,  I,  and  a quadrature  (out-of-phase)  com- 
ponent, Q.  Then,  a coordinate  transformation  is  performed  on  the  com- 
ponents to  put  them  in  a form  better  suited  for  computing  the  amplitude  ratio 
and  phase  shift.  Finally,  the  Bode  plot  parameters  are  computed  in  the 

amplitude  ratio,  | jjr~|  , and  phase  shift,  02  - 0^ , computer.  A detailed 

description  of  these  computations  is  included  in  the  following  paragraphs. 
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Ideally,  the  Incoming  sensor  signals,  E2(t)  and  E^t),  would  both  be  sine 
waves  at  the  test  signal  frequency,  but  phase  shifted  and  of  different  ampli- 
tude. In  this  case,  the  Fourier  filters  shown  in  Figure  A-4  would  not  be 
necessary,  as  the  amplitude  ratio  and  phase  shift  could  be  readily  identified 
from  the  raw  signals.  However,  in  practice  the  sensor  signals  are  rich  in 
harmonic  content  and  contain  excessive  noise  which  must  be  removed  before 
the  amplitude  ratio  and  phase  shift  data  can  be  identified.  Therefore,  the 
sensor  signals  are  processed  by  Fourier  filters  which  automatically  reject 
harmonic  components  and  noise  from  the  signals. 

A block  diagram  of  a Fourier  filter  is  shown  in  Figure  A- 5.  The  sensor 
signal  is  analyzed  in  the  filter  by  multiplying  it  by  a sine  and  cosine  refer- 
ence and  integrating  each  product  over  a whole  number  of  cycles.  That  is, 
the  following  integrals  are  computed: 

NT 

E(t)  sin  Wjr  dt 
NT 

E(t)  cos  Wjt  dt 

where 

N = Number  of  cycles  over  which  integration  is  performed 
(a  whole  number) 

T = Time  per  cycle  = 2tt/u)j 

Wj  = Test  frequency  in  radians  per  second 

It  should  be  noted  that  the  test  signal  is  assumed  to  be  maintained  at  a con- 
stant frequency  in  this  development,  whereas  the  frequency  is  actually  time- 
dependent  in  practice. 


108 


The  incoming  sensor  signal,  E(t),  can  be 
the  test  frequency,  , : 

E(t)  = Aq  + Aj  sin  Ujt  + . . . An 

+ B-  cos  u,t  + . . . B„ 
li  n 

+ (noise) 


expanded  in  a Fourier  series  in 
sin  nidjt 

cos  nujt  (A-8) 


If  we  ignore  the  noise  term  in  this  expression  and  substitute  the  resulting 
E(t)  into  the  above  integrals,  we  obtain: 


JL 

NT 


NT 


NT 


I 


E(t)  sin  Ujt  dt 


NT 


I 


E(t)  cos  Wjt  dt 


(A-9) 


since 


1 

NT 


1 

NT 


±_ 

NT 


NT 

sin  Wjt  • sin  nu^t  dt 
NT 

cos  Wjt  • cos  nwjt  dt 


NT 


sin  Wjt  • cos  nwjt  dt 


1 for  n = 1 
0 for  n # 1 

1 for  n = 1 
0 for  n ^ 1 

0 for  all  n 
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NT 

cos  Ujt  ♦ sin  nwjt  dt  = 0 for  all  n 

Thus,  the  Fourier  filter  rejects  the  harmonic  content  of  E(t)  and  identifies 
the  fundamental  components  of  the  signal,  Aj  and  These  constants  are 
then  multiplied  by  sin  u^t  and  cos  Ujt  to  give  the  in-phrse  and  quadrature 
components  of  E(t): 

In-phase  = Aj  sin  Ujt 
Quadrature  = cos  w^t 

An  analysis  of  the  noise  rejection  capability  of  the  Fourier  filter  (presented 
in  Reference  6)  indicates  that  the  filter  also  effectively  removes  noise  from 
the  sensor  signs  1,  E(t).  The  results  presented  in  the  reference  show  that 
noise  at  frequencies  other  than  the  test  frequency,  Uj , is  filter  out,  and  the 
noise  contribution  which  does  get  through  the  filter  (noise  at  frequencies  near 
Uj)  can  be  identified  from  visual  inspection  of  the  results. 

To  get  the  in-phase  and  quadrature  components  in  a form  better  suited  for 
the  computation  of  amplitude  ratio  and  phase  shift,  they  are  converted  to 

Aj  sin  Ujt  + Bj  cos  Ujt  = E sin  (u^T  + 0 ) (A-10) 

by  the  transformation 


E = 


(Ax2  +B12)1/2 
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9 

Then,  in  the  final  step,  the  amplitude  ratio,  ( — | , and  phase  shift, 

<t>2  - * between  the  Channel  2 signal,  Eg  sin  (u^t  + 0g),  ^ the  Channel  1 

signal,  Ej  sin  (Ujt  + are  computed. 

TAPE  RECORDER  - ANALOG  COMPUTER 

An  Ampex  tape  recorder  and  the  APL  analog  computer  were  used  to  record 
sensor  signals  during  the  engine  frequency  response  tests.  A block  diagram 
of  the  sensor  signal  recording  process  is  shown  in  Figure  A-6.  The  analog 
computer  was  used  to  remove  the  DC  bias  on  the  sensor  signals  and  to 
amplify  the  AC  components  before  the  signals  were  recorded  on  magnetic 
tape. 

Thirteen  signals  were  recorded  on  the  data  tape.  They  are  identified  below: 

Channel 

1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 


Description 


DC  level 

Triangle  wave  BAFCO  servoanalyzer 

riangie  wave  reference  signals 

Square  wave 

Actuator  input  (uf,  uA  , uBLD,  uIGV ) 

8 

Actuator  position  (W^,  Ag,  BLD,  IGV) 
Rotor  speed  (N) 

Compressor  discharge  pressure  (Pg) 

Mach  number  sensor  (AP/Pg) 

Turbine  discharge  pressure  (Pg) 

Burner  temperature  (T^) 

Blank 
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Channel  Description 

12  Turbine  discharge  temperature  (T&) 

13  Fuel  nozzle  pressure  differential  (Pfn  - Pg) 

14  Compressor  pressure  at  bleed  opening  (PgLD) 

Reference  signals  for  tape  playback  were  recorded  on  the  first  three  chan- 
nels. Sensor  signals  were  recorded  on  the  other  11  channels,  except  for 
Channel  11  which  was  left  blank. 

All  of  the  record  and  reproduce  amplifiers  except  for  Channels  1,  2,  and  3 
were  carefully  calibrated  for  ±2.  0 volts  full  scale  before  recording  the  data. 
The  gain  factors  on  Channels  1,  2,  and  3 were  adjusted  according  to  the 
requirements  of  the  BAFCO  servoanalyzer. 

In  addition,  the  frequency  response  characteristics  of  the  analog  computer 
amplifiers  and  the  tape  recorder  amplifiers  were  also  measured  before  the 
engine  tests  were  begun.  It  was  determined  that  these  components  would 
not  affect  the  engine  frequency  response  data. 
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ENGINE 


Ml! 


DC  REJECT  AND 
AMPLIFICATION 


Figure  A-l.  Engine  Test  Configuration 
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B.  FUEL  SYSTEM  TRANSFER  FUNCTION 

Figure  A -2.  Fuel  System  Instrumentation 
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A.  ON-LINE  MODE 


B.  TAPE  MODE 


Figure  A -3.  Data  Analysis 
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Figure  A -4.  BAFCO  Servoanalyzer  Block  Diagram 


Figure  A -6.  Sensor  Signal  Recording  Process 


APPENDIX  B 

MODEL  IDENTIFICATION  PROCEDURE 


INTRODUCTION 

The  computer  program  TFNS  used  to  identify  transfer  function  models  from 
experimental  frequency  response  data  is  discussed  in  this  appendix.  The 
principal  topics  discussed  include  a theoretical  treatment  of  the  identification 
procedure,  a description  of  the  algorithm  flow  chart,  and  a description  of  the 
input  data  cards  used  to  run  the  program. 

The  algorithm  presented  is  based  on  the  identification  method  proposed  by 
Yutaka  Suzuki  in  Reference  7.  Suzuki' s method  provides  a systematic  means 
of  finding  the  transfer  function  model  of  lowest  order  which  best  approximates 
the  experimental  frequency  response  data.  This  method  does  not  depend  on 
knowledge  of  the  order  of  the  system;  the  transfer  function  of  lowest  possible 
order  is  automatically  identified.  Details  of  the  method  are  discussed  in  the 
following  section. 


THEORY 

Before  beginning  the  actual  theoretical  discussion,  the  notation  used  in  this 
section  is  introduced.  The  symbol  is  used  to  represent  the  experimental 

frequency  response  data.  This  information  can  be  expressed  in  either  of  two 
ways;  as  real  and  imaginary  components  of  the  frequency  response  locus,  or, 
alternatively,  as  amplitude  ratio  and  phase  shift  data.  The  equivalence  of 
these  two  representations  is  illustrated  by  Equation  (B-l): 

Gfju^)  = RL  + jl.  = A.  cos  6.  + jA.  sin  0.  (i  = 0,  1,  . . . r)  (B-l) 
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where 


R = Real  part  of  frequency  response 
I = Imaginary  part  of  frequency  response 
A = Amplitude  ratio 
6 = Phase  shift 

For  convenience,  the  first  notation,  real  and  imaginary  parts,  is  used 
throughout  the  appendix.  The  subscripts,  i,  on  the  variables  indicate  that 
the  frequency  data  are  discrete;  i.  e. , the  frequency  data  are  composed  of 
ordered  pairs  (R^,  I.)  measured  at  (R  + 1)  different  frequencies  u.  (i  = 0, 

1,  ...  r). 

The  transfer  function  model  to  be  identified  from  the  <perimental  data, 
is  represented  by  the  following  rational  function; 


where 

m = Order  of  numerator  polynomial,  N 
n = Order  of  denominator  polynomial,  D 

& 

The  subscript,  a,  in  this  notation  stands  for  "approximation.  " In  general, 
the  orders  of  the  numerator,  N , and  the  denominator,  D , will  not  be  equal, 
but  all  physical  systems  are  characterized  by  the  inequality,  n * m. 

The  complex  variable,  s,  has  been  used  in  Equation  (B-2)  to  avoid  confusion. 
Rewritten  in  terms  of  the  variable,  ju,  the  equation  is; 


Ga(s>  = DlSi 
a 


Na(s)  _ b0  * V * V2  + V™ 


1 + a^s  + &2b  + • • • 3 j®11 


(B-2) 
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+U 


2 2 
+ ...)  + ju(b.,  - b«u  +...) 

G/j,.\  _ «■ ~ _ u & l o 

c ^a'*1  (1  - agio  + . . . )+ju(a1  - agu  + . . . ) 


Nft(ju) 


(B-3) 


The  objective  of  the  identification  procedure  is  to  determine  a transfer  func- 
tion, G (ju),  which  accurately  approximates  the  experimental  data,  G(jw). 
This  objective  entails  solving  two  separate  problems.  First,  it  is  necessary 
to  determine  what  order  [values  of  the  constants  n and  m in  Equation  (B-2)] 
the  transfer  function  model  must  be  to  obtain  a good  approximation,  and, 
secondly,  once  the  order  has  been  established,  the  unknown  coefficients  a^, 
a2'  " * an'  b0'  bl'  “ * bm  must  1)6  calculated  by  some  technique. 

The  identification  method  suggested  by  Suzuki  contains  a systematic  approach 
for  solving  both  of  these  problems  simultaneously.  The  procedure  is  com- 
posed of  two  iteration  loops:  the  outer  loop  addresses  the  problem  of  finding 
the  order  of  the  model  (n  and  m),  whereas  the  coefficients  a1#  a2,  ...  an, 
bQ,  bl'  * * * bm  are  calculated  in  the  inner  loop. 

To  start  the  procedure,  the  user  makes  an  initial  guess  for  the  variables 
n and  m,  called  nQ  and  mQ.  Then  the  coefficients  aj,  a2,  ...  a-,  bQ,  b^, 

. . . bmQ  corresponding  to  the  transfer  function  model  of  order  nQ,  mQ  are 
calculated  in  the  inner  iteration  loop.  Once  these  parameters  have  been 
found,  the  frequency  response  of  the  model,  G (ju),  is  compared  with  the 

cl 

experimental  data,  Gfju^.  If  the  difference  is  small,  the  calculation  is 
stopped.  However,  if  the  difference  is  large  (indicating  that  the  model  is  not 
a very  good  approximation  to  the  experimental  data),  the  procedure  is  not 
stopped.  Instead,  the  order  of  the  transfer  function  model  is  increased  by  1 
(nj  = nQ  + 1,  mj  = m^  + 1)  and  the  search  for  the  coefficients  a^,  a2,  ...  anQ, 
anj,  bQ,  bj,  . . . bmg,  bmj  is  begun  ane  v.  Note  that  this  time  there  are  two 
additional  coefficients,  a^  and  bm^,  which  must  be  found  because  the  order 
of  the  transfer  function  model  has  been  increased  by  1.  The  procedure  is 
repeated  until  a transfer  function  model  of  suitable  accuracy  has  been 
obtained. 
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Since  the  order  of  the  transfer  function  model  is  increased  systematically  in 
the  outer  iteration  loop,  the  procedure  results  in  the  transfer  function  of 
lowest  possible  order  being  found.  This  is  a highly  desirable  outcome, 
signifying  that  minimum  complexity  has  been  achieved. 

The  important  calculations  included  in  the  two  iteration  loops  are  considered 
in  the  following  paragraphs. 

In  the  inner  iteration  loop  the  unknown  parameters  a^,  ag.  . . . an,  bQ,  b^, 

. . . b in  the  approximation  G (ju)  are  calculated  to  minimize  the  following 

ill  a 

error: 


r 

E = £ |G(ju.)-  Ga(jUi)|2  (B-4) 

i=0 

This  error  is  called  the  real  error. 

Several  numerical  techniques  could  be  used  to  find  the  set  of  coefficients 
which  minimize  this  error.  However,  most  of  the  techniques  require  many 
iterations  because  of  slow  convergence  in  the  vicinity  of  the  extremum.  To 
avoid  this  difficulty,  the  following  method  is  used. 

First,  weighting  functions,  0.,  are  added  to  the  error  equation  and  the  terms 
rearranged  to  give: 

r 

En=  Y.  *1  liyiM,)  GflUj)  - Na(jUi)|2  (B-5) 

i=0 

The  set  of  parameters  a^,  a2,  ...  an,  bg,  b^,  . . . bm  which  minimizes 
Equation  (B-5)  is  determined  by  solving  the  system  of  first-order  simul- 
taneous equations  derived  by  differentiating  E with  respect  to  each  of  the 
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unknown  parameters  and  setting  the  result  equal  to  zero.  That  is,  the 
following  system  of  equations  is  solved  for  a^,  a2,  ...  an,  b0,  b^,  ...  bm: 


0 


= 0 


= 0 


In  matrix  notation  these  equations  can  be  expressed  as: 


Ap  = c 


(B-6) 
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b0 

'"o' 

bl 

Ti 

b2 

S2 

b3 

T3 

• 

0 

and  c = 

0 

• 

0 

a2 

u2 

a3 

0 

0 

• 

0 

• 

• 

• 

Jm 

xh=E  #i“ih  iui 


i=0 

r 


Sh  = Z 0iwih  RiAui 

i=0 

r 

t.  = y 0.w. h i Aw. 

i=0 

uh  = f «rth(Ri2  + Ii2)i“l 

i=0 

(h  = 0,  1,  2 ...) 


These  equations  are  derived  in  Reference  7.  The  unknown  coefficients  a^, 
a2,  . . . an.  b0,  bj.  . . . bm  are  contained  in  the  vector,  p. 

The  weighting  coefficients,  0.,  are  chosen  to  be: 

0,  - -y  (1=1.  2,...  r)  (B-7) 

|Da(H>!2 


so  that  when  they  are  substituted  into  the  error  equation  [Equation  (B-5)"l, 
it  reduces  to  E [defined  in  Equation  (B-4)]. 


However,  since  D (ju)  is  unknown  until  the  system  of  equations,  Ap  = c,  is 

a 

solved,  an  iterative  technique  must  be  used  to  obtain  0.  which  approximates 
1/|D  (jw.)l  as  closely  as  possible.  The  following  method  is  used: 

cl  I* 


1 

(k-1)  12 

W I 


(B-8) 
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Namely,  at  the  kth  iteration,  0^)  [s  determined  from  the  results  of  the 
(k-l)**1  iteration.  By  approximating  0.  in  this  manner,  it  gradually  approaches 


This  iteration  on  0.  forms  the  basis  of  the  inner- loop  iteration.  Parameters 

(k-l)  (k-l),  ' (k-l)  (k-l)  (k-l)  _ obtained  from  the 

1 2 n 0 1 m (k) 

(k-l)in  iteration  are  used  to  estimate  the  weighting  coefficients,  0.'  for  the 
k**1  iteration.  These  coefficients  are  substituted  into  the  system  of  equations, 
Ap  = c [Equation  (B-6),  which  is  then  solved  for  new  estimates  of  the  param- 


eters a 


(k)  „ (k) 


(k)  . (k) 
» D0 


‘1  , ~2  , . . . «n  , w0  , bx  . - - - -m 

tinues  until  the  weighting  coefficients,  0. , have  converged  to  a limit. 


<k) 


(k) 


This  procedure  con- 


Whether  or  not  0^  converges  is  a key  point  In  this  algorithm.  Convergence  of 
0.  implies  that  a minimum  of  the  error  function,  E,  in  Equation  (B-5)  has  been 
found.  If  0i  does  not  converge,  the  identification  procedure  will  not  work. 
Although  a sufficiency  condition  for  convergence  of  0.  has  not  been  discovered, 
in  practice,  convergence  has  been  obtained  in  every  case. 


The  normalized  error  criterion,  E 
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E. 


= — y 

r + 1 L 


i=0 


1 - 


D.(k)(jw.) 


a 


D_  (ju.) 


a 


(B-9) 


referred  to  as  the  equation  error,  is  introduced  to  monitor  the  convergence 
of  0^.  When  the  equation  error,  E^,  becomes  less  than  a prescribed  error 
tolerance,  c^,  the  inner-loop  iteration  is  stopped,  since  this  implies  that  0^ 
has  converged. 


After  convergence  of  0^  has  been  obtained,  the  real  error,  EN,  defined  in 
Equation  (B-5)  is  computed  in  the  outer  iteration  loop  to  ascertain  the  accu- 
racy of  the  approximation,  Gft(ju).  The  value  of  EN  is  compared  with  a pre- 
scribed error  tolerance,  cN,  to  determine  what  action  is  to  be  taken. 
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If  En  is  less  than  or  equal  to  cN,  the  algorithm  is  halted.  This  outcome 
indicates  that  a transfer  function  model  of  sufficient  accuracy  has  been 
identified. 

The  other  alternative,  EN  greater  than  c^,  indicates  that  further  computation 
is  necessary.  This  outcome  implies  that  the  assumed  order  of  the  approxi- 
mation is  not  large  enough  to  accurately  model  the  experimental  frequency 
data.  Therefore,  in  this  case  the  orders  of  the  numerator,  N (ju),  and  the 
denominator,  D (ju),  of  G (ju)  are  increased  by  1 and  control  is  transferred 
back  to  the  inner  iteration  loop.  The  procedure  is  repeated  until  a transfer 
function  model  with  EN  £ cN  is  identified. 


PROCEDURE 

A flow  chart  of  the  computer  algorithm  is  presented  in  Figure  B-l  and 
discussed  in  the  following  paragraphs. 

First,  the  constants  n and  m,  and  the  orders  of  D (ju)  and  N (ju)  at  the 
beginning  of  the  computation,  are  read  in.  Although  these  constants  can  be 
set  to  1 in  any  case,  it  is  often  possible  to  roughly  estimate  the  order  of 
G (ju)  from  the  experimental  frequency  response  data.  Needless  computations 
are  avoided  by  selecting  reasonable  values  of  n and  m at  the  beginning.  If  the 
selected  order  is  not  high  enough  to  approximate  the  experimental  data,  the 
computer  algorithm  will  automatically  increase  the  order.  This  is  often  the 
case. 

Next,  the  equation  error  and  real  error  tolerance  parameters,  c and  c^, 
are  read  in.  The  constant  c^  determines  how  hard  the  computer  must  work 
to  find  a transfer  function  model  of  a given  order.  This  constant  controls  the 
inner  iteration  loop  in  the  program.  When  the  weighting  coefficients,  0. , 
have  converged  to  within  the  tolerance  specified  by  c^,  the  inner- loop  itera- 
tion is  stopped  and  the  accuracy  of  the  transfer  function  model  is  checked. 
Experience  indicates  that  a value  of  about  0.  01  for  c^  will  give  good  results. 
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The  real  error  tolerance,  c^,  determines  how  accurate  the  transfer  function 
approximation,  Ga(ju),  must  be.  This  constant  controls  the  outer  iteration 
loop  in  the  program.  After  the  inner  loop  has  converged,  the  integral  square 
error,  EN,  defined  in  Equation  (B-5)  is  checked.  If  the  error,  EN,  is  found 
to  be  less  than  c^,  the  computation  is  stopped,  since  this  implies  that  a good 
approximation  has  been  obtained.  However,  if  the  error,  EN,  is  larger  than 
cN,  a good  approximation  has  not  been  obtained.  In  this  case,  the  orders  of 
N (ju)  and  D (ju)  are  increased  by  1 and  the  computations  are  begun  anew. 

It  has  been  discovered  that  a value  of  about  0.  01  for  cN  will  produce  satis- 
factory results. 

After  the  error  tolerance  constants  c^  and  cN  have  been  read  into  the  program, 
the  experimental  frequency  response  data,  G(ju.),  is  read  in.  Provision  has 
been  made  in  the  program  to  input  this  information  as  Bode  gain  and  phase 
data.  Nyquist  polar  data,  or  Nyquist  rectangular  data.  The  data  is  converted 
to  rectangular  coordinates  for  internal  use  in  the  program. 

Next,  the  weighting  coefficients  0i  are  initialized  to  a value  of  1 and  the  inner 
iteration  loop  counter  k is  set  to  1. 

Then,  the  coefficients  in  the  matrix  A and  vector  c of  the  equations  Ap  = c are 
calculated  and  the  equations  solved  for  the  vector  of  unknowns,  p. 

Following  this,  the  real  error,  E^,  given  by  the  following  expression: 

r 

En  * Y,  h I'W  G(j“i)  • Na°“i)l2 
i=0 

is  calculated  and  the  result  stored  for  future  comparison  with  the  error 
tolerance  parameter,  cN. 
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(U) 

Then,  the  denominator  polynomial,  D (ju.),  is  evaluated  at  the  u.  (i  = 0, 

1,  ...  r)  corresponding  to  the  experimental  frequency  response  data.  These 

numbers  D (ju.)  will  be  used  later  on  to  update  the  estimates  of  0.  if  an  addi- 
a i i 

tional  inner  loop  iteration  is  necessary. 


At  this  time  the  inner  iteration  loop  counter  k is  checked  to  see  if  its  value  is 
1.  A value  of  1 indicates  that  this  is  the  first  time  through  the  inner  iteration 
loop.  If  k equals  1,  the  next  calculation  cannot  be  performed,  so  control  is 
transferred  to  the  end  of  the  inner  iteration  loop  where  the  weighting  coeffi- 
cients 0.  are  updated. 


If  k is  not  equal  to  1,  the  equation  error  E . is  computed  as 

‘ 0 


na(k>  (fot> 

Da(k‘1)  <JU;> 


and  compared  with  the  equation  error  tolerance  parameter,  c^  If  is 

greater  than  the  allowable  tolerance,  c .,  the  inner  loop  is  n^+  converged. 

9 

In  this  case  the  weighting  coefficients  0.  are  updated  according  to  the  relation 


1 

|Da<k)  <J“i>|2 


and  control  is  transferred  back  to  the  beginning  of  the  inner  iteration  loop 
denoted  Station  2 in  the  flow  chart. 

However,  if  E . is  less  than  the  tolerance,  c .,  the  inner  loop  iteration  on  0. 

9 9 1 

has  converged  and  control  is  transferred  to  a check  on  the  real  error,  EN> 

E^,  which  has  been  computed  in  a previous  step,  is  compared  with  the  real 

error  tolerance  parameter,  cN#  If  EN  is  less  than  c^,  the  computation  is 

stopped,  since  this  implies  that  the  transfer  function  model  obtained  is  a good 
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approximation  of  the  experimental  frequency  response  data.  However,  if  EN 
is  greater  than  c^,  the  approximation  is  not  good  enough.  In  this  case,  the 
orders  of  the  numerator,  N (ju),  and  denominator,  D (ju),  in  the  approxi- 
mation  G (ju)  are  increased  by  1 and  control  is  transferred  to  the  beginning 

a. 

of  the  outer  iteration  loop.  Station  1 in  the  flow  chart.  The  iteration  proce- 
dure continues  until  an  accurate  transfer  function  approximation  has  been 
obtained. 


TFNS  PROGRAM 

The  identification  algorithm  discussed  in  the  previous  sections  has  been 
assembled  into  a Fortran  program  called  TFNS.  In  addition  to  performing 
the  calculations  associated  with  the  identification  algorithm,  the  program  also 
contains  a section  which  computes  the  frequency  response  of  the  transfer 
function  approximation.  These  data  are  crossplotted  with  the  experimental 
frequency  response  data  on  a single  Bode  plot  to  permit  rapid  visual  verifi- 
cation of  the  accuracy  of  the  transfer  function  model  obtained.  A listing  of 
the  Transfer  Function  Identification  Program  is  presented  in  Tables  B-l 
through  B-12.  The  input  data  cards  needed  to  run  the  program  are  identified 
below. 

The  input  data  deck  consists  of  two  groups  of  cards,  a program  control  group 
and  a frequency  response  group.  Each  group  ends  with  a card  with  an 
asterisk  (*)  in  Column  1.  The  cards  in  each  group  are  identified  below. 


Program  Control  Group 

These  cards  have  a flag  character  in  Column  1 and  parameters  in  10- character 
fields  thereafter.  These  cards  may  be  in  any  order,  and  there  may  be  more 
than  one  of  each  type  in  the  group.  However,  only  data  from  the  last  card  of 
a particular  type  in  the  group  are  considered  in  that  run. 
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Card  With  "C1*  in  Column  1 — This  card  inputs  parameters  pertaining  to  the 
linear  least  squares  algorithm.  If  a fit  is  to  be  performed  in  a run,  A "C" 
card  must  appear  in  the  Program  Control  Group: 

1)  Denominator  Starting  Order,  Columns  2-10  --  Right- justified 
integer. 

2)  Numerator  Starting  Order,  Columns  11-20  --  Right-justified 
integer.  Should  be  s denominator  starting  order.  The  dif- 
ference between  numerator  and  denominator  order  is  preserved 
as  the  orders  are  increased  by  the  program  for  better  fit. 

3)  Equation  Error  Tolerance,  Columns  21-30  --  Free-field 
floating-point.  This  determines  how  hard  the  computer  should 
work  for  a fit  at  a particular  order.  Recommended  value  = 0.  01, 
but  can  be  tightened  up  to  achieve  lower-order  fit. 

4)  Real  Error  Tolerance,  Columns  31-40  --  Free-field  floating- 
point. When  the  computer  has  ground  out  a transfer  function 
at  a particular  order  which  satisfies  the  equation  error  toler- 
ance, it  checks  a discrete  approximation  of  the  integral  squared 
error  over  the  entire  curve  against  the  real  error  tolerance. 

If  the  real  error  tolerance  is  satisfied,  the  algorithm  halts. 

If  not,  numerator  and  denominator  orders  are  increased  by  1. 

5)  Frequency  Scale  Factor,  Columns  41-50  --  Free-field  floating- 
point.  This  is  used  to  center  the  range  of  the  logarithms  of 
input  frequencies  around  zero.  If  this  scale  factor  is  not  used 
when  input  frequencies  are  greater  than  10  radians  or  less  than 
0. 01  radian,  then  ill- conditioning  in  the  matrix  used  in  the  linear 
least  squares  algorithm  will  result.  For  example,  if  the  lowest 

• frequency  is  1 radian  and  the  highest  is  100  radians,  then  a 
suitable  scale  factor  would  be  0. 1. 


129 


6)  Term  Rejection  Criterion,  Columns  51-60  --  Free-field  floating- 
point. Sometimes  the  algorithm  is  forced  to  compute  a higher- 
order  numerator  than  is  needed  for  a good  fit.  When  this 
happens,  it  comp  tes  a polynomial  coefficient  which  is  very 
small.  Inclusion  of  such  a term  causes  ill- conditioning  which 
affects  the  accuracy  of  the  computation  of  the  zeros.  To  pre- 
vent this,  there  is  logic  to  zero- out  high- order  coefficients 
which  are  too  small  to  affect  the  transfer  function.  If  the 
numerator  polynomial  is  represented  as 

a0  + al  * + a2  + “““  + am  * 

and  if  there  are  r frequency  data  points,  (i  = 0,  1,  . . . r), 

then  the  term  a is  tested  as  follows.  If 
m 

Y.  Iam(j“i,m 

i=l 

r m . 

i=l  k=0 

is  less  than  term  rejection  criterion,  then  reject  am,  and 
reduce  numerator  order  by  1. 

Card  With  "O"  (the  letter)  in  Column  1 — This  card  controls  the  output  Bode 
plot.  It  must  be  included  if  a Bode  plot  is  desired  for  a run.  The  Bode  plot 
always  gives  four  decades: 

1)  Number  of  Points /Decade,  Columns  2-10  — Right-justified 
integer. 
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2)  Number  of  Points  Evaluated  For  Plot  Routine,  Columns  11-20  — 
Right- justified  integer  must  be  * number  of  points /decade. 

3)  Plot  Starting  Frequency,  Columns  21-30  — Pree-field  floating- 
point. 

4)  Gain  Scale,  Columns  31-40  --  Free-field  floating-point: 

- If  gain  scale  > 0,  plot  scale  = ± (gain  scale)  db. 

- If  gain  scale  > 0 or  blank,  plot  scale  chosen  auto- 
matically from  plot  values. 

Card  With  "Z"  or  MP"  in  Column  1 --  It  may  be  desirable  to  neglect 
extraneous- looking  poles  or  zeros  when  evaluating  the  computed  transfer 
function  for  a plot.  This  is  done  with  a root  deletion  card.  The  roots  to  be 
deleted  are  punched  left-justified  in  A4  fields  from  left  to  right  on  the  card. 
Each  field  begins  with  a Z or  P,  meaning  zero  or  pole,  followed  immediately 
by  a two-digit  number  corresponding  to  the  number  of  the  desired  pole  or 
zero  indicated  on  the  output  of  the  previous  run.  For  instance,  if  pole  3 and 
zero  2 from  the  previous  run  are  to  be  deleted,  then  the  card  would  appear 
as  follows: 


P03  Z02 

t t 

Columns  1 5 

The  poles  and  zeros  may  be  listed  on  the  card  in  any  order  whatsoever,  but 
the  first  blank  A4  field  terminates  the  scan  of  the  card.  Thus,  up  to  20  poles 
and  zeros  may  be  deleted  in  one  run.  Note  also  that  if  a number  refers  to  a 
conjugate  pair,  both  members  of  the  pair  are  deleted,  and  the  order  of  either 
the  numerator  or  the  denominator  is  reduced  by  2 by  such  a deletion.  The 
output  of  a deletion  run  will  print  "DEL”  in  the  appropriate  entry  of  the  polt- 
zero  listing.  Also,  deletions  can  be  made  on  a fit  run,  if  pole  and  zero 
numbers  are  known  a priori. 
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End  program  control  group  with  card. 


Frequency  Response  Cards 


These  cards  have  input  data  points  on  them,  one  per  card.  The  cards  must 
be  in  order  of  ascending  frequency.  The  frequency  can  be  in  He  or  radians 
per  second,  and  the  units  need  not  be  consistent  from  one  card  to  the  next. 
However,  they  must  still  be  in  order,  as  if  they  were  in  consistent  units. 

In  other  words,  10  radians  must  appear  before  5 Hz. 

These  cards  have  a common  format: 

• Column  1 - Data  Type  Flag: 

- Blank  = Bode 

- B = Bode 

- P = Nyquist  polar 

- R = Nyquist  rectangular 

• Columns  15-29  - Frequency,  Free- Field  Floating-Point: 

- If  H in  Column  15,  frequency  is  in  Hz. 

- If  R in  Column  15,  frequency  is  in  radians  per  second. 

• Default:  If  blank  or  frequency  entry  starts  in  Column  15 , 
frequency  is  in  radians  per  second. 

• Columns  30-44  - Gain,  Magnitude,  or  Real  Part;  Free-field 
floating-point.  Defaults  vary  according  to  Column  1 entry, 
but  D in  Column  30  indicates  decibels,  and  M in  column  30 
indicates  magnitude. 
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• Columns  45-59  - Phase,  Angle,  or  Imaginary  Part:  Defaults 
vary  according  to  Column  1 entry,  but  D in  Column  45  indicates 
degrees,  and  R in  Column  45  indicates  radians. 


Defaults  Chart 


Col.  1 

Frequency  Field 

Part  1 

Part  2 

Blank 

rad/sec 

Gain  (db) 

Phase  (deg) 

B 

rad/sec 

Gain  (db) 

Phase  (deg) 

P 

rad/sec 

Magnitude 

Angle  (rad) 

R 

rad/sec 

Real 

Imaginary 

As  indicated,  defaults  on  units  in  frequency,  gain,  phase,  magnitude,  or 
angle  can  be  overridden  by  punching  the  appropriate  character  (H,  R,  D,  M) 
in  the  first  column  of  the  appropriate  field. 

End  frequency  response  data  group  with  card. 
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Figure  B-l.  TFNS  Flow  Chart  (Concluded) 
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Table  B-l.  Transfer  Function  Identification  Program  --  Program 
to  Identify  S- Domain  Transfer  Function 

tt 


1 0000P000 

t Oooooool 
1 0000000* 

4 OOOOOOOS 

5 0000000^ 

4 OOOOOOOS 
7 0000000* 

• OOOQOOO^ 
9 OOOOOOIO 

10  OoOoOoil 

11  OOOOOOIZ 

12  00000013 

13  00000019 

14  OoOQOOlS 

15  OOOOOOl* 

14  OOOOOOI* 
17  00000020 

15  00000021 

19  00000022 

20  00000023 

21  00000029 

22  00000025 

23  0000002* 

24  0000002? 

25  00000030 

26  00000031 

27  00000032 

28  00000033 

29  0000003* 

30  00000035 

31  0000003* 

32  00000037 

33  00000040 

34  00000041 

35  00000042 

36  00000043 

37  00000044 

38  00000045 

39  00000046 

40  00000047 

41  00000050 

42  00000051 

43  00000052 

44  00000053 

45  00000054 

46  00000055 

47  0000005* 

48  00000057 

49  00000060 

50  00000061 

51  00000062 

52  00000063 

53  00000064 

54  0000006® 

55  0000006* 

56  00000067 

57  00000070 

58  00000071 

59  00000072 

60  *0ooo0o73 


T8  jOtNTlrv  S-DOMAIN  TRANSFER  FUNCTION 

riTjjNa  prcouency  response  data  points  . . . 

COMMON  N ( 50 ) < DELW  ( SO ) * RESP 1 50 ) # PM  I ( 50 ) # A ( 4 1 # *2 ) 

COMMON  S ( 43 ) , T ( 43 ) « U (43 ) # LAMBDA ( 43  ) 

COMMON/KUUOf/M.  kWUOC 

DIMENSION  P(42)#L0RD(4)  ... 

D I MENS  JON  FREQ IN  ( 1Q0 ) 4PART 1« 100  > , PART2 { 100 ) 

DIMENSION  LI*AaE<20)#jPT(20)#JZT<20) 

COMPLEX  3 A ( 5Q ) » TNUM#  Tp£N0M#  Cl # FR  AC# CPOL Y#  j W# TFA# CMLVF 
COMPLEX  OA<50>«OAOIDCSOL. 

complex  Zeros (22># poles (22) 

COMPLEX  "ESP 

real  lambda 
INTEoER  oType 
logical  C0MP,0uT 

DATA  NjTER/ioO/ 

DATA  NMAX/?o/#MMAx/SO/#NgEC/4/ 

DATA  ICR/5/#LP/6/#LBODE/lMB/#LPOLAR/iHP/#LRECT/lMlt/ 

DAT  A j L00° ( I i # I • 1 » 4 ) / 2HST , 2HN0#  2HRD#  2MTH/ 

DATA  OTYPE/InB/,  . . 

DATA  LC0Mp/ihC/*LSTAR/1H#/#LBLANK/1H  /#L0UT/1H*/ 

DATA  LP0LE/lMP/iLZER0/lH7/ 

call  restart 

nR.ATAN( 1*0)/4B» 

ROil./OR 
PI»4»0*ATaN( J»0) 

rerun  sense  SWITCH 

200  C*UL_SSWTj-wi  HON) 

I r ( IJ0n«Eq,2)ST0P 
COMP.. FALSE. 

OUT.jFaLSe. 

DO  91  J-1,20 
JpT { I ) .0 
JZT ( I ) .0 
91  CONTJNUE 

JPTP.O 
JZTP.O 

read  in  Parameters  and  roots  tb  be  deleted 

90  READi  IcR'iooi)LlMAQE 

DECOOE(1#1002#LIMAQE)LFLAG 

read  in  starting  denominator  brdeR/Startinq  numerator  order# 

ERROR  TOLERANCE. REAL  ERROR  TOLERANCE#FREqUENCY  SCALE 

FACTOR, TERM  rejection  criterion. 

IF (LFLaQ. Eo.LCOMP) DECODE (60# IOOO.LIMAGEiND.Nn.CS.CT/FSCAlE.TOLMAq 
COMPiCOMP.OR. (LPLAQ.EQ.LCOMP) 

OUT.OUT.Or.(lFLAG.EO.LOUT) 

READ  IN  No,  OF  POINTS/DECADE#NUMBER  OF  POInTS.PLOT  STARTING 

frequency, desired  gain  scale,  if  gain  scale  .le.o*  gain  scale 

CHOSEN  AUTOMATICALLY. 

IF (LFLAG»E3«l9UT) DECODE (40# 1000# L IMAGE  )NPDEC#NPOINTS.WLOt«»oGAINX 
IF(LFLAG*NE.LP0LE.AND.LFLAG.NE.LZER9)G0  TO  9g 
DO  93  1.1*20 

DECODE ( 3# 1003# L IMAGE ( I ) ) LABEL# IPT 
IF(LABEL*EQ.LBLAN<)G0  TO  92 
IF(LABEL*NE.lPOLE)GO  to  94  ' 

JPTP'.JpTPjj 
JPT ( jPTP ) . JPT 
GO  T»  93 

94  JZTPidzTPji 

JZT( jZTP) . IPT 
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61 

63 

63 

64 

65 

66 

67 

68 

69 

70 

71 

n 

73 

74 

75 

76 

77 

78 

79 

80 
81 
83 

83 

84 

85 

86 

87 

88 

89 

90 

91 
93 

93 

94 

95 

96 

97 

98 

100 

101 
103 

103 

104 

105 

106 

107 

108 

109 

110 
Hi 
113 
113 
116 
118 
116 
117 
HI 
119 
190 
111 


90 


8? 


83 

81 

80 


84 


84 
83 

85 


Table  B-l.  Transfer  Function  Identification  Program  — Program 
to  Identify  S- Domain  Transfer  Function  (Continued) 

00000076 

00000075 
00000076 
00000077 

O000O100 
O0O0O10I 
O0O3O1O2 
OOOOOIO2 

00000106 

O0OOOIO5 
OgOgOlO6 

8orn01n7 
OO0O11O 
OoOoOlll 
00000113 
00000113 
OoOoOl 14 
00000115 
00000116 
00000117 
0000°12° 

OoOoOl 2 1 
00000122 
00000123 
00000124 
00000125 
00000126 
00000127 
00000130 
00000131 
00000132 
00000133 
00000134 

00000135 
00000136 
00000137 
00000140 
0000°l4l 
00000142 
00000143 
0000°146 
00000145 
00000146 
OoOoOl 47 
0000°150 
00000151 

ooooois* 

OOOOOIB3 
00000156 
O0OOOI88 
OoOoOlS* 

OOOoOlS-7 
OoOoOl 60 

Oooooui 
0000016* 

00000163 
OoOoOl 64 
OoOoOl 68 

000001** 

OOOQOJ** 

OOOOOITO 


93  CONTINUE 

92  IF(LFLA0*NE»lSTAR)GO  TO 

C***  ..SORT  JPT  AND  JZT  ARRAYS 
lF(JPTp*LE.i'AQO  TO  80 
UPTPhI.jPTP-1 
DO  81  1*1# JPTPM1 
IP1- 

00  82  J»lpi # jPTP 
IF( JpT( J) •OCtJPT ( l ) )G0  TO 
ITEMp.jPTn) 

JPT(|).JPT(J) 

JPT( jlilT^HP 
CONTINUE 
CONTJNUE 

!F(JzTptUE.i)QO  TO  85 
JZTPmI-jZTP-1 
DO  83  ’.I4JZTPHI 
lPl«IM 

DO  8J  j»lPl«uZTP 
IF(JzT(J).3E.JZT{ll)00  to 
.ITEMP.j2T{  J) 

JZT(p.jZT(J» 

JZT(j)ilTEMP 
CONTINUE 
CONTINUE 
CONTINUE 

IFi.n0T«C0MPiQ0  TO  100  . 

WRITEaP43060)ND4NN,CS.CTfFSCAtE,WLBw,NPOEC»NP8lNT8iTOUMAQ 

read  In  frequency  response  data  and  convert  from  External 

type  to  Internal  typE(REctangularj 

WR!TE(lP*20Ql) 

CAIL_CAR0r0(  !CR,LP|M|MMAx> 

SUBTRACT  time  delay  from  raw  data 

if (Sense  §wjtcm  2>  55#56 
wRITE(84  3ggo') 

R£AD( 8*30011  TDELAY 
WR1Te<LP*3003)  TDELAY 
DO  56  J • 1 # M 

TEMPO  f CaBS(«eSP( in  . ... 

TEMPp  ■ ATAfI3(AlMA0<PESP(in»REAL(RESP(I)n»RD 
TiMPp«VEMpp^(  I)«tDELAYt57»293 
TEMPp»TfMpP#0R  ... 

RESP ( I » • T£MPfl«CMPLX  < COS ( TEMPP > » S I N(  TEMPP n 
CONTINUE 
C*'*«’SCAL£  DOwN  FREQUENCY 
DO  6p  J»1*M 
W(  d,w<I).fscalE 

c*«* * • initialize  phIiCompute  delw#and  find  din'a’minatrR  real  err*r» 

ERPlRDyO* 

PHI<1  J»l» 

DELW(l)«<W(l)4W(2>)/t* 

TEMP.C^8fiRE8Pll)) 
er^0rD|EProro4Tcmp«timp«delw  m 
mmi«m»j 

DO  S, I»?AMM1  . . 

OCUWt 

PHKJltl*  .. 

T|MP»Cf»J( Rfl?1 1 ) 1 
ERROR!'  ,ERrOR06TE«P4 


55 


56 


l•TCMP•OC^H(l) 
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Table  B-l. 


122  0o0o0l7l 

123  00000172 

12*  0o0o0i?3 

125  0000017* 

126  00000175 

127  00000176 

128  000001?7 

129  00000200 

130  00000201 

131  0000020? 

132  00000203 

133  0000020* 
13*  O00CO2O5 

135  00000206 

136  0o0o0207 

137  00000210 

138  OjOoOail 

139  OoOOC2l2 

1*0  00000213 

1*1  0o0o0?l* 

1*2  00000215 

1*3  00000216 

t*f  00000217 

1*5  00000220 

1*6  O0O9O22I 
1*7  0000022 2 
1*8  00000223 

1*9  0000022* 

150  00000223 

151  0000022* 

152  00000227 

153  00000230 
15*  00000231 

155  00000232 

156  00000233 

157  0000023* 

158  00000235 

159  00000236 

160  00000237 

161  000002*0 
162  000002*1 
163  000002*2 
16*  000002*3 

165  000002** 

166  000002*5 

167  000002*6 

168  000002*7 

169  00000250 

170  00000251 

171  00000252 

172  00000253 

173  0000025* 
17*.  00000255 

175  00000256 

176  00000257 

177  00000260 

171  OOOQ0261 
179  OoOQO«62 
ISO  0000026* 
111  0000026* 
112  00000265 


Transfer  Function  Identification  Program  — Program 
to  Identify  S- Domain  Transfer  Function  (Continued) 

5 CONTINUE 
PHI ( M) al • 

0ELWjM)»W(M).W(H-l) 

TEHP.C*BS(RESP(*in 
ERRORDiERR«RD*T£MP#TEMP»DELW ( M > 

ID*MAX»2#NMAx*i 

ci-<T»*o») 

*.«  I TE  < LP*2013 ) 

begin  outer  loopi incrementing  n> 

6 «»0 
NR.Nn*n0*{ 
nc.nr+1 
NDPlaND+1 
Npi.NN^l 
Np2«NN»2 

BEGIN  inner  l08P(  INCREMENTING  k> 

7 

IF<K«GT*NjTER)ST0P  66666 

CALCULATE  a- MATRIX  and  C»vECT0R 

CALL  LDACjM,ND,NN)  . 

SOLVE  A*P«C  FOR  P 

CALL  TOlNyRt ISOL<IDSOL#NR,.NC#A, IDAMAX*kwA#DeT) 

IF( < ISOL^IOSOL) »GT»2)ST0P  77777 
CALL  M0VE(A(iiNC)/P»NR#NP2) 

CALCULATE  VALUE  OF  COMPUTED  TRANSFER  FUNCTION  AT  FREQUENCY 

DATA  POINTS  AND  THEN  EVALUATE  ERROR  FUNCTION. 

EN'.O. 

DO  8 1 a 1 * M 
JWaCMPLX<0*»W<!))  _ 

TNUM,CP8LY(P(1MNPI*  JW) 

TDENBMiCPOLY ( P ( NP2 ) # NDP I « JW ) 

DA( I j iTOEnOM 
GA<  liaTNUM/TOENOM 

temp,cabS(ResP(I)-gA{ n > 

TEMP.TEMP.TEmP.OELWH  ) 

ENaENaTEMp 

8 CONTINUE 

en;en/er*ord 

IF(K.Eq«1)G0  TO  10 

calculate  approximate  equation  error 

ECaO. 

DO  9 i;i«M 
FRACaOAt I )/DAOlD(  I J 
ECaEc*CABS(Cl-FRAC) 

9 CONTINUE 
ECiEC/Ft-OATiM) 

WRITE{LF»2002)EN#K£NO#NN#EC 

!F,Ec.LE*C8)08  TO  H 

C» • • • .MOVE  DA  To  DaOlD  AND  CALCULATE  PM! 

10  PO  12  !>1«M 
DAOLpj I)«DA( !j  . 

TEMP,CAB*«OAOl.Dlln 

PHI(J)*1vTEmF/TEMP 
12  CONTINUE 

00  TO  7 

£.*• .tCHECx  FOR  REAL  ERROR  CONVERGENCE  AND  IF  MAX*  ORDER  REACHED* 

11  IF<EN»LE»CT|flO  TO  16. 

IF J NO»|Q*NMAx)dO  TO  13 
NNeNN*l 

ND#N0*1  ... 

IF ( NO  »fiT • ( M*I )) 8T0P58555 
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Table  B-l.  Transfer  Function  Identification  Program  --  Program 
to  Identify  S- Domain  Transfer  Function  (Continued) 


it:  8$8R* 

I >9  00090270 

lift  00000271 

187  00000272 

188  00000273 

189  0000027* 

190  0000027* 

191  00000278 

19*  00000*77 

193  00000300 

19#  00000301 

195  0000030? 

196  .00000303 

197  O0ODO3O* 

198  00000305 

199  00000306 

200  0000030^ 

201  00000310 

202  00000311 

203  00000312 

20#  00000313 

205  0000031* 

20#  00000315 

207  0000031* 

208  00000317 

209  00000320 

210  O0O0O32I 

211  00000322 

212  00000323 

213  0000032^ 

21#  0000032* 

215  0000032* 

216  00000327 

217  00000330 

218  00000331 

219  0000033? 

220  00000333 

221  0000033* 

222  00030335 

223  0000033* 

22#  00030337 

225  00000340 

226  0o0o03#l 

227  000003#2 

228  000003#3 

229  000303#* 

230  0000°3#5 

231  O0O0O36* 

232  000003#7 

233  00000350 

23#  00000351 

235  0000°352 

236  00000353 

237  0003035* 

238  00000355 

239  0003035* 

2**0  00000357 

2#i  00000360 

2#2  00000361 

2#3  00000362 


is 

i* 

c«  •• •• 


#1 


WRlJf <*0*2003 ) NO# NN#FN 

CONTINUE 

REDUCE  8N9ER  Of  NUMERATOR 

WRJTEUf#*669> 

NEN0*NRi#N0fl 

WRITE (U^^gOOS* I > # » 

WRITE JLf *200*1  If (I)#l«Nf»#NEND) 
nnoRq«nn 


NOORO*NO 

IFfTpLMAG.EQ.OOQO  *5 
fMAGaO 1 
DO  #J  I»1#M 
JW»Cl/gMPLxifi»#W(I)) 

DO  #2  J«1#NPJ 

JW»JW#CNf|tX(0»|W(I)) 

PMAQ'PMAG«CAbS ( CMPUX I P( J I #0* ) • JW ) 

CONTINUE 

CONTJNUE 

PMAG.PMAG/FLeAT<M> 

DO  *3  !»ltNP) 

IRiNpl-I*l 

TMAGiO. 


DO  ##  J*1#M 

UW«CmPlX(Q.#WI j) >•»( IP>1) 
TMAG’.TMAG*CABStCMPLX(P(  !R»#0O#JW> 

##  CONTjNuE 

TmAG.TmAG/P|.8AT(M) 

RATI8»TMAg/PmAG 
IF(RaT10.gE.T0LMAG)G8  TB  #6 
Pi  I R j • 0 • 

NNPRd»NNOro.j 
#3  CONTINUE 

#6  CONTINUE 

PRINT  scaled  transfer  function  and  oc  gain 

WRITeiUP»200#) 

WR  ITe ( UP# 2005 ) (P(I)#I;1#NP1) 

#RITe<LP#2006)  (P«I># I«NP2,NENDj 
#5  CONTINUE 

DCQAjNiPl J) 

#RITE(LP'2Ci6)DCGAIN 
FIND  ZEPO$ 

CAUL  R80T(pii},NN8RD#zER0S#NZ#A#-LP) 

c* • * . .find  poles 

CALL  R80T(P(MP?) #NDORD#POLES#NP/A#»LP> 

unscale  transfer  function  and  frequency 

TEf'PiFSCALE 
DO  6j  I •2>NPj 
P(I).TemP»R(I) 

TprPiTpMP.FSCALF 
61  continue 

NP3»nP2+1 


TEf-PiFsCAuE 
DO  63  J»Np3#nEnD 
Pl I ) .TEmP#P( I ) 

TENP.TEMP.FSCALE 
63  CONTINUE 

DO  6j>  I»1#N. 
w<I).wiD/FSCaLE 
62  CONTTNuE 

PRINT  UNScaleD  TRANSFER  FUNCTION 
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Table  B-l.  Transfer  Function  Identification  Program  — Program 
to  Identify  S- Domain  Transfer  Function  (Continued) 


244 

0000036? 

245 

0000036* 

266 

00000365 

247 

0000036® 

100 

248 

00000367 

c* • • • • 

269 

00000370 

250 

00000371 

251 

00000372 

258 

00000373 

253 

0O0QO376 

254 

00000375 

255 

00000376 

256 

00000377 

257 

00000400 

258 

00O0O401 

259 

O0O0O402 

260 

00000403 

261 

0000040* 

71 

268 

00000405 

263 

OoOoOao® 

264 

OoOoOaO7 

68 

265 

00000410 

266 

0000O411 

267 

00000412 

70 

268 

00000413 



269 

0900041* 

270 

00000415 

271 

00000416 

278 

0000041 7 

273 

00000480 

274 

00000421 

275 

00000422 

276 

00000483 

277 

00000426 

278 

00000485 

279 

0000°486 

280 

00000427 

281 

00000430 

73 

288 

O0O0O43I 

283 

0000°432 

284 

00000433 

.69 

285 

09000436 

286 

00000435 

287 

00000436 

72 

288 

0o00°437 

C« • • • • 

289 

09000440 

c» • • • • 

290 

OoOcOAAl 

291 

00000442 

' 

298 

00000443 

293 

00000446 

294 

00000445 

295 

O9O3O446 

296 

00000447 

297 

00000450 

298 

00000451 

20 

299 

00000452 

300 

00000453 

301 

00000456 

30 

308 

00000455 

303 

00000456 

40 

304 

OoOoOaS7 

15 

WRITE<LP'20l4> 

WRlTE(LP'8PQ5)<P(lM  IjIaNPI) 
WR1TE<LP'8006)<P<ImI.NP8#NEND> 

CONTINUE 

PRINT  UNSCALEO  ZEROS 

iord;nn»ro 

IF(Nn8rO*QT.») iord.a 

WR I TE ( UP* 201 1 1 NNOro * LORD ( IBRD) 

WRlTg(LP'20l5> 

IPT.J 

D8  To  I»l#Nz 

IF (C8MP)ZER0Sm*ZER8S(  Ii/CMPLX(FSCAlE#0.> 

PREA|_*REAuZfROS(I>  > 

PIHAa»AlMAQ(ZER8S{ I) > 

IF  ( A(jS(PlMAj  )«QT<ltOE«7)Q8  T8  71 
WRITE (lP# 2016) IfPREA^ 
liB  TO  6® 

FREa,CABS(ZER8S(  m 
DAMP,-PREal/FREQ 

WRITE  < LP»2016)  IiPREAL*PlMAG<DAr;P#FREu 
IF(JzT(IPT).NE.!)G0  TO  70 
WRlTE^^On) 

IPT.JPT+I 

CONTINUE 

PRINT  UNSr>LE0  POLES 

iordInooRo 

If  ( Nd8RD»(jT • 4 ) IBRD.4 
wRITE(LP'?C12)N08RD#L8RDi IBRD) 
wRITg  < LP'2015 ! 

IPT-1 

08  72  I ■ 1 # Np 

IF(C9Mp)P0LES<I)«P8LES(I)/CMPLx(FSCALE#0t» 
PREAl»REAl(P8LES( I> ) 

PlMAQ«AlMAG(p8LES( I ) ) 

IF ( AbS( plMAg ) »3T»l »0E»7 ) qO  TB  73 
wR!Te(LP»2016)  XPREAL 
GO  Tp  69 

FREQ.CABSjpBLESt I) ) 

DAMpi-pRE*L/FREO 

WRITE  <LP*2016)  nPR'EAL»PlMAG/OAMp,FREg 
IF ( JqT ( IpT) .NE. I )Q8  TB  7? 

WRITe(LP»2017) 

IPT. IPT*1 
CBNT  J NuE 

PRINT  but  FREQUENCY  response  8F  computed  TRANSFER  FUNCTION 
at  Input  data  points 

IF( •nOT*C0MP)GB  TB  101 
wR!TE(LP'?007) 

D8  15 

IF(6tYPE«EQ.LP8LAr)58  TB  go 
IF(8TYPE‘EQ.IPECT)GB  TO  30 
ARGli2C*'ALBGlO(CApS(GA(  I ) > ) 

ARG2iATAN?(AIMAG(GA( I ) ) / RE *L f G A ( I ) )) »R0 
GO  To  4C 
ARGl.CABS) 3A ( I ) ) 

ARG2i ATAN?( AJMAG(GA{ I ) ) jREAL(GA{  I ) ) )#RD 
GB  TO  AO 
ARG1.REAU(GA( I ) ) 

ARG2'.AIMAG(3*(  n > 

WRITE(LP'P00A)w( I)»ARG1/ARG2 
CONTINUE 
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Table  B-l.  Transfer  Function  Identification  Program  — Program 
to  Identify  S- Domain  Transfer  Function  (Continued) 


305  00030460 

306  O000°46l 

307  00000*62 

308  00000*63 

309  0000046* 

310  00000*65 

311  00030*66 

312  00030*67 

313  00030*70 

3l*  00030*71 

315  00000*72 

316  00030*73 

317  00030*7* 

318  00030*75 

319  00030*76 

320  00030*77 

321  00030500 

322  00030501 

323  00030502 
32*  00030503 

325  0000050* 

326  00000505 

327  00030506 

328  0003050? 

329  00000510 

330  OoOoOSll 

331  00030512 

332  00030513 

333  0003051* 

33*  00000515 

335  0o00°5l6 

336  0003051? 

337  00000520 

338  00030521 

339  00000522 

3*0  00000523 
3*1  0003052* 

3*2  00030525 

3*3  00000526 

3**  0000052? 

3*5  00030530 

3*6  .00000531 
3*7  0000053? 

3*8  00000533 
3*9  0000053* 

350  00000535 

351  00000536 

352  0000053? 

353  000005*0 
35*  000005*1 

355  000005*2 

356  000005*3 

357  000005** 

358  000005*5 

359  000005*6 

>60  OOOOOS*? 
>61  00000550 

>62  00000551 

>63  0000055* 

>6*  00090553 
>65  ' 00000556 


UP  array  Or  FREQUENCY  RESPONSE  Te  BE  PLOTTED# 

deleting.  DEStPFD  poles  And  zeros 

lcl  CONTINUE 

lrr  »nOT»OuTjqO  TO  51 
*HlTp <LP'2oi*> 

ExPO.A^eGiO(t.LOW) 

E INT  vL« FLOAT  j NOEC) /FLOAT ( NPO INTS> 

DO  5p  I«l«NP9lNTS 
ExP8.ExPO;E!nTvL 
FRE3IN,  I > i jo -••EXPO* 

FHZ«FREQIn< J */?#/PI 
jw.cmplx  <6*#  fReqin,  in 
TnlM'.CPOLyF,  jw# ZEROS# NZ#jZT) 

TDENOMiCPOLYF,  JW#P8LES#NP#JPT, 

tfn.dcqaIn.tnum/tdenom 

PART) ( I >*20.*Al8010(CABS|TFN)  ) 

PART?) !>«ATAn2(AIMAG,TFN»#REAL(TFN) >.RD 
C*  • • #•  ADD  TImE  DELAY  B*CX  INTO  DATA 
IF(SeNSE  S-.TTCH  2)  57#58 
17  PART? ( J ) *PART2|  I ) -FREQ IN  1 I ) .TDELAY.57.293 
5g  CONTINUE 

WRITE(LP'2019)FREQIN( I)#FHZ#PART{  [ I )#pART?i  I > 

50  CONTjNuE 
ISC-0 

IF(QGAINX.|„E.O.)tSC.i 

51  CONTINUE 

add  time  delay  back  Into  raw  data 

IDSeNSE  SwITCw  2)  *7#*8 

*7  DO  *8  I • 1 # M 

tenpq;cabs(Resp( m 

TENPp;aTAN2(aHAG(RESP(  In«REAL(RESP(  I)  n.«D 
TEMPp.-[EMpp.w,  I)«TDELAY*57*293 
TEMPp.TEMpp.oR 

RESP(  I J«Tempg«CMPLX(C8S(TEMPP)»SIN(TEMPP) ) 

*8  CONTINUE 

CALL  BODPlT)  FREQ  IN#PART1#P/  RT2#I,i#WL0W#NDEC#NPDEC#LP#1HG#1^P#  ISC# 
•NPOInTS#ABS{QGAINX) > 

READ) IcR»2016)  • 

PAUSE 
GO  TO  200 

1000  FORMaT(1X#I9,Iio#*E10.3) 

1001  FORMAT, 20A4) 

1002  FORMAT,  Al) 

1Q03  FORMAT, Al£Ig, 

2000  FORMAT,  iHmRHDENOMlNfTOR  0RDER.#17X#  112/ 

ZlX#  ENUMERATOR  ORDER.#  19X#  112/ 

A1X#25HEQUATI0N  ERROR  TOLERANCE. #10X#Ei2*5/ 

B1X#2iHREA(.  ERROR  T8LERANCE.#1AX#E12»5/ 

C1X#23HFREqUEnCy  SCALE  FACTOR*# 12x#fl2. 5/  _ . 

D1X#29HB*Se  FREQUENCY  FOR  BODE  PLOJ*#  6X#El 2*5/ 
ElX#25HpOlNTS/DECADE  FOR  BODE  PL0T.#7x# IJ2/ 

F1X#3*HT0Tal  POINTS  PLOTTED  ON  pODg  PLgT.#lx#U2/ 

G1X#25HTERm  REJECTION  CRITERION*# i OX# El 2*5// ) 

2001  FORMAT ( lMg# 3JX#8HRAW  DATA#*0X#16HRECTANGULaR  DATA/) 

2002  FORMAT (lOx#El*,5#2]10#3X« I2#*X#El2.5) 

2003  FORMAT  I lXt3Q» 1H.)#17HMAX  ORDER  REACHED# 2 1 10#E12*6) 

200*  FORMAT ( lMQ#gn(  1H#', #3BhTRunCATEP  TRANSFER  FUNCTION! SCALED,  j 

2005  F0RMAT|iH0#1iHNUMERAT0R-.#2X*(6Ei5#5'/, 

2006  FP^MaT(I.Hq#13HDENOMINATOR— #(6El5.5)) 

2007  FORMAT|//1H0420(1H*)#2*HFINAL  FREQUENCY  RESPONSE/) 

2001  FORMAT,  IX,  El2*5#10X#2(ix#E12#5)) 
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378 

379 

380 
3*1 

382 

383 


Table  B-l 


00000555 

00000558 

00000557 

00000560 

00000561 

0000056? 

00000563 

0000056* 

00000565 

00000566 

00000567 

000C0570 

00000571 


0O000572 

0O000573 

0000057* 


00000575 

0o0o0576 


. Transfer  Function  Identification  Program  --  Program 
to  Identify  S- Domain  Transfer  Function  (Concluded) 

2099  20 i ih*  ) #3*hcbmputed  transfer  funct i on < scaled)) 

2010  FORMaTiIHq^oJIH^j/SHDC  GAlN»*El2«5) 

2011  FORMAT (lHjjjodH#) *9HRBBTS  8F  *!2*A2*i7H  BRDER  NUME*ATBRi) 

2012  .FORMAT) 1Hq*2q ( JH* ) *9HRB6TS  BF  *I2*A2*i9H  BROER  DENOMINATOR* ) 

2013  FBRMAT(lH3i89x»17HITERAT]BN  HISTORY/ 

• 1AX/2HEN/ 15X# 1HK«9X* 1HN*  t AX* 2HEC/ ) 

201*  F0RMaT<1Hq,?0(iH*),33HFInAL  TRANSFER  FUNCTIBn(UNSCALED) ) 

2015  FBRMAT(lH0,13x,*HREALi8Xi9HIMAG!NARV#3X*7HDAMPlNa#BX#»HFREoUENCY/) 

2016  FORMAT (5X,|g,5x/A(E12»5)l 

2017  FORMAT (lH;,7xi3H0EL> 

2018  FORMAT ( lH{j  55X/21HPL0TTED  BUTPUT  VALUES/ 

A9x<13HfREq( RAO/SEC  )>10X<8HFREfl( HZ )#l0X»8HGAlN(DB)#i0X#10HPHASE(DE0 
B)/) 

2013  FBRMAT(4<iox*Ei2.5)l 
300c  FBRMaT(2)(<?2h1mPUT  TIhE  DELAY  (SEC)) 

3001  F0RMaT(E15.a> 

300?  F0RMaT<//,I0xm5HTIM£  DELAY  AT  *El0*A*9H  SECONDS*//) 

END 
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Table  B-2.  Transfer  Function  Identification  Program  — 
Subroutine  Bode  Plot 


1 

OoOoOooO 

SUBROUTINE  BODpL7 ( FRE^IN,  QAIMN/PhASIN#  ICA3EzMAxCASE»HL6/./nDFC» 

? 

OoOoCoOl 

• NPC'EC/IM'LGAtN/LPHASzNSCALzNPOINT/SCuAlN', 

3 

00000002 

COMMON  H<50>  jDEI.*<50>iRESP(50>»PHl(50t>A(i,T»%?) 

6 

00000003 

COMMON  S(63>#T|61) . U( 63 > . L*H9D A( 43 > 

5 

OoooOoo* 

COMMjN/kLuGE/H 

KLUGE 

6 

Oo0a0oo5 

dimension  1 9uf ( 120 i 

7 

ooooOoo6 

DIMENSION  rREOTS(l)»aAjN!N(MAxCASE/l),PHAS!M|MAXCASE»n#  I*BRk(11) 

/ 

8 

OOOOOOO7 

♦PRJNTS(ll) 

9 

00000010 

REAL  l£NG 

10 

OoOoOoil 

COmPlEx  flESP 

KLUGE 

11 

OOOOOOI? 

REAL  LAHBoa 

12 

00000013 

DATA  taL*N<,/*H  /»  ISLASH/1N1/#  I t/6Hl  / , I Ml NUS/4H.  / 

13 

OoOOOoi* 

DATA  LSTAR/irt./<LPLUS/lH*/ 

KLUGE 

16 

00000015 

c 

PLpTS  GAIN  AN0  PHASE  (IN  THE  HORIZONTAL  DIRECTION) 

PLOT  33 

15 

0o0o0ai6 

c 

VS.  FREQUENCY  (In  THE  VERTICAL  DIRECTION)  • 

PLOT  60 

16 

OoOoOOl7 

RD.6>5»/ATan(  1 .c) 

KLUGE 

17 

0000°020 

Tn8PJ»8«*ATAN(10) 

It 

00000021 

LEAG.NpDEc 

19 

00000022 

c 

NSCAL.O  SlQNlrlES  THAT  the  DEFAULT  Input  vALUESI-SCGAtNz*SCGAlN) 

20 

0000°023 

C- 

ARE  USED  F8R  The  GAIN  SCALE/ 

PLOT  60 

21 

0000002'* 

c 

nscal  . T SIGNIFIES  that  the  gain  scale  IS  ADJUSTED  AUTOMATICALLY. 

PLOT  70 

22 

00000025 

IlENq»nDEc»npDEC*1 

23 

00000026 

CGAInx'.SCqAIn 

26 

00000027 

IF (NSCAL*E3.0)G8TB*00 

25 

0o00°030 

c 

SEARCH  FOR  THf  HAXIMUM  GAIN. 

PLOT  290 

26 

00000031 

3MAX‘.-1.0e*35 

27 

00000032 

CMIn;i.oE+35 

28 

00030033 

08  360  1 PAINT  . 1,  NPSlNT 

PLOT  320 

29 

O000003* 

Q . GAINlN  I ICASE. IP8INT) 

PLOT  330 

30 

00000035 

IF  UMAX  .LT.  0)  QMAx  . 3 

PLOT  360 

31 

00000036 

IF  (OMIN  .3T . 0)  QMIn  . Q 

PLOT  350 

33 

00000337 

36-  CBnTinjE 

PLOT  360 

33 

00000060 

AdSMIN  , A5S  ( QM IN ) 

PLOT  370 

36 

00030061 

IF  UMAX  «LT.  A8SMIN)  3MAX  • ABSMIN 

PLOT  360 

35 

OoOoOo*2 

CALL  SC*P  (QMAx#  3GAINX) 

PLOT  390 

36 

00000063 

60'.  mqaKx  a IF  1 x (QGAINx  ♦ 0.5) 

PLOT  6C0 

37 

00000066 

3DIFF  . OQAJNX  / 6.0 

PLOT  610 

38 

00000065 

QBOUnO  • 1 . J ?5  • 3G A I NX 

PLOT  620 

39 

O0O0O066 

c 

PRINT  HEADING, 

PLOT  500 

60 

00000067 

hRUe  (I6,  5?0i  LQAIN/  LPHAS 

PLOT  510 

61 

00000050 

52r  FBRMaT  (1H1/37X/15H  PLOT  8F  GAIN  (zAi,l3*>  And  PHASE  (/ Aj, 

PLOT  520 

62 

00030051 

1 1SH)  vs.  FREQUENCY/) 

PLOT  530 

63 

00000052 

c 

PRIM  the  QAls.  Sg*LE. 

PLOT  560 

66 

00000053 

D8  560  I i l#  9 

PLOT  550 

*5 

00000056 

56c  PRINTS  (I)  i - QGAINx  ♦ 3DIFF  • F^OAT  (1*1) 

PLOT  560 

66 

00000055 

6RITE  (I*/  580)  < PRINTS! I )/  I . 1/  9) 

PLOT  570 

67 

OoOflOoS6 

580  F0RMaT<1X,7MRAD/SEC/6X/2HHZ/2Xz6H  GA I N! /9( F6. 1 / 6X )) 

68 

Oo00°057 

XPT.i 

KLUGE 

69 

00000060 

ipeisT  * i 

PLOT  590 

50 

00000061 

ILINE  • 0 

PLOT  600 

51 

00000062 

61c  IlINE  ■ IlIne  ♦ 1 

PLOT  610 

52 

00000063 

C 

RESET  IBuR  TB  *LL  BLANKS. 

PLOT  620 

53 

0003006* 

D8  6*0  I i 1 » 120 

PLOT  630 

56 

00000065 

66C  IBUF  ( I ) • I BLANK 

PLOT  660 

55 

Oq00°O66 

C 

print  horizontal  axes. 

PLOT  650 

56 

00000067 

IF  (<Il!nE..nE.  1)  .and.  (ILINE  .NE.IlENg'))  go  TO  750 

57 

0000°070 

D8  6x0  I * POi  110 

PLOT  670 

58 

00000371 

680  IBUF  ( I ) a I w I NUS 

PLOT  680 

59 

00000072 

C 

pRjnt  tick. MArxs. 

PLOT  690 

60 

00000073 

IF  ( ILINE. NE.l)  GO  TO  69o 

61 

00000076 

IC.I5 

62 

00000075 

ID.IO 

63 

00000076 

69c  IF  (ILINE. NE.ILENG)  QO  TO  700 

66 

00030077 

IC-21 

65 

O0OOO10O 

ID. 9 

66 

O0O0O10I 

70o  DO  730  I • T/  9 

67 

00000102 

73c  I8UF  ( 1C  ♦ I • ID)  ■ II 

PLOT  730 

68 

0o0fl0l03 

C 

PRINT  THE  VERTICAL  AXES. 

PLOT  760 

69 

OoOoOlO* 

750  IBUF  1 20)  • I T 

70 

00000105 

IBUF ( 1 10 ) . I I 

71 

OoOoOlO6 

IBUF  (65)  ^ ISLASH 

PLOT  760 

72 

OoOoOlO7 

C 

calculate  the  log*RJ™m  of  the  plotted  frequenct. 

PLOT  770 
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Table  B-2.  Transfer  Function  Identification  Program  — 
Subroutine  Bode  Plot  (Concluded) 


73 

OOOoOllO 

F«EQLO»AL9Gl6<Hl.en|*|FLeATi:'  InE-ih/lENq 

7* 

OoOoOlll 

c AVERAGE  THE  3*1  s AND  PHASES  OF  THE  POINTS  TB  rE  PLOTTED 

PLOT  793 

75 

00000112 

C PM  THIS  JNE  l nE. 

PLOT  800 

76 

00000113 

COUNT  ■ 0.0 

PLOT  813 

77 

OoOoOll* 

GAIN  ■ 0*0 

PLOT  823 

78 

OoOoOl l5 

PHAS  • • 18O.O 

PLOT  833 

79 

0000011* 

GAINSM  • 0*0 

PLOT  8*3 

80 

0oco0ll7 

PHASSH  ■ o»0 

PLOT  853 

Si 

00000120 

860  IF  ( IP8tNT  .gT.  NP51NT)  GO  TO  930 

82 

00000121 

IF  ( ALOG1o< FREolN( IPO  I NT  1 ) .GT.  1 FREuLG* .5/LENQ) ) GO  TO  930 

PLOT  873 

83 

0000012? 

GAINSM  < GAINSM  ♦ GAIN1N  { ICASEi IPOINT) 

PLOT  883 

8* 

00000123 

PHASSM  ■ PM*sS«  ♦ PHASIN  (ICASEi IPO I NT » 

PLOT  893 

85 

0000012* 

count  . count  * i*o 

PLOT  903 

86 

00000125 

IPOInT  ■ IPOINT  * 1 

PLOT  913 

87 

0000012* 

GO  To  860 

PLOT  923 

88 

00000127 

93 r.  IF  (COUNT  .Ed*  0*0)  GO  TO  1150 

PLOT  933 

89 

OoCoOlSO 

GAIN  • GAINS*  , COUNT 

PLOT  9*3 

90 

00000131 

PHAS  ■ PH*SS*  / COUNT 

PLOT  953 

91 

00000132 

C TE6T  WHETHER  the  GAIN  FALLS  .ITHlN  THE  BOUNDS  i**5i  **5). 

PLOT  963 

92 

09000133 

IF  (GAIN  «3E > eOSOUND))  GO  TB  10*C 

PL»T  973 

93 

0000013* 

IC  ■ 13 

PLOT  983 

9* 

00000135 

990  ENC0DE|6M009i!W8RK(l,|5A!N 

95 

00000138 

1000  F0PM*T(F6.j, 

96 

00000137 

DEC00E<6Air'2i#!H8R*(  1 , 1 laUF ( IC) 

97 

000001*0 

1020  F0RM*Tt6*l', 

98 

000001*1 

GO  To  1100 

PL9T1333 

99 

OoOoOl *2 

1C*C  IF  (GAIN  .LE.  OBOJND)  GO  TO  1 8C 

PLPT13*D 

100 

000001*3 

IC  • UC 

PhOT1353 

101 

OoOoOl** 

GO  T»  990 

PL9T1363 

102 

0o0o0l*5 

C PLOT  jhE  gain. 

PLOT137J 

103 

OoOoOl *8 

10*0  IBUFc  IFIX(!,5.*GAIN/OBOUNO*65.5)).LCAIn 

10* 

0oco0l*7 

C TEST  whether  the  phase  FALLS  wITHIf,  T«  BOUNDS  (.360#  01  • 

PL0T1393 

105 

OO000J50 

110r  IF  (PHAS  ,GE.  (.361.1)30  TO  ,1110 

106 

00000151 

PHAS.PHAS^tah.O 

107 

00000152 

GO  To  ilOo 

108 

00000153 

1110  IFiPhAS.Le. 0.0)30  TO  n?-> 

109 

0000015* 

PHAS.PHAS-?60.f' 

Ho 

00030155 

GO  To  111(1 

Hi 

00000156 

c PueT  the  phase. 

PLOT  1123 

112 

09030157 

112-  IBUF  ( I F 1 x <".?5  . PHAS  ♦ 110.51)  ■ LPHAS 

PL0T1133 

113 

00000160 

C PLPT  data  POINTS 

11* 

00000161 

IF ( APT .GT .m | gO  TO  1122 

kLUGE 

US 

00000162 

1123  IF(AaS(ALOGlr(w(KPT)  ).FRE0LG)*GT.(«5/lEVj)  )3B  TO  112? 

.kLUGE 

116 

OoOoOl 63 

XGA 1n«20 •• ALOGl 0 < CABS (RESP (KPT) ) ) 

KLUGE 

117 

0903016* 

IBUF ( IF  I X ( *5 . »xGA I N/QB0UND*65.b ) ) .LSTAR 

KLUGE 

118 

00000165 

XPHAs.ATAn?( aI*AG(RESP(KPT)  )#REAL(RESP(<PT) ) ).RD 

KLUGE 

119 

O0O0OI68 

1112  IF ( XPHAS *gE . ( *360. ) ) GO  TO  1113 

KLUGE 

120 

OQC00167 

XPHAS*xPHaS*169. 

KLUGE 

121 

0000°170 

GO  To  Hi? 

kLUGE 

122 

O0C0O1 71 

1113  IF ( XpHaS «LF .1* ) 30  T9  1121 

kLUGE 

123 

03030172 

xPHAS.xPHAS.36". 

kLuGE 

12* 

OoOoOl 7 3 

GO  To  11 13 

KLUGE 

125 

0903017* 

1121  IOl’F ( I F I X ( .?s*»PHAS*ll  '.j)  )»LPUUS 

126 

Ooo:0i75 

<PT.<PT*1 

127 

OoOjOl 76 

GO  To  1123 

KLUGE 

128 

03030177 

1122  CONTINUE 

KLUGE 

129 

00000?00 

C PRJNT  THf  (,INE. 

PL0T1 1*3 

130 

0o0o020l 

115'  ohga. 1 - . • .frequs 

131 

0000020? 

OMGAHEiOMgA/THBPI 

132 

0000020? 

WRITES  IH#116i)8MGA,9M3AH7#  ( I3UF(  I )#  I.PO.l  1 A) 

133 

0000020* 

1160  FORMAT  ( 1 X, F7 • 2#  lX#F7«2»3x# 99*1 ) 

n* 

Ooco0?o5 

IF  (ILINE  .lt.ILENG)  30  TO  61 

135 

0900020* 

C PRINT  The  pw*sE  SCALE. 

PLOT  1183 

136 

■00000207 

1199  00  l?0o  I . 1.  11 

PL0T1193 

137 

00000210 

l20r  IwOR<  ( I ) i - 396  ♦ 36  • I . _ ■_ 

PL0T1203 

138 

OOC0O21 1 

HR I Te  (IN.  1220)  (IwORKiJI#  I • 1#  11) 

PL0T1213 

139 

0o0o0?l2 

1 22  c FOPMaT  <7x,7M  PHASE  S#3X#11(I*#5XI) 

PL0T1223 

1*0 

0o000?l3 

RETURN 

PL0T1233 

1*1 

0000021* 

END 

PL0T12*3 
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Table  B-3.  Transfer  Function  Identification  Program  — 
Subroutine  Card  Read 


1 OOOOOOOO 

2 O0O0O00I 

3 OoODOOO? 

6 0o0o0003 

5 OOOOOOO* 

6 00000005 

7 00000006 

8 O0O3O0O7 

9 OOODOOIO 

13  00030011 

11  00030012 

12  00030013 

13  OOCOOOI* 
n Oooo0oi5 

15  OoOoOoi* 

16  00030017 

17  00000320 

18  00000021 

19  00030022 

20  00030023 

21  0000002* 

22  00030325 

23  00000326 

26  00000327 

25  00030030 

26  00030031 

27  00000032 

28  00000033 

29  0000003* 

30  00030035 

31  00000036 

32  00000037 

33  00000040 

36  00000061 

35  00000362 

36  00000063 

37  0003006* 

38  00000065 

39  00000066 

60  00000367 

61  00000050 

62  00000051 

63  00000052 

66  00000053 

65  0000005* 

66  00000055 

67  00000056 

68  00000057 

69  00000060 

50  00030061 


SUBROUTINE  CA«D«D(UUl*LUBiM/MHAx) 

C«»*.  •SUBROUTINE  T8  READ  FREQUENCY  RESPONSE  DATA  8fF  CARDS 

and  the  covert  data  to  internal  typeIreciaNqulaR) 

COMMON  w(50) #DELW<50>#RESP{50>*PHI(50WAf  61*62) 

COMMON  St63)«T(63)iU(63)#LAMB0A(63) 

DIMENSION  vAlUE ( 3 ) / L ABEL ( 3 ) /L I MAqE ( 20 ) 

COMPLEX  rESP 
REAL  LAMBDA 

DATA  LB9DE/iHB/*LP0LAR/lHP//LRECT/lHR/iLEND/iH»/#LH2/lHH/ 
DATA  LbU</1H  /iLDEQ/lHD/#LRAD/iHR//LDB/lHD/»LMAQ/iMM/ 
PI«6.*aTAn( 1.0) 

DR»P 1/6 • /65 • 

M.O  . 

1 h«M*j 

IE(M.f,T.MHA^jr,0  to  2 
READ( LUl AjoOQ ILIMAqE  _ 

DECODE (1*1001*LIMAqE<1) >l  in 
IF(LIN.EQ.LEN0)RETURN 

call  STRlNGiLlMAG£'VALUE*LA6EL) 
w<M);vALUgai  . . 

IF(LABEU<  n .EQ*LHZ)W<M>»6(M)»2..PI 
IF(LIN.EO.lPOLaR)gO  to  10 
!F(LJN,EQ.LRECT)gO  TO  20 

BMAG.VALUEI2)  ■ 

IF(LABEL(2) .EQ.LBLK.OR.LAbEL(2).EQ.LDB)BMAQ'1O»**(BMAg/20*) 
BANG* VALUE ( 3 1 . ‘ 

IF ( LABEL ( 3 ) .EQ.LDEg*0R.LABEL(3).EQ.LBLK)BANG.BANq*DR 
RESPi  M)  •bmaq#cmplx(  cbs  (Bang)  *s ini  bang)  > 

GO  To  30 
10  CONTJNUE 

PMAG.VALUft?) 

if  <label(  2)  .eq*ldb)Pmag*io»«»(pmag/2o.  j 
PANGiVALUE(3) 

IF ( LABEL  < 3 ) .EQ*LDEg.0R.LABEL(3).EQ.LbLK)PANQ.PANq»DR 
RESpj  M ) ■pMAQ«CHpLX  (COS(  pANG)  'S  IN{  pANQ ) ) 

GO  To  30 

20  RESP(M).CmpLX(vALUC(2)#VALUE(3)) 

30  iriLu9.LT.oiG9  T9  1 

WR!TE(LU9«2000)LlN»(LABEL(!>*vALUE(I>«Ui»3>«W(M>«RESfHM) 

go  to  1 

2 IF<LTN.EQ.LEnD)RETURN 
READ<LUl*ioOQ*LlMAoE 

DEC0dE<1*i001*LJ*a0E(1>  M-JN 

IP  ( L IN.NE,  LEND)  WRITEc  LUOi  200  DLIMAgE 

GO  To  2 

1000  FORMAT ( 20A6 )' 

1001  FORMAT(Al) 

2000  FORMAT  < IX,  Aj,3x*3(ix«Al/lj(*E15>7>*5X*lH«jSx*3(l)(lClS*7)) 

2001  FORMAT (IX* IihExTRA  CARD.I20A6) 

END 
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Table  B-4.  Transfer  Function  Identification  Program  -- 
Subroutine  Change  Signs 


1 OOOOOOOO 

2 OoOOOOOl 

3 00000002 

* 00000003 

5 0003000* 

6 ooooOoo® 

7 OOOOOOOA 

8 0103000? 

9 00000310 

10  Ooooooil 


SUBROUTINE  CMSiARR/N*  1»*AT> 

DIMENSION  aRR(J) 

SUBROUTINE  T8 ’CHANGE  SIGNS  IN  PROPER  PATTERNS 

IP(  IpAT.Eq.  ^return 
!BEGilPAT/4>M 

lNCRi2-H0()(  JpATiP) 

DO  1 I«IBEG,N«INCR 
1 ARR{  1 )i«A(?R(  j) 

RETURN 

END 


Table  B-5.  Transfer  Function  Identification  Program  -- 
Subroutine  Evaluate  Complex  Polynomial 


1 

00000000 

complex  Function  cpoly<coef#n*s) 

2 

OoOOOOOl 

DIMENSION  cOEF<n 

3 

OoOOOOO? 

COMPLEX  S 

% 

OoOoOoo3 

evaluate  a complex  polynomial  with  real  coefficient* 

3 

OOOOOOO* 

USING  hORneRS  rule* 

0 

ooooboos 

CPOLY«CMP|lxiC»EP(NMO.» 

7 

OOOOOOO* 

ip(N,eq»1) return 

• 

OOOOOOO? 

nmi»n-j 

9 

ooooOoio 

DO  1 Ul*NMi 

10 

Ooooooil 

IR.NmI.IM 

11 

OOOOOOI? 

CP0Ly.CP®LV*S*CMPLX(C0EFcIR)/0O 

12 

OoOoOoi3 

1 CONTINUE 

13 

OOOOOOI* 

return 

1 A 

OoOoOoiB 

ENO 

Table  B-6.  Transfer  Function  Identification  Program  -- 
Subroutine  Evaluate  Factored  Polynomial 


1 

OoooOooo 

C8MPlEx  FuNCTlAN  CPALYF(jw»R89TSiNfiRD#  JRt'> 

2 

OoOOOOOl 

c** 

<• .EVALUATE  A FACTORED  PBLYnBMJAl 

3 

OoOoOOO? 

DIMENSION  JPT(1) 

A 

OOOOOOO3 

CBMP|_Ex  rAATS<1 )/J.# JWSO.pRBD/FACT 

5 

OoOoOOOA 

PRPDi < 1 * A 0 * » 

6 

O000O005 

IP^I 

7 

OoOoOOO6 

jwS3'.Jw»Jw 

8 

OoOoOOO7 

D8  1 I.l'NARD 

9 

OoCoOOlO 

IF ( I • Eq • JRT { IP))G8  T3  J 

10 

ooooooil  ■ 

IF  ( AgS^  AImAQ  ( RA8TS(  I ) » ) «l_E • 1 *0E»07  ) <3J  TB  1 

11 

OoOoOoi? 

BMEGiCABS(R8AT5( I) > 

12 

OoOoOoi3 

BmE3SQ»0Meg*AMFG 

13 

OoOoOoi A 

zeta;reaL(Rbats(  ih/bmeg 

1A 

OoOoOoiS 

FACT;j*S0/AMFGSQ-Jx#2.»ZETA/8MEG.(  1 ..6») 

15 

O000O016 

G8  Ta  k 

16 

OoOoOoi7 

3 

fact;-jw/raatS( i >♦( 1 .#o> ) 

17 

00000020 

4 

prbd;prad,fact 

18 

OoOoOo?! 

GB  ta  ; 

19 

0000002? 

2 

Ipi IPAl 

20 

00000023 

1 

continue 

21 

0000002* 

cpplyf;ppao 

22 

OoCoOo?^ 

return 

23 

0oco0c.?6 

end 
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Table  B-7 


Transfer  Function  Identification  Program 
Subroutine  Load  A- Matrix  and  C- Vector 


1 oocooooo 

2 O9O0O90I 

3 0oco0302 

4 00000P03 

5 OOOOOCO* 

6 OOOOOOOS 

7 OOOOOOO* 

s QoooOoo7 

9 00000010 

13  03030311 

11  000000l2 

12  OoOoOol 3 

13  0000031* 

14  0o0o00l5 

15  0003001* 

16  03000317 

17  03030020 

18  00030021 

19  O30*'0d22 

20  00000c<?3 

21  0303002* 

22  0303032^ 

23  00030026 

24  O3O3O027 

25  09000330 

26  00030031 

27  00030332 

28  09030333 

29  0903003* 

33  09030335 

31  O9O9O336 

32  00030337 

33  09000340 

34  O000O34I 

35  09030342 

36  03000043 

37  O9O3O3-* 

38  090J0345 

39  O3C0O346 

43  09030047 

41  O303O05O 

4P  09030351 

43  00030052 

44  03030053 

45  OoOoOoS* 

46  00030055 

47  03030056 

48  03030057 

49  09030C60 

50  OO00O061 

51  09030062 

52  03030363 

53  O9O3OO6* 

54  090j0365 

55  00000366 

56  O9O0O367 

57  03000370 

58  00030371 

59  0?000072 

63  09030:73 


SUBR*UT!*E  LOAC(H#kD«NN) 

C8MM9N  «<50)*DFLW<50)aRESP(50)*PHI<5c  >*A<41  *42) 
C8PH9N  S(43) .T{43),U(43)»LAMBDA(43) 

REAL  LAyBD* 

C8*PlEx  rES» 

LBGIcAl  EvE'J#800*L^3*SM*uM»TM 

subroutine  t«  l«ad  a. matrix  and  c-vectsr 

n»aD 

npi«nn*i 
NP?»«gN*2 
MS I Z f ■ nD+nn* { 

!LC8^M.mSize4i 
fSP2.MSlZE*2 
08  1 Iil'MSP? 

riRSTiCftMPUTE  uAMBDA*S*T,U 

ExENiM80( J*2) *E3»0 
UA^BpA  ( ! ) ■ ">. 

S ( I ) aO • 

T(  I ) »0  • 

U( I).0. 

ce  2 j«i*m 

Sa“A3«CABS(RESP( J) ) 

SjmAj.SCMaG.SQmAG 

TemP.Ph! ( J)»DEL*(U)»W(  JU»( 1*1 ) 

IF(EvEN)T(  I ) iTEMP*AlMAG(RESP( J))*T(  I ) 

IF|EvE\)Ga  T«*  ? 

LAMB9A( I ) tLAMBDA ( I ) 4TEMP 
S( I)fS( I)trrMp.REAL(RESP,  j) ) 
i ( n.ui  1 utemp.sqmag 

2 CSMTlNuF 
1 C8NT I NjE 

FILL  MAIN  AND  UPPER  DIAGBNAL  8F  a a*o  THEN  TpANSPBSE 

08  4 I a 1 » “S I 7E 
08  5 Ua I * *S l 7E 

LMh«  I »Lf  •NO’  . A' D.  jtLE*NPi  •AND.mBO)  Uj»2)iE.*' 

Sm»I  .LE*npT.  aND.J.GT.NPI.AND.mODi  I*(  j-NPJ  ) #2)  *EU*1 
TM.I  ,LE«Npi  .ANC*J»GTtNPl.AND.M80(  !* ( j»NPl) * 2) *EQ*0 
I'M.I.GT.Np} . aND.j.gT.Npi  .AND.M9D(  ( I ♦ J I * 2 ) ♦ E)  • C 
*( I A j)»0* 

if(Lmu)A<  1 1 ji'.la^bdai 

IF  C U- 1 A(  I,  JU'Jl  ( I-NPU2)4(  j-NPi-i ) ) 

if(Sv)ai  i,jj;s(i4<  j-npdi 
if{Tv,)ai  Lj);t(  !♦<  j-NPin 
*( j# 1 i.ai I* J) 

t C6N  T j NuE 
u C8k  T 1 NuE 

CHANGE  S 1 GVS  I*  PR8PER  PLACES 

FIRST  paPtjtjSn 

08  6 jil'MPl 

ipat;mpdi j-j . 4 1 ♦ 1 

CALL  ChS(am,J)»NP1*IPAT) 

* Cef-Tl^f 

PARTITION 

D8  7 ji^P?.MSI?E 

• ipAT’.Meoi  j-NPi;2#4>;i 
CALL  ChS( A ( J » J J a NPl * JPAT) 

7 C8NTtNjE 

c.**. .third  paRtitibn 
ce  B jil*NPl 
iPATiMODI  j-1 ,4ul 
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Table  B-7.  Transfer  Function  Identification  Program  — 
Subroutine  Load  A -Matrix  and  C-Vector 
(Concluded) 


61  0002007* 

62  03030075 

63  OjO00076 

6%  00030077 

65  O0O3O10O 

66  O0O3OIOI 

67  00000102 

68  00000103 

69  OOCOOIO* 

70  OOOOOIO® 

71  OOOOOIO* 

72  00000107 

73  OOOOOHO 

7*  00030111 

75  00000112 

76  O0C0O1 1 3 

77  0003011* 

78  OoOoOll5 

79  00000116 

80  00000117 

81  00000120 

82  00030121 


CALL  CHS(A<np2,J>,n,IpAT> 
f COMfNjE 

FOURTH  PARTITION 

DO  9 

iPATiMODI  j-NPi;2#<0*l 
CALL  ChS( A(NP2.j)«»# IpaT) 

9 CONTINUE 

LOAD  C-vEcTOP 

DO  lj  I«l,NPj 
EVEN, MODI I>2)»Ep*0 
ODD*  .NOT  ‘EVEN 

if i odd  > a i i , ilcaum) ist  1 1 

IF(EvEn)A{ I, ILCOLM).T(I) 

Ir  CONTINUE 

DO  lj  I*NP2#mSI2E 
EvFN.MoDI |-NP1,2)»EQ»6 

odd*.not*even 

IF|OoD)A(  1, iLCALMjio. 
1F(EvEn)A(  I,ILC0LM)*U(I.NIP1*1) 
li  continue 
RETURN 
ENC 


Table  B-8.  Transfer  Function  Identification  Program  -- 
Subroutine  Move 


1 

ooooOoco 

SUBROUTINE  -0VFC AIN, abut, N,  IBEGIM 

2 

Oooooool 

DIMEvSjON  a I *J ( 1 ) , ABUT  ( 1 ) 

3 

00000002 

DO  1 I ■ 1 * N 

A 

ooodooo3 

ABUT  ( I » » A I Ni  ( I) 

5 

OoODOoOA 

1 

CO^  JNjE 

6 

OqooOoO5 

c. 

....insert  1*  in  proper  place. 

7 

0000O006 

IENOjN;i 

8 

OaooOoo7 

TEMP.l, 

R 

0000O010 

DB  2 I.IBEGIN/ IEND 

10 

00000011 

TEMPO*aOUT( I ) 

11 

00000012 

ABUT] I ) *tEHP 

12 

00000013 

temp, tempo 

13 

OOOOOOIA 

2 

CONTINUE 

16 

OoOoOolS 

return 

15 

00000016 

END 
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Table  B-9.  Transfer  Function  Identification  Program  — 
Subroutine  Root 


1 

OoOoOooO 

SUpR«UTlNE  ROOTiPaNORDaROBTSiNROOTSaSCRaLUO) 

z 

OoOoOOOl 

. 

COMPLEX  *00T§<i) 

3 

OOOOOOO? 

DIMENSION  P( ! ) ;SCR(NBRD#1 > • 

4 

OCOOOOO3 

if ( n^ro«qt • 1 i go  to  10 

5 

OoOoOOO* 

ROOTS! })«C*PLX<-P< 1 >/P(2)aO») 

6 

OOOOOOOS 

NRPOTS* 1 

7 

OOOOOOC* 

GO  Ti  5 

8 

OoCcOoO7 

10 

N0POP1.NORO*? 

9 

OoOoOolO 

nORO^I^Ord.i 

10 

OoOjOoi 1 

DO  1 lil'NORDMl 

11 

OoOoOoi3 

DO  1 j'.l'NORD 

12 

00000013 

SCR< 

13 

OpOoOSl* 

IF ( J.Eq«I  I + h >SCR<  1#J)«1. 

14 

OoOoOOl5 

1 

CONTjNyE 

15 

OoOoOOl* 

DO  2 M'NSRD 

16 

OoOoOoi7 

SCMnORDa  j)'.-P(  I >/P( N0RDP1 ) 

17 

00000020 

2 

CONTINUE 

IS 

00000021 

IP<NSR0.Le.?)G9  to  11 

19 

0000002? 

CAUL  H£SSen<nORDaSCRaN8RD> 

20 

00000023 

11 

CALL  GRC*LL< N0RD<SCRi ROOTS# NR90TS/N0RD) 

21 

0000002* 

krpots.nRbots/? 

22 

000000?5 

p 

IF<LuO.LT.O)RETURN 

23 

00000026 

kRITE(lUO,20oO> 

24 

O0C0OO?7 

CO  3 1.1'NPOOTS 

25 

00000030 

PREAl«RE*L(RpOTS(  1 ) ) 

26 

OoOoOo3l 

PIMAc,.A1MA3(R09TS(  I ) ) 

27 

0000003? 

IF  ( ArS(  pln*0)  «fjT  «1  «0E«7 ) qB  TO  4 

28 

00000033 

t-RITF(LU0#?00l  )pREAL 

29 

0-000003* 

GO  T-}  3 

30 

00000035 

4 

FRE3iCABS(R0f»TS(  1)) 

31 

O0000036 

DAf-pi-pREAL/FREQ 

32 

00030037 

i»RITe(LU0#?001)PREAL#p1Mag#DAmp#FREQ 

33 

0^0o0o40 

3 

CONTINUE 

34 

00000041 

RETURN 

35 

ooooOo*? 

2000 

FeRMAT(lPp,nX»4HRE4L#8X#9HIMA3lNARY#3X#7HDAMpING#5Xi9HFREuUENCY/> 

36 

00000J43 

2001 

rOP,lAT{l2x*4<Ei2»5)  > 

37 

OOC0004* 

EMC 

Table  B-10.  Transfer  Function  Identification  Program  -- 
Subroutine  Scale 


1 

O0O3O00O 

Subroutine  sc*p  o^ax*  qscalxi 

SCAP 

10 

2 

OoOjOool  ' 

C ESTABLISHES  The  L'PpER  BOUND  OSCALX  OF  THt  SCALE* 

SCAP 

20 

3 

4 

ooooQoo? 

0o0o0oo3 

DI^EnS J8N  SCALxdl) 
DATA  NSC*LX/il/ 

SCAP 

160 

5 

6 

OoOoOoC* 

OoOoOOOS 

IF  <-)NAX  .LE.  1*0)  G6  TS  260 

SCAP 

180 

7 

00030006 

fact  • o*i 

SCAP 

190 

8 

0000°C07 

20o 

fact  • 10, n , FACT 

SCAP 

200 

9 

OOOOOOIO 

D8  240  ISCALX  < 1,  NSCALX 

SCAP 

210 

10 

OoOoOOll 

DSCAlX  ■ FACT  ; SCALX  (ISCALX) 

SCAP 

220 

11 

00000012 

IF  OHAX  <lE • '.SCALX ) GO  T9  363 

SCAP 

230 

12 

00000013 

2*C 

CONTINUE 

SCAP 

240 

13 

OoOlOol* 

Gd  To  geo 

SCAP 

250 

14 

OnosOol? 

26- 

fact  • 1*« 

SCAP 

260 

15 

00000016 

asE«  • i*j 

SCAP 

270 

16 

OO0JO31 7 

28-' 

fact  • CM  • Fact 

SCAP 

26  0 

17 

00000020 

09  340  ISCALX  • If  NSC*LX 

SCAP 

290 

18 

O'.000j?l 

ISCAlY  • NSCALy  ♦ 1 - ISCALX 

SCAP 

300 

19 

0c0o0o22 

2SCALX  • 3'JE* 

SCAP 

310 

20 

oocoO:?3 

ONE"  • FACT  , SC*L»  (ISCALY) 

SCAP 

320 

21 

0oo:0^?4 

IF  ( ;MAX  ,"iT.  •v'lEW)  30  T9  360 

SCAP 

330 

22 

0"'C00025 

34- 

CONTINUE 

SCAP 

3*0 

23 

Gd  To  ?80 

SC*P 

350 

24 

OoooO:?7 

36f 

RETURN 

SCAP 

360 

25 

0?0o003c 

END 

SCAP 

370 
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Table  B- 11.  Transfer  Function  Identification  Program  — 
Subroutine  Logical  Function  Search 


1 

OoOoOooO 

LOGICAL  EunCTION  SEARCH, LABEL# ACCEPT) 

2 

OoOoOod 

INTEGER  ACCEpT<5') 

3 

OOOOOOO? 

do  i r.i's 

A 

OOOOOOO3 

SEARcH'.LAbEL.EQ. ACCEPT,  I) 

5 

OOOOOOO^ 

IF(SEARCH)RETURN 

6 

00000005 

i 

CONTINUE 

7 

OOOOOOO6 

RETURN 

8 

OOOOOOO7 

END 

Table  B-12.  Transfer  Function  Identification  Program  -- 
Subroutine  String 


1 oooooooo 

2 OoOoOoo! 

3 OOOOOOO* 
A OOOOOOO3 

5 OoOaODO* 

6 00000005 

7 0000000$ 

8 OOOOOOO7 

9 OOCOOOIO 

10  OOOOOOU 

11  OoOoOol3 

12  OoCoOol3 

13  OoOoOol4 

1*  Ooooool5 

15  OoOoOo-l6 

16  0000001 7 

17  00000020 
19  OOC0OO21 

19  00000022 

20  00000023 

21  OoOfOo?* 

22  00000025 

23  O0O0O028 
2A  00000027 

25  00000030 

26  00000031 

27  0002033? 

28  00000033 

29  0q0c003* 

30  00000035 

31  00000036 

32  00CJ0O37. 

33  OcOOOoAO 

3A  000000^1 
35  0o0o0?*2 


SUBROUTINE  STR1 no ( L I HAOE#  VALUE# LABEL > . 

DIMENSION  ACCEPT (5) # LABEL ( 3)  #L  IMAGE ( 20 ) 

DIMENSION  VALUE! 3) 

integer  aCcept 
LOGICAL  Pind#SEARCH 

DATAi  ACCEpTj  n#I»l#5)/lH  ,lHD#iHR,lHM,lHH/,LBL</iH  / 
DECBpE;3#ioo6'Ll"AGEm  JLABEUl) 
riND.SEARcM(L*BEL(l>#ACCEPTj 
IF iFIND)0eCODE(17#3O00*L IMAGE! A) ) VALUE!  1) 

IF <F]ND)00  TO  1 
LABEL!  D'LBLK 

DECODE ( 1 7#  30Ql »L I MAGE ! A ) » VALUE ( 1 ) 

1 DEC0pE(2#iC01'Ll*AGE!8> )LAUEL,2) 

FIND. SE*RcHILABEL(2># ACCEPT) 

IF(F  IND)0eC8DE|  16#  3002# L IMAGE ( 8)  )VALUE!2) 

IF(FIND)Go  TO  2 
LABEl(2)»LBU< 

DEC0DE(l6t3oo3#LlMAGEi8,>  »VALUE(2) 

2 DEC00E(1#1OO2»LIMAgE!12))LABEL(3) 

FINDiSEARCHj LABEL! 3)# ACCEPT) 

IFlFlNQjDECflCEf  15#  300A«LIMAGEll2>)  VALUE  (3> 

IF (FjndjReturn 

LA8El<3)»LBL< 

DECODE! lB,3on5#L IMAGE! 12) ) VALUE! 3) 

RETJRN 

1000  F0RMaT<2X, Al  , 

1001  form*T(1x.ai  , 

1002  FORMaT(AI) 

3000  FORMAT ( 3X.  E 1 A »7) 

3001  FORMAT  (2X,Fl 5*7) 

3002  FORMAT, 2X,E1A'7) 

3003  FORMAT, 1X.E15*7) 

300*  FORMAT,  1X#Ha«7) 

3005  FORMAT, El5«7) 

END 
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